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I 


Abstract 

This report is intended to present the total picture of the 
work which has been done by Dr. Farber's group at the University of 
Florida in arriving at credible explosive yield values for liquid 
rocket propellants. 

The results are based upon logical methods which have been well 
worked out theoretically and verified through experimental procedures. 

Three independent methods to predict explosive yield values for 
liquid rocket propellants will be described. All three give the same 
end result even though they utilize different parameters and procedures. 
They are: 

I. The Mathematical Model 

II. The Seven Chart Approach 

III. The Critical Mass Method 

A brief description of the methods, how they were derived, how 
they were applied, and the results which they produced will be given in 
the text. 

The experimental work used to support and verify the above methods 
both in the laboratory and in the field with actually explosive mixtures 
will be presented. 

In the last two parts of this report the methods developed will 
be used and their value demonstrated in analyzing real problems, among 
them the Destruct System of the Saturn V, and the early configurations 
of the Space Shuttle. 



II 


The depth to which the material will be presented here depends 
upon whether and to what extent some of the material has been presented 
in previous publications. 

There exist about two dozen reports and papers which were presented 
by Dr. Farber's group and publications in the Annals of the New York 
Academy of Sciences, various Space Congress Proceedings from both the 
East and the West Coast, Transactions of the International Cryogenic 
Engineering Society, Proceedings of the Department of Defense Explosives 
Safety Board Annual Meetings, The International Fire Institute Bi-Annual 
Reports, etc. 

In addition, some of the material has been presented both formally 
and informally to NASA groups at NASA Headquarters in Washington, D.C. , 
to NASA groups at Kennedy Space Center, to NASA groups at the Marshall 
Space Flight Center, to groups at the Air Force Rocket Propulsion 
Laboratory, at Edwards Air Force Base, to members of the Aerospace 
Corporationq to the Ballistics Research Laboratory at Aberdeen Proving 
Grounds, and others . 

It is believed that this report presents the only rational approach 
in existence today to predict explosive yield values of liquid rocket 
propellants and other characteristics. 

This work, carried out over a period of nine years, is presented 
consistent with the numerous papers and reports produced during this 
period. It is believed that this produced the clearest presentation 
even though it created some cosmetic, notation and numeration problems. 
Quality photocopies can be obtained from the master report residing 
with NASA. 
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Part I 

Methods of Analysis for the Prediction of Explosive Yields from Liquid 
Rocket Propellant Explosions 

With the large scale use of liquid propellants in the space pro- 
grams of this country for the Saturn Rockets, Skylab, projected 
Space Stations, and the Space Shuttle it became imperative to have know- 
ledge of the characteristics of these liquid rocket propellants . 

The yield from liquid propellant explosions as a result of missile 
failures is of extreme importance in assessing the hazards to astronauts, 
launch- support personnel, launch- support facilities and surrounding com- 
munities. Since explosive tests of large liquid propellant rockets are 
not practical because of the costs and hazards involved, prediction 
methods must be used in estimating the expected explosive yields. 

For this reason three independent methods have been developed by 
Dr. Farber and his group at the University of Florida, and all three 
give essentially the same predictions even though they use different and 
independent parameters to do this. 

The methods are briefly described below with the detail left to 
other publications dealing specifically with each one of the methods. 
Only the results from these methods pertinent here will be mentioned. 
Many c idler results have been obtained describing the behavior of the 
liquid rocket propellants 

A. The Mathematical Model 

1 A OA AO 

A mathematical model 1 was developed by Dr. Farber about 



the year 1964, with very limited information available, to describe 
some of the characteristics of liquid rocket propellants. Because of 
the limited data and infoimation at hand at that time, it was difficult 
to evaluate the validity of this model. But since the model was designed 
to be quite flexible, it was assumed that new information could easily 
be incorporated into the model. 

Considerable information has become available since the development 
of the model, and it has only supported the original work. Data 
obtained by the principal investigator and this group of the Uni\ ' y 
of Florida by instrumenting two 25,000 lb LOX/RP explosive experiments'*, 
and one 200 lb. LOX/RP cold flow and explosive experiment carried 
out at the Air Force Rocket Propulsion Laboratory at Edwards Air Force 
Base, California, which established the Yield Function - Spill Function 
relationship. The data reported by Project PYRQ as soon as it became 
available was also used to check the results predicted by the model. 

This information, including the inert mixing experiments carried out 
in our laboratories, increased the confidence in this model and estab- 
lished it as a very useful tool to describe the liquid rocket propellant 
behavior. 

The basic idea for the model is that in the case of an explosion of 
liquid rocket propellants only the portions which are actually mixed 
(fuel and oxidizer) can take part in the explosion. Furthermore, since 
this must be the case, a definite functional relationship must exist be- 
tween the mixing of the fuel and oxidizer, the Mixing Function, and the 
expected explosive yield, the Yield Function. 
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To make it easier for the reader, a brief description of the model 
will follow so the references wi 1 ! only be needed if greater detail and 
insight is needed. 


Yield Function 

The Yi^ld Function, as used in this work, is defined as the fraction 
of the theoretical maximum yield which is actually obtained from: 


7 


Y 

Y Theor. Max. 


CD 


This yield can be expressed as TNT equivalent yield on an energy, 
or other basis, but care must be exercised in predicting damage, since the 
pressure-time trace for liquid propellants is considerably different 
from that of TNT, expecially in the near field. 

Mixing Function 

The Mixing Functic. as used in this report, is defined as the 
fraction of the total volume of propellants mixed at any time x , 
multiplied by some modifying factors. 

x ■ Wf < 2) 

x •» Mixing Function 

■ Propellant Volume Mixed 
Vp • Total Propellant Volume 



4 


Ep = Turbulence Factor 
Fg = Boiling Factor 
Fp = Freezing Factor 


Inert laboratory experiments utilizing such fluids as water and oil, 
hot wax and water, hot oil and water, and water, and kerosene, 
etc. established the factors F^, Fg, and Fp. 

It was shown that in the early stages of mixing these factors have 
a value near one, and thus the Mixing Function is essentially the normal- 
ized nixing volume. This latter fart was also established by the 25,000 
LO^/KP and the 200 lb. LO^/RP explosive experiments. 


Yield Function - Mixing Function Relationship 

In the development of the mathematical model it was assumed that the 
relationship between y and x can be expressed as 


y - 



(3) 


where b, c, and d are constants. This expression containing three con- 
stants is very flexible. 

Again the explosive experiments of 25,000 lb. LO^RP and the 200 lb. 
LC^/RP cold flow and explosive experiments proved the above relationship 
valid. 


The Mathematical Model 


With the relationship between Mie ' 


on (y) and the Mixing 
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Function (x) established, the mathematical model can be formilated, 
resulting in a statistical function which is capable of incorporating 
the above y-x relationship, and is able to provide for the valid esti- 
mating procedures of the parameters involved. 

The statistical function is a modified Dirichlet bivariate surface 
with four parameters a, b, c, and d. It is 

c . , dT fa ♦ b + cl d-1 ,, d^a-1 b-1, d .c-l ... 

ffx.l') * KiffOi) f(<j x (1 - X ) y (x -y) (4) 


where r is the Gamma Function. The only restrictions on this function 
are that: 

y > 0, x > 0, y < x d , d f 0 


To fully define the above function, it is necessary to evaluate 
the parameters a, b, c, and d on the basis of the particular y-x rela- 
tionship describing the physical phenomena. This can be done by following 
statistical estimating procedure. 


Defining 



( 5 ) 


four simultaneous equations can be written for the four parameters 
a, b, c, and d. 
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In v • ¥ (b) - ¥ (b + c) (6a) 
In v = In (b) - In (b + c) (6b) 
In u = ¥ (a) - ¥ (a + b + c) (6c) 
In u * In (a) - In (a + b + c) (6d) 


where a bar over an expression indicates the average value of all 
available values 

In indicates the natural logarithm (base e) 

¥ is Euler's Digamna Function 

From the mathematical model, the modified Dirichlet bivariate sur 
face, a wealth of information can be extracted. Some of these are: 


A. 


Probability Distribution of the Yield , 


PyOO “ 



f(x,y) dx 


( 7 ) 


From this probability distribution, the average yield value can be 
found as well as the confidence limits, indicating that a certain per- 
centage of all yield values lies below the selected ,r ield value. 


B. Probability Distribution for the Mixing Function, P„ 



P x (x) - J f(x,y) dy (8) 


This distribution can by analyzed the same way as in part (A) above. 
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The function for a is a distorted S curve as seen in the refer- 

1 4 24 42 104 

ences * * * * , which also show how this function is obtained from 

the best available information. 

The only experiments, except the ones in the laboratory, which were 
fully instrumented to obtain the Yield-Spill Relationship were the two 
25,000 lb. LO^RP explosion tests and che 200 lb. LC^/RP cold flow and 
explosion experiment. The Yield -Mixing relationship v.hich was verified 
in these experiments was assumed to hold true in the remainder of the 
experiments and failures. 

The effect of parameter a on the yield is rather small in the range 
of today's liquid propellant rockets, and thus the predictions can be 
considered very good since large variations in the parameter produce 
only small changes in the predicted explosive yields. 

Figure 1 presents the results from the mathematical model. The main 
curve presents the average explosive yield values as a function of the 
propellant weight involved. Also shown is the 95 percent confidence 
limit. For the Saturn V, the prediction is that in the case of a catas- 
trophic failure the average yield would be 3.5 percent and that in 95 
times out of 100 such cases it would be less than 9 percent. 

The mathematical model, as used here, included all types of pro- 
pellants as well as all modes of failures. 

It is clear that the mode of failure, as well as the propellant type, 

2 3 15 

has a distinct influence upon the actual yield obtained ’ ' . If only 

a particular type of failure or a particular combination of propellants 
is to be investigated, then only that information should be used for 
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Figure 1 Estimated Explosive Yield as a Function of Propellant 
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analysis in the mathematical model, and the values of the parameters 
a, b, c, and d will change. The average explosive yield value should 
be more accurately determined in such cases and the confidence limits 
will be found, in general, closer to the average values. 

Thus for large-scale liquid propellant rockets it may be desirable 
to control the mode of failure with a properly designed destruct system 
to give a minimum explosive yield. With such a procedure, the explosive 
yield value can be lowered and the explosion yield prediction accuracy 
increased. 

B. The Seven Chart Approach 

An approach predicting the most probable yield, the most probable 
mixing, probability distributions, confidence regions, confidence limits, 
and so forth, by means of a mathematical model was presented in section A. 

The method described there accomplished the ultimate goal of leading 
to a valid prediction procedure of yield, mixing, etc., of liquid pro- 
pellant explosions; it did not provide an insight into the physical 
phenomenon occuring. 

The approach in this section suggests a more fundamental method to 
this problem by considering the physical phenomena in detail. This 
approach, can, through understanding of the physical processes and 
phenomena, provide the information necessary to control these processes. 

The method presented here and earlier^’^’^ is referred to by the 
authors as the "Seven Chart Approach" since the procedure can be sum- 
marized in seven charts, constituting a complete, well -planned program. 
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outlining the necessary steps to be followed. 

Furthermore, the "Seven Chart Approach" uses presently available 
information regarding these poorly understood phenomena producing the 
liquid propellant explosion yield. It points out where more theore- 
tical and experimental work is needed, and what information it should 
provide. In this manner the ideal balance is obtained between theory 
guiding the experimental work and the results from the experiments 
modifying the theory. For these reasons this procedure is able to reach 
the desired goals along a more direct route in the shortest possible 
time and at minimum cost. 

Previous theoretical and experimental investigations suggest that 
the actual phenomena producing the yield in liquid propellant explo- 
sions can be divided into groups which lend themselves to separate 
study. These may be both theoretical and small-scale experimental in- 
vestigations . 

In this method the problem is divided into three groups of phe- 
nomena which can be studied separately but when combined allow the 
desired prediction. The three groups are: 

I. The Yield Potential Function 

II. The Mixing Function 

III . Delay and Detonation Times 

These allow the incorporation of the basic characteristics of the 
particular propellants involved, of the missile design configuration, 
and of the mode of failure. 

The Yield Potential Function (I) is basically controlled by 
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chemical kinetics, the Mixing Function (II) by the principles of 
hydrodynamics modified by heat transfer, and the Delay and Detonation 
Times (III) by characteristic functions for some propellants such as 
hypergolics or by random processes for others. 

The separate studies can be combined by taking the Yield Potential, 
when expressed as a time function, and multiplying it by the Mixing 
Function to obtain the expected yield at any time after the start of 
the failure or after the mixing has begun. This Mixing Rinction will 
be different for different modes of failure and missile configura- 


tions 


15,75 


The actual expected yield can be determined by superimposing the 
Delay and Detonation Times upon the above obtained expected Yield 
Function, either as a fixed value where applicable or as a statistically 
most probable value with proper confidence limits. These Delay and 
Detonation Times are characteristics of the propellants such as hyper- 
golics or cryogenics, modified by the propellant quantities, missile 
configurations, modes of failures, and so forth. 

The total procedure can be summarised, with the seven 
charts^ 104 supplying the necessary Liformation, as the rela- 
tionship 


where 


y - (y p - *) T * 

y expected yield 

y Yield Potential at time t A 

P 

x Mixir Function at time t* 

t* Most probable detonation time 


( 10 ) 



13 


For the development of the seven chart method, conditions were 
assumed so that quantitative results could be calculated for cases 
which were investigated experimentally and for which results are 
reported in the literature. This gives more meaning to the procedure 
suggested and allows comparison of results obtained by the "Seven 
Chart Approach" with actual test results. The approach would be the 
same if other initial conditions, propellants or configurations were 
used. 


The Yield Potential Function 

The Yield Potential Function for any combination of propellants 
as a function of time can be obtained from theoretical considerations 
in four steps as follows: 

Maximum Theoretical Energy Release (Chart 1) 

The maximum amount of energy which can be released from any par- 
ticular liquid propellant fuel -oxidizer mixture can be calculated em- 
ploying the basic laws of chemical kinetics. 

Figures 2A and 2B (in greater detail) show the results from such 
calculations for a three -component propellant mixture of Lp^/LHo/RP'l • 
The upper curve in these figures is the result of this three - 
component mixture with the ratio of to RP-1 held constant. In 
arriving at the numerical values, it was assumed that all Li^ always 
reacts, and as much of the RP-1 as can be supplied with LX^. Atmos- 
pheric oxygen could also be included if desired without any particular 
difficulty. 
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Figure 2A 
(Chart I) 
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Figure 2B 
(Chart I Detail) 


OxidUer to Fuel Ratio 
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The lower curve is the result of a two component mixture 
IX^/RP-l again here without atmospheric oxygen contribution. This 
curve is applicable to a two component mixture or could be considered 
the condition after all the II^ of the three-component mixture has 
evaporated. 

Thus any three -component LC^/If^/RP'l mixture will have its 
starting point on the upper curve and will, due to evaporation of both 
IH2 and the LC^, follow a path from the upper to the lower curve toward 
the origin; this is if the reaction does not occur somewhere along 
this path terminating this process. The actual path depends upon the 
changes in the relative quantities of each component present. Two 
such paths are shown in Figure 2A and in more detail in Figure 2B 

How they are calculated will be explained later, but it might be 
mentioned at this time that they are for a mixture which was actually 
used in field experiments * * * * » 

One path assumes that the system is thermally isolated from the 
surroundings and the other path assumes that maximum thermal inter- 
action between the system and the surroundings occurs. 

The two paths are not as different as might be expected, indicating 
that the effect of the surroundings is minor. 

Yield Potential as a Function of Oxidizer to Fuel Ratio 

(Chart IT 

The explosive yield of the liquid propellants will depend not only 
on the quantity of energy released, but also upon the rate at which 
this energy is released. Because of lack of information as to the 
variation in the reaction rates as a function of the propellant 
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composition, it was assumed for these calculations that the reaction 
rate remains essentially constant throughout the LC^/ fuel ratios 
under consideration here. 

With this assumption which can, however, be replaced by reaction 
rate information as soon as it becomes available, and the information 
of Figure 2, the yield potential can be calculated and normalized in 
terms of the theoretical maximum. The results are presented in 
Figure 3. 

Mass-Fraction Time Relationship for IH^ and LO, (Chart 3) 

To be able to determine the actual paths as previously discussed 
and shown in Figures 2A, 2B, and 3, it is necessary to know the time 
variation of the I^/IX^ ratio. This is easiest to obtain from cal- 
culations of the quantities of IH^ and LC^ present at any time. 

The calculations are more or less standard, involving the prin- 
ciples of thermodynamics and heat transfer, but are very long and te- 
dious. They involve simultaneous heat balance and heat transfer rela- 
tionships with the proper heat transfer coefficients which allow, 
through step- by- step iterated calculations, the estimation of the 
quantities of cryogenics vaporized, escaping, or again condensed in the 
mixture, the quantities of fuel and oxidizer frozen and portions re- 
melted, and so forth. Some simplifying assumptions were made* wherever 
it appeared advantageous, in reducing the large amount of computations 
without appreciably affecting the results. Where quantities were 
encountered which had the same order of magnitude, but opposite sign 
and were relatively small, they were sometimes cancelled against each 


other . 
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These actions helped tremendously in reducing the scope of the 
necessary computations. 

Contact area variations based upon mixing studies both at the 
University of Florida in connection with the study of explosive hazards 
of liquid propellants, and information found in the literatme were 
used in the heat transfer equations together with the best available 
heat transfer coefficients to obtain the mass-fraction for IJ^ and 

u> 2 . 

The information needed and used, aside from that supplied by 
supporting studies at the University of Florida, are referred to in 
Table I; but only the results from the actual calculations can be pre- 
sented in Figures 4A and 4B. 

Yield Potential - Time Relationship (Chart 4) 

In the method for the calculation of the yield potential - oxi- 
dizer to fuel ratio relationship, time x was the common variable 
used; therefore it is possible to put a time scale on the paths as 
shown in Figure 2B. 

With these time scales on the paths of Figure 2B, these curves 
can be re-plotted, givinp the yield potential versus time relationship 
as seen in Figure 5. 

These curves represent the theoretical maximum yield which 
could be obtained at any time x from the above propellants due to 
the quantities of the constituents which are present at that time. 

One curve again represents the yield potential for the isolated system 
the othei; maximum theoretical thermal coupling with its surroundings. 
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Table I 


Li st of Literature References Used in Support of the Cal dilations 
for the Results Presented in Figures II-1A through I I-3B 


SUBJECT REFERENCES 

Average Chemical Formulas for Kerosene, RP-1 37, 43, 44, 46 

Average Heat of Combustion for Kerosene 37, 38, 45 

Heat of Combustion for Hydrogen 38 

Propellant Proportion used in Heat Transfer 
Calculations 6 

LAN/RP-1 Contact Area versus Time Data for 

LOg/RP-l Analogy 22 

Film Coefficients for LAN/RP-1 Interface 22 

Film Coefficients for LN /LH„ Interface and 

LN_/LH. Contact Area versus Time Data for 

LOj/LHfj Analogy 22 

Latent Heat of Evaporation for H_ and Specific 

Heat for GG 2 c 39 

Specific Heat for L JP-1 to simulate RP-1 28 

Latent Heat of Evaporation for O 2 39, 28 

Latent Heat of Fusion for 0 2 39 

Approximation of Latent Heat of Fusion for RP-1 46 

Approximation of Specific Heat of Solid Kerosene 40, 41, 42 









Figure 5 Yield Potential as Tine Function 
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Since the curves of Figure 5 give the yield for propellants when 
perfectly mixed to produce maximum yield, these results must be modified 
by the Mixing Function (The ratio of the mixed amounts to the total 
amount present at any time x ) . 

The Mixing Function (Chart 5) 

While the yield potential function as calculated above for a 
specific case established the actual quantities of the various consti- 
tuents present and the maximum theoretical yield, if all these consti- 
tuents are mixed most effectively, it does not give any infoimation as 
to the degree of mixing of the constituents. 

For example, at time x = 0 when the constituents just begin to 
mix, none of them are actually mixed and therefore an explosion could 
not be produced. Thus the Mixing Function is zero while the yield 
potential Junction is near its maximum. The product of the yield po- 
tential and the Mixing Function at this time gives the true or expected 
yield. 

The Mixing Function is essentially a hydrodynamic function, however, 
complicated by high rates of heat transfer. This makes the analytical 
approach difficult, and at least to start with, an experimental ap- 
proach for determining this function is more promising. This is true 
especially since questionable assumptions are not involved. 

Four methods have been developed in connection with the over -all 
systematic approach presented in this section. These four methods 
allow the detailed study of the mixing process and phenomena producing 
the Mixing Function of liquid propellants and have been used with 
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great success. In preliminary studies often applying two methods to 
the same experiment, these methods have independently produced results 
which are in excellent agreement . 

The detailed description of these methods and the various methods 

15 75 

of analysis and the results obtained by them are presented elsewhere * . 

Since they implement the approach suggested in this part of the report, 
they are briefly mentioned here and will be discussed in greater detail 
later. They are: 

1. Film Analysis* A high speed photographic technique, giving 
by use of mirrors a three-dimensional picture of the mixing process on 
the same film frame. Special analysis of these frames as to mixing 
profile, mixing volume, and turbulence factor allow the determination 
of contact area and degree of mixing. 

2. Wax Cast Analysis^ : By use of hot wax and cold liquids, the 

mixing process can be "frozen" at different stages of the mixing by 
varying the hot and cold temperatures. The "frozen" state of the mix- 
ing process can then be studied at leisure at any tune later. These 
casts can be analyzed for outside area by projection or coating methods; 
they can be serially sectioned to give the total contact area, tur- 
bulence factors (total contact area over profile area) , and sc on. 

3. Vibration Mixing Analysis* ^: This method consists of mounting 
a particular configuration on a vibrating table, simulating the various 
propellant components by partic les of different color, size, density, 
shape, etc., and after removing partitions, partially or completely 
shaking the system. The components will mix and the degree of mixing can 
be periodically checked at desired locations. Evaporation and other 
losses can be simulated by removing programmed quantities or numbers of 
particles at desired locations and prescribed intervals. 
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4. Thermocouple Grid Analysis* ^; This method of analysis employs 
a three-dimensional grid of fine thermocouples with each junction being 
monitored continuously. The traces give information regarding the mixing 
front, the degree of mixing at a particular point, the degree of tur- 
bulence at a point, the point or points of ignition, the time delay from 
the start of mixing (or time of failure) to ignition, the propagation 
of the reaction front, the propagation of the shock front, the separa- 
tion of the shock front from the reaction front, and so forth. 

Results from the above methods can be correlated and compared 
easily by simultaneously applying the different methods of analysis to 
the same experiment. These experiments provide information needed for 
the better understanding of the mixing phenomena of liquid propellants. 
They provide data as to the statistical reproducibility in seemingly 
identical experiments, the variation due to test configurations, and 
so on. 

The Thermocouple Grid Analysis method is the most powerful since 
it directly relates the mixing phenomena and the yield obtained all in 
one and the same experiment. It is, however, considerably more expen- 
sive than the others. Instrumentation for high-speed monitoring of 
the individual junctions is expensive and the reduction of the data ob- 
tained time-consuming. 

However, this Thermocouple Grid method is capable of taking measure- 
ments in liquid propellant mixtures from the start of failure up to and 
after ignition. If the grid is extended beyond the original boundaries 
of the propellant configuration, information can be obtained as to 
fireball growth rate, extent, temperature, shock wave strength, shock 
wave velocity, and so on. 

Further and more detailed discussion of these four methods of 
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analysis which can provide the Mixing Function-time relationship is 
left to a later section when they are used for rather detailed studies 
of the Saturn V Destruct System, early configurations of the Space 
Shuttle, etc. 

Only one of the results is presented here to show the complete 
procedure of the Seven Chart Approach corresponding to the spill 
experiments used as examples for ccnparison with the calculated numer- 
ical results. It is the mixing function presented in Figure 6, in 
this particular case obtained by method 3, the vibration mixing analysis. 
Since this method has no absolute time scale, a number of runs were 
made adjusting the amplitude and frequency so that easily measurable 
changes were observed in reasonable time intervals (about five seconds) . 
Since from theoretical considerations the maximum should occur at 
about seven seconds, this time was ascribed to the maximum point of 
the mixing curve. In this manner the absolute time scale was established. 

The reproducibility of this curve was still within plus or minus 
4 percent. The reproducibility became better as the mixing violence 
or quantity of mix increased. This fact was observed in all experi- 
ments whether simulated on the vibration table or with real liquids 
using the other methods. 

Delay and Detonation Times 

Probably the least understood phenomena of the ones discussed in 
this section are those controlling the delay and detonation times. 

Both these quantities will be discussed in considerable detail 
in the next section where they are evaluated, and detonation times 
calculated, based upon the new hypotheses proposed and referred to as 
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"Fireball Hypothesis" and "Critical Mass Hypothesis". 

The delay time is defined as the time from the start of failure 
to ignition, or the time from the start of mixing to ignition, which- 
ever is preferable. 

The detonation time is the time from ignition until the reaction 
reaches the boundary of the original propellant configuration. 

For the purpose of this presentation, actually measured delay 
times from the experiment are taken and statistically analyzed so as 
to establish the most probable yield value as well as the 95 percent 
confidence limits. More data is needed to establish these quantities 
with greater accuracy. 

For the test used for comparison here, the average delay plus 
detonation time was 3.3 seconds and the standard deviation for fixing 
confidence limits, 1.1 seconds. 

Expected Yield Function - Time Relationship (Chart 6) 

Having discussed the three groups of phenomena playing a role 
in producing the yield of liquid propellant explosions, the results 
obtained above can now be combined. 

Taking the Yield Potential Function calculated and presented in 
Figure 5 and then the Mixing Function presented in Figure 6, and com- 
bining them by multiplying corresponding time ordinates, the expected 
Yield Function is obtained. This result is presented in Figure 7 which 
shows the yield which could be expected at any time t if detonation 
did occur at that time. Only the curve for the isolated system is 
presented here but the other is obtained easily in the same manner. 
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This expected Yield Function has a plus or minus 4 percent variation 
in yield value due to this variation in the Mixing Function. 

The expected Yield Function has a characteristic shape starting 
at zero, increasing with a dip or double hump to a maximum value and 
then decreasi lg again. The dip or double hump is due to the initial 
proportions of the propellant components. 

The maximum is considerably below 100 percent since, with any 
appreciable quantity o«. liquid propellants, perfect mixing is almost 
impossible to achieve and furthermore, due to the time elapsed between 
the start of mixing, and best mixing, the y'eld potential falls below one, 
because of evaporation losses of the propellant components. 

Expected Yield (Chart 7) 

The last step in this series, to obtain the expected yield for 
liquid propellant explosions, is to superimpose the information de- 
veloped above to obtain the combined results. 

Figure 8 shows the results, the final step in this systematic ap- 
proach, with the expet -d Yield Function of Figure 7 modified by the 
delay and detonation times fixing an interval within which, statis- 
tically, 95 percent of all expected yield values must lie. 

The highest value for the expected yield fraction predicted for 
this test series, using 95 percent confidence limits, should be about 
0.43, the lowest 0.13. All values predicted by this approach should 
fall between 0 and 0.50. 

Another theoretical method for predicting the ignition time is 
presented in the next section and agreement between the theoretical 
and experimental results is excellent. 
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The seven steps of this procedure divide the complicated problem 
into three very distinct parts, which when combined give the following 
desired results: 

1. Maximum Theoretical Energy Release 

2. Yield Potential as a Function of Oxidizer to Fuel Ratio 

3. Mass-Fraction - Time Relationships 

4. Yield Potential - Time Relationships 

5. The Mixing Function 

6. Expected Yield Function - Time Relationship 

7. Expected Yield 

The three distinct parts are the determination of the yield po- 
tential, the mixing process analysis, and the ignition and detonation 
phenomena. 

The procedure to determine the expected yield for a particular 
problem is the most detailed of the ones presented here. It gives 
specific answers but considerable effort is required for this analysis. 
It will be used extensively in the last two sections of this report 
for the oaturn V Destruct System Analysis and the analysis of some of 
the earlier Space Shuttle Configurations. 

C. Critical Mass Method 
Introduction 

When large quantities of liquid fuels and oxidizers are brought 
together during experiments, or accidentally, the results are liquid 
propellant explosions. 
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Many different phenomena can triggei the explosion, such as 
flames, sparks produced by striking or breaking metal structures, hot 
materials or hot spots produced by slow chemical reactions of fuel and 
oxidizer, the breaking of crystals which are formed when one of the 
liquid components freezes the other and which are broken mechanically 
or by thermal stresses, or by static electricity which is a result of 
internal friction and which may produce a spark discharge. 

Many more possibilities could be cited but these additional sources 
do not change the basic picture. 

From the above it can be seen that as soon as there is contact or 
mixing of the fuel and oxidizer, ignition is possible if an ignition 
source is available. If not, the mixing process may proceed with more 
and more of the fuel and oxidizer mixing until an ignition source 
appears, either through external or internal action. 

This section discusses the inherent production of an ignition 
source due to the mixing of oxidizer and fuel, specifically for cryogenic 
rocket propellants. 

85 

It has been well documented over many years that wher luid layers 
move across each other that electric double layers are produced resulting 
in electrostatic charges and in very high voltage differences, if the 
liquids are good dilectrics so that the charges cannot leak off rapidly. 

If a gas bubble or other medium of lower dielectric constant is 
interposed between the highly charged liquid layers, electric break- 
down can occur (a spark jumps) which may well act as the source of 
ignition. 
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Since this electrostatic build-up is always present within the 
mixing regions of fuel and oxidizer of dielectric liquid rocket pro- 
pellants, electric discharge is inevitable when the voltage difference 
has reached the necessary value for the breakdown to occur across one 
of the bubbles produced by the boiling of the constituent with the 
higher boi 1 **‘g point. 

This study is rather general but was applied here to mixtures of 
and propellant combinations of main interest at this 

time. 

It will be shown that voltage build-up increases with the quan- 
tities or masses or propellants involved and when the ignition proba- 
bility reaches a value of one (or becomes a certainty) the quantities 
which have actually been mixed will be referred to as CRITICAL MASS. 

Fluid Plug Model 

Mich work was done with bulkhead type failures since they are 
likely ones and one of the more gentle methods of bringing the com- 
ponents together. This type was selected in connection with this work 
since large amounts of experimental data are available involving cryo- 
genics for either fuel, oxidizer, or both. The bulkhead failure mode 
was then used to fozmulate the Fluid Plxj Model. 

When one component of a fuel-oxidizer combination fa* , s into the 
second, this can be likened to a plug of fluid falling into the other, 
Tigure 9. 

When the fluid plug of one component enters the other component 
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heat transfer at the contact surface area, which can be calculated, 
results in large amounts of vapor being generated. The amount of 
heat transfer constantly changes since the contact area changes, and 
the resulting amount of vapor generated as a result varies constantly 
making the fluid plug bob up and down. 

The variable amount of vapor generated with time and the varying 
fraction of this vapor which is entrained at any one time in the fluid 
plug, provides the mass of the plug and the farce acting on it. The 
resulting motion is an oscillation with varying and generally decreasing 
amplitude due to the constantly varying driving force and the constantly 
varying mass of the fluid plug. 

The fluid plug motion can be expressed mathematically as 

a- g eg) (11) 

where 

. - Hk.iJ . . p g.!a> 

g g 

« 5 zA + f i* 

g g 

- r <V f V (12) 

and £F is the force acting on the fluid plug, m is the constantly 
varying mass of the fluid plug, y is the displacement from the 
instantaneous equilibrium position, and t the time. 

The density of the liquid in the plug is , the volume of the 
liquid in the fluid plug is V^, , the density of the vapor in the 
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fluid plug is p , and the volume of the vapor in the fluid plug is 

o 

Vg, , with g being the gravitational constant. V is the total 
volume of vapor generated and f is the fraction of it which is 
entrained in the fluid plug at time t . 

The volume Vp of the fluid plug can be expressed, if a cylin- 
drical configuration is taken, in terms of the fluid plug radius Tp , 
and fluid plug height H . 


V = irr H 
P P 


(13) 


Considering the heat transfer at the surface of the fluid plug 
q = h X AT At 

' Vs h fg (14) 

with q the heat transfer, h the heat transfer film coefficient, 

J. the average area of heat transfer during the interval At, AT the 
temperature difference, and h f g the latent heat of vaporization. 

From equation (14), the amount of vapor generated can be calculated. 

„ _ h T AT At 


’g 


p g h fg 


(15) 


Using the above equations and relationships, the equation of motion 
for the fluid plug can be written 

ZF « -pg mr p 2 y 


d 

“37 


K 

L 2 H - f i* t?\ i 

• 

k 

g 


[x 


( 16 ) 
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y * <Kx) (18) 

In the above equations is the density of the liquid sur- 
rounding the fluid plug, a Q the instantaneous equilibrium position 
or center of oscillations during dT and <f> is a function which has 
to satisfy the motion as dictated by the heat transfer and the re- 
sulting vapor generation. 

The above equation is very complex and for this reason simplifi- 
cations are made by assuming that H is constant for all the cyclic 
increments, and that £ , the plug density, is made to vary to com- 
pensate for this assumption. 

Solutions giving excellent agreement with experimental work can 
be made by further simplifying the conservation of momentum equation to 



( 20 ) 
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which has solutions of the form 

y ■ c * cos urr + c ' ' sin tux 


with 


to ■ 


Wb^p 2 

w*v 


1 

7 


( 21 ) 

( 22 ) 


c ' and c ' ' being arbitrary constants , and to the instantaneous fre- 
quency of oscillations of the fluid plug. 

The above equations can describe the actual motion of the fluid 
plug by successive application to fractions of cycles or successive 
time intervals. This can be indicated by 


d 2 y, 


dr 


d 2 y 2 


dr! 


+ Vi 


+ u 2 y 2 = 0 


d 2 y, 


n 


dx_ 


♦ to v * 0 
ir n 


n 


T 0 1 T i T 1 


T 1 T t 2 


( 23 ) 


T_ i < T < t_ 

n-l — — n 


with the solutions of 

y^ * Cj cos u)jT + Cj* sin to^x 


y 2 * c£ cos to 2 x + c£' sin <o 2 x 


( 24 ) 


y • c ' cos iiit + c" sin to x 
7 n n n n n 


and the boundary and initial conditions for each cycle increment. 
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The final conditions of one cycle increment fora the initial con- 
ditions for the next. 

t * tj , y - -Jj, a* - »i ■ 



H. 


l 


(25) 


where the subscript "i" refers to the instantaneous values for the 
cycle increment under analysis. 

From the above the total displacement (submergence) of the 
plug is 


Y i - a i * n 


p pi 

— - -in V 


(26) 


And the volume mixed at time is 


V. 

l 



(27) 


And the contact area at time is 


A i * (lTr p. + 2lir p i Y i ) F T 


(28) 


where F,p is a turbulence factor which may be added if the contact 
surface is very turbulent or stirred up. 

In using the Fluid Plug Model, the equations can best be solved 
by trial and error making the nece sary assumptions and then varying 
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them until the loop, which has to satisfy the physical phenomena, 
closes. The motion of the fluid plug produced by the density changes 
in the fluid plug must agree with the heat transfer producing these 
density changes which in turn is a function of the fluid plug motion. 

The detailed steps involved in the method are: 

a) The determination of r , H, u Q from the physical configur- 
ation and failure mode. 

b) The selection of At or fraction of cycle to be analyzed 


c) The calculation of V 


^LCyRP 

l ' L< Y LH ; 


®i 


7776 Btu/hr-ft - °F 
1140 Btu/hr-ft 2 - °F 


d) The selection of f (or calculation based upon bubble rise 
velocities) . 

e) The calculation of to. and a^ . 

f) The use of the initial and boundary conditions to calculate 
the amplitude of instantaneous oscillation for the At under 
consideration. 

g) Calculation of the penetration. 

h) The determination of the volume mixed V . during the 
interval under analysis. 

This procedure is then repeated as often as needed or desired to 

obtain the Mixing Function V . versus t 

mix 

The analysis described above, using the Fluid Plug Model for 
a; proximate ly half cycle increments is presented in Figure 10. 



Volume Mixed, cubio inches 
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Fluid Plug Height S inches 
Radius 1.5 in. 


40 U 



Time, Seconds 

Fig. 10 Fluid Plug Mixing Function 
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Also plotted on this graph is the experimental result of the 
Mixing Function for LN 0 /RP-1 . More discussion on the experimental 
work will be found in a later section. 

E lectrical Analogy 

As has been discussed in the literature for years, motion of 
fluid layers will, through the relative motion and friction, produce 
charged regions and set up potential differences within the fluid if 
these charges are prevented from leaking off too rapidly. 

Since the motion of the fluid plug is described by the Fluid 
Plug Model developed, this motion should make it possible to describe 
also the electrostatic charge and voltage generation. Substitution 
in the equations of motion of the equivalent corresponding electrical 
quantities for the mechanical ones will provide a relationship between 
the charges generated and the motion of the fluid plug, the mixing 
and the heat transfer. 

The frequency of the motion in the basic equations remains the 

same. 

While the physical motion, is however, of oscillatory nature, 
the charge generation and build-up is cumulative, or in other words, 
independent of the direction of the motion of the fluid layers. 

The expression becomes 

d 2 Qi 2 

5“ ♦ o.f Q. - 0 


dt 


( 29 ) 
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when the same basic assumption is made for simplification purposes 
as before. The general solution is again of the foim 

Qi * Dj cos WjT + E^ sin u.x (30) 

Q. is the instantaneous charge generated and and are 
arbitrary constants, uk is the one determined from the physical 
motion of the fluid plug. 

These equations allow the calculation of the charge generated 
during each time interval or cyclic increment, with the constants 
determined from either experiment or physical and electrical properties 
of the fluids. 

The value of the accumulated charge, which is the quantity which 

eventually provides the ignition source through generation of the 

required potential, is obtained from 

Q - Z |Q.| (31) 

i 1 

where Q is the total accumulated charge calculated from the sum- 
mation of the absolute values of the charges generated for the cyclic 
increments. 

Figure 11 presents two curves of charge generation calculated in 
this manner from, the fluid plug model . The lower curve is calculated 
for a fluid plug 3 inches in diameter and 5 inches high. It corres- 

T 

ponds to the use of 140 in of in a 504 bulkhead failure arrange- 
ment (d/D = 0.5), where d = diameter of the hole and D * diameter of 
the tank. 

The upper curve represents twice the quantity or a fluid plug 
twice as high. It is shown that the doubling of the mass produced 



Charge Generated, Coulombs. x 10 
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Fig. 11 Theoretical Charge Generation, Fluid 
Plug Model. 
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essentially twice the charge. Experiments over two orders of magni- 
tudes demonstrated a linear relationship. 

In all the experiments carried out with both organic (hydrocarbon) 
fuels and as well as inorganic liquids such as distilled water 
and LN^j it did not seem to make any difference whether the mixing to 
produce the electrostatic charges was carried out in steps or the total 
mass was brought together all at once. This is only true if the leak- 
age was negligible since the step by step mixing took much more time 
to produce the same charges. 

Since the fluid layers act as plates of a capacitor, and the 
electrode screens are imbedded in the fluid layers, the relationship 
between electrostatic charge, capacitance C, and voltage V can be used 

Q = CV (32) 

With the Q calculated from the equations of motion and the 
electrical equivalencies, and C from the physical configuration, V 
can be obtained from equation (321 . 

The voltage build-up obtained in this manner is presented in 
Figure 12. It is presented here as a function of mass which is the 
form in which it will be used later for determining the Critical Mass. 
The experimental verification is presented in the next section. 

Critical Mass 

The hypothesis proposed here, and believed to have been verified, 
is that in the absence of any external ignition sources a mixture of 
fuel and oxidizer will, through the nature of the mixing process, pro- 
duce electrostatic potentials. These will reach vuues at which 



O Average experimental value for 
many experiments 

Electrode Spacing 1.25 inches 




Theoretical, based upon Fluid 
Plug Model 
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breakdown through sparking occurs supplying ignition sources for the 
explosive mixture. 

As a consequence of the above, it is impossible to mix unlimited 
amounts A cryogenic rocket, propellants since the voltage build-up 
becomes larger as the quantity mixed becomes larger, and finally with 
a certain limiting quantity mixed, ignition becomes a certainty since 
breakdown will occur. This quantity is referred to by the writer as 
CRITICAL MASS. 

Basically what this means is that an explosion can occur at any 
time with any quantity of explosive mixture, but the probability of 
it occurring increases and finally becomes a certainty at a particular 
quantity, and it is impossible to mix quantities larger than this 
limiting value of Critical Mass without an explosion occurring. (The 
quantity of 'ritical Mass, however, varies with the energy of the mixing 
process. This will be discussed later in the report. 1 Ignition is 
certain to occur as soon as the Critical Mass is reached. 

According to the literature, it takes 0.02 millijoules to ignite 
hydrogen, and under unusual conditions it has been ignited by 1300 
volts^ Usually it takes higher potential differences, and vol- 
tages of 14,000 to 20,000 volts are quoted®^ Seme authorities, 
according to these sources, consider it necessary to have an electric 
field strength of 76,000 volts per inch before spt king can occur. 

Using the last lumber, * . id require about 19,000 volts for 

discharge across a 1/4 inch bubble, the most prevalent size in the 
experiments . 
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Further, using the knowledge that 4 volts are produced at an 
average, for each 200 ml of LJ^ and recognizing that the physical - 
electric- 1 properties needed here of LN^, IJ^, and LO^ are similar, 
it can be estimated that for a bulkhead type failure the Critical 
Mass, or mass which when mixed will produce this voltage, is about 
2300 lb. for LC^/U^ and approximately 2800 lb. for LO^/RP-l. 

Trusting that the above results are correct and applicable, one 
can estimate the maximum expected explosive yields for such failures, 
obtaining about 3 percent for the S-IVB PYRO experiment and about 11 
percent for the 25,000 lb. LC^/RP-l bulkhead type explosion experi- 
ments. These yield values are in excellent agreement with the reported 

119 

experimental values 

Variation of Critical Mass with Mixing Energy 

When LO^ awl RP-1 or LC^ and LI^ are brought together the dif- 
ference in temperature produces rather rapid heat transfer. This heat 
transfer produces large amounts of vapor which through its sudden dis- 
placement of liquid produces rather violent mixing of the fuel and the 
oxidizer. 

With the propel lai'. combinations under these conditions, at least 
the boiling will always be present, and a fraction of this energy will 
produce agitation thus contributing to the mixing process. 

More gentle mixing is not possible, and mixing energies smaller 
can only be produced when the propellant components are brought to 
temperatures less different from each other. 
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If the mode of failure is changed so that the propellants are 
driven together with greater energies, the mixing will take place 
faster. The charge and voltage generation will also occur faster 
but does not relatively increase as much as the mixing rate. 

From these considerations it can be seen that if propellants can 
be mixed very gently, say drop by drop, infinite quantities could be 
mixed, which however would take infinite time, without producing 
ignition. On the other hand, if the propellants could be brought 
together with infinite energy, they could be completely mixed in zero 
time, again allowing the mixing of infinite quantities before ignition 
can occur. 

This analysis then indicates that the Critical Mass, when plotted 
against the Mixing Energy, which is supplied by both the mode of 
failure and the boiling process, approaches infinity both as the 
mixing energy approaches zero and also as it approaches infinity. 

Thus, the plot of Critical Mass versus Mixing Energy is a distorted 
"U" curve, and all the real values which can be encountered will fall 
below this curve. 

From mathematical analysis, the Critical Mass function can be 
divided into two characteristic branches, the one to the left of the 
minimum and the one to the right. The minimum point represents the 
value of Critical Mass which is the maximum value which could be pro- 
duced ; th the boiling energy of the bulkhead type failure. Since 
this value is easily calculated for any quantity of propellants, all 
other cases are normalized with respect to this energy value. 
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From experiments with non-cryogenic mixtures, it is seen that 
the time it takes for the voltages to be leached which could produce 
ignition decreases with an increase in the mixing energy; conversely 
it approaches infinity when the mixing energy approaches zero. The 
relationship below closely describes this behavior. 

x*E * c 3 (33) 


where t* is the Critical Time, or time it takes to build up the 
voltage necessary for ignition, and E is the Mixing Energy. 
is an arbitrary constant. 

From the Fluid Plug Model the mass mixed at any one time can be 
expressed as 


M 


c^ sin urt 


„ „ 3 3 _ 5 5 

C^OJfT - C^W T C^U) T 

— 3T + T! 


= c 4 wt 


(34) 


and therefore 


M* * c.wr* 
4 


and from (33) 


c 4“j 


( 35 ) 
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This equation is seen to give the largest values for Critical 
Mass at low values of Mixing Energy. Since with cryogenic propellants 
under normal conditions, the Boiling Energy is the minimum which 
produces the mixing. Equation (35) is not the most significant con- 
tribution. 

Again from the Fluid Plug Model and laboratory experiments with 
cryogenic fluids, it could be seen that twice the energy supplied for 
the mixing process produced essentially twice the mass mixed for the 
same time interval, three times the energy produced three times the 
mass mixed, etc. 

This can be stated as 

M = c 5 Et (36) 

From the previous work at the University of Florida with 
LN^/PP'l and the mixing studies carried out by FYRO with ** 

was found that 

M * c 6 t 2 (37) 

Setting equations (36) and (37) equal, the relationship between 
and E is obtained. 

c 5 E c 6 t2 

t = ^ E (38) 

c 6 

Substituting this last obtained relationship into equation (36) 
gives the significant M and E functional relationship for Mixing 
F* orgies larger than the Boiling Energy. 



SI 


C 5 

M = c,. E E 
5 c 6 

- c 2 E 2 (39) 


The two relationships developed (35) and (39) can now be com- 
bined for a general equation relating Critical Mass and Mixing Energy. 


M* 




(40) 


For the analysis here, the Critical Mass - Mixing Energy rela- 
tionship was normalized with respect tu the Boiling Energy. 



Evaluating the constants c^ and c^ from the Critical Mass 
detennination based upon the charge and voltage generation, the 
Critical Mass Function takes on the foim: 



and can be plotted. Figure 13. 

Figure 13 presents a plot of Explosion Mass (the amount of mass 
in any particular case which actually takes part in producing the 
explosion) versus the Mixing Energy Ratio (Actual Total Energy pro- 
ducing the mixing divided by the boiling energy) . On a plot of this 
sort all explosions of liquid propellants can be recorded for which 
enough information is available, such as yield values, and in addition 



1 PYRO 278 

25,000 lb. L0 2 /RP-1 

2 PYRO 282 

25,000 lb. L0 2 /RP-1 

3 PYRO 301 

94,000 lb. W 2 /RP-1 

4 PYRO 285 

25,000 lb. L0 2 /RP-1 

Drop Test - 44 fps / 

5 S-IVB Test 

lo 2 /lh 2 / 

Other points plotted represents / 

experiments and failures, both with / 

L0 2 /LH 2 and LOg/RP-l. / 


CRITICAL MASS 


... * _ « . . Total Energy Producing the Mixing 

Mixing energy Ratio, Epergy p r oftdri by Boll In F On l y 


Fig. 13 Explosion Mass - Mixing Energy Relationship 
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the upper bound of these values can be plotted as the Critical Mass 
Curve. 

All actual masses taking part in an explosion of liquid rocket 
propellants must fall below this curve. 

4 

It is seen from Figure 13 that all known values of explosions 
for liquid rocket propellants for vdiich data is available fall, when 
plotted in this manner, below the Critical Mass Curve and only a few 
values approach the curve. 

Many actual cases have been analyzed, but only a sample of them 
are plotted here. For these cases the boiling energy for a bulkhead 
failure with equal quantity of propellants was calculated for refer- 
ence and then the actual boiling., potential, and kinetic energies 
available for the mixing process were considered. 

It is believed that this new concept of Critical Mass for liquid 
propellants is very useful, indicating that there are self- limiting 
phenomena in the mixing and explosion processes of liquid rocket pro- 
pellants and that there exists an upper limit to the size of explosion 
which can be realized. 

For total quantities of propellants less than the Critical Mass 
explosive yields of 100% are possible, although the absolute size 
of the explosion is smaller than that produced by the Critical Mass. 
Furthermore, maximum explosive yield values for liquid propellant 
explosions depend upon the mode of failure which, in part at least, 
provides the energy available to bring the propellants together. 

With large quantities of propellants it is physically impossible 
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to get all the propellants involved in an explosion since ignition 
occurs at the latest as soon as the Critical Mass is reached in the 
mixing process, while the remainder of the propellants take part 
only in rapid burning and the foimation of the fireball. 

The hypothesis of this self-limiting process and the fact that 
electrostatic charge and voltage generation provide this physical 
limitation, has been shown in the laboratory. It has further been 
verified by checks against the explosive yields of large propellant 
quantities which were produced both by tests and missile failures. 

The mathematical Fluid Plug Model presented here demonstrates 
how the Critical Mass can be predicted from theoretical analyses for 
a desired mode of failure. 

The results and prediction;* based upon the Critical Mass 
Hypothesis agree with actual experiments and actual failures and also 
with results obtained by the Mathematical Model rnd the Seven Chart 
Approach. 

D. Fireball Hypothesis 
I ntroduction 

Sections A, B, C, above give some of the characteristics of 
liquid rocket propellant explosions. Section A presents the mathe- 
matical model which is a statistical treatment of the overall char- 
acteristics culminating in prediction of the explosive yield that 
one can expect. Section B, the Seven Chart Approach, gives details 
of how the fuel and oxidizer come together, mix for certain modes of 
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failures, and the expected yield for any time following the failure. 
Section C, the Critical Mass Method, sets the upper limit for liquid 
rocket propellant explosions indicating that the mixing process has 
built-in self-limiting characteristics so that very large explosions 
with large liquid propellant quantities are not possible. 

The present section takes the phenomena from the time of ignition 
through the formation of a detonation front at the reaction front, 
which in turn produces a shock front. The reaction front and the 
shock front behavior which produces a fire! all and a following com- 
bustion product cloud is discussed. 

After ignition occurs, the reaction of fuel and oxidizer produces 
a reaction front which develops into a detonation front producing a 
shock wave which has a maximum velocity due to confinement of the 
surrounding liquid. The velocity then decreases with the decrease in 
confinement as the fronts come closer to the surface of the container. 
When the shock front, pushed by the reaction front, reaches the edge 
of the mixture and enters into the inert liquid, it moves ahead of 
the reaction front and is reflected back and forth between the wall 
and the reaction front. Each pressure impulse can actually deflect 
the walls while this occurs. When finally the reaction front and 
shock front are at the surface together, they suddenly break through 
with an increase in velocity due to the transfer of momentum to a 
lighter fluid and emerge as a shock front and a reaction front which 
forms the edge of the fireball. The shock front is very rapidly 
attenuated as it moves through the atmosphere. The reaction front 
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also is attenuated rather rapidly forming the typical hesitation hump. 

For the purpose of this discussion, the problem can be divided 
into 5 parts which together form the "fireball hypothesis". These 
5 parts are: 

I. The region where ignition produces phenomena that develop 
in co the detonation phenomenon. 

II. The region where the reaction front and the shock front travel 
together through the propellant mixture. 

III. The inert liquid region where the shock front and the reaction 
front separate and where the shock front is reflected back 
and forth between the walls of the container and the reaction 
front. 

IV. The liquid propellant-air interface, actually a region where 
the liquid boundary begins to move and where the velocity of 
the reaction front forming the fireball and the velocity of 
the shock wave increase sharply due to the entrance into a 
lower density medium. 

V. The region in which the shock wave travels through the at- 
mosphere and where the fireball grows and develops separately 
far behind the shock wave. 

Figure 14 is a sketch of the hypothesis showing the five regions 
for one case of confinement and yield, etc. These regions are dis- 
cussed in detail in the following pages. 

It might be well to mention that the scales chosen for presenting 
the five regions in Figure 14 are different for each so as to be 
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able to show the variations occurring in each region. Region 4 is 
actually very small; regions 1, 2, and 3 make up the physical space 
of the liquid propellants, the relative size being a function of the 
explosive yield obtained from the propellants. Region 5 is by far 
'nrgest. 

Also on a time scale, the phenomena in regions 1, 2, 3, and 4 
will happen in a matter of milli-seconds and the flash in micro- 
seconds, while those of region 5 ordinarily stretch over many seconds. 

This last fact is also the reason why measurements are available for 
a good part of region 5, while they are not for regions 1, 2, 3, and 4. 

The other sections report work covering another phase of the 
problem of liquid propellant explosion hazards, mainly the Mixing 
Function - an important factor in the prediction of explosive yields, 
and propose a method which is believed to be able to give experimental 
information in regions 1, 2, 3, and 4. 

The Five Regions 

Region I - The Region Where Ignition Produces Phenomena that Develop 
Into the Detonation Phenomenon. 

Some time after liquid propellants are brought together, either 
intentionally or through failure, ignition may occur. At this time 
of ignition a certain fraction of the total propellants involved if 
mixed and ready for reaction. Since mixing continues with time, but 
also evaporation losses occur where cryogenics are involved, the 
ignition time has a pronounced effect upon the explosive yield^*^ *75,104 
The time of ignition may be essentially a constant for hypergolics or 
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or a random function for cryogenic propellants. 

The reaction between the propellant components thus initiated 
will progress rather quickly among those molecules that are ready 
to react. This may be expressed simply as 



(43) 


N * number of molecules ready to react 
t * time 

k * reaction factor 

Since the molecules that are ready to react at time,T , must 
be in "contact" with each other, they can be thought of as located 
on a "contact area" or "surface" produced by the mixing process. Thus 
the above equation can also be expressed as 

gf - kA (44) 

A = contact area 

This contact area, or a surface proportional to it, can be mea- 
sured experimentally for a particular missile configuration and mode 
of failure as a function of time. For the purposes here and from a 
theoretical point of view, it could be considered spherical for a 
first approximation, thus assuming the reaction progresses in all 
directions, which will certainly be true in the initial stages. 

Thus we can write 

47rr y 2 2 

A - — - 4mr (45) 
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r = equivalent radius of actually spherical surface 

r = radius of actual total contaet area of which only a 
fraction y reacts 

y = yield 

Little information is available on the reaction factor k in 
equation (44). Let us assume that it is a function of time and that 
it can be expressed as 

k = k 0 t a (46) 

where k and a ave ranstants. 
o 

With this inf tion and the above assumptions. Equation (44) 
can be integrated. 

(A) For a = 0 in equati :*n (46) , equation (44) can be integrated, 
and utilizing equation (45) one obtains; 



(47a) 


(47b) 


Where v is the velocity in Region I and r 9 the corresponding 
y 2 

radius at which this velocity is reached. 

Plottinp these parametric Equations (4' a) and (47b) as v 
versus r, the part of the r .reball Hypothesis falling into Region I 
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is obtained (Figure 15A). 

(B) Repeating the operation of part (A) but with a * 1 , 
making the reaction factor a linear function of time instead of 
leaving it a constant as in (A) we obtain 



r 2 



(48a) 


(48b) 


Plotting these parametric equations as in (A), Figure 15B is obtained. 
Other functions can be selected for the reaction factor and with it 
other velocity rise rate curves can be obtained. The decision as to 
which relationship represents the true case best (for particular 
propellants) will have to wait until experimental information, either 
direct as velocity measurements or indirect in terms of contact area- 
time measurements for some of these cases, becomes available. 

The y?eld ; y, was considered constant or at least an average 
value in the above calculations. Any other function can be selected 
in this analysis as soon as there is some justification foT it. 

It is believed that the above approach gives better insight 
into the happenings in Region I and allows for the expression of 
these happenings in a satisfactory manner. 
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Region II - Hie Region Where the Reaction Front and the Shock Front 
Travel Through the Liquid Propellants Together. 

In Region II, if the propellants were properly mixed, uniform 
or smooth propagation of the reaction front and the shock wave 
would occur. Since this is most likely not the case, the traces 
mist be considered average curves since actually they would have 
small steps, like a stairway, superimposed upon them. These finer 
points can be added later on, after the overall hypothesis is de- 
veloped and assuming that enough information is available to do this. 

If the physical system is very small or if the system is es- 
sentially unconfined, thus not capable of supporting pressure gra- 
dients , then the waves would travel with the velocity of sound as 
soon as these veloc ties are reached. Under these conditions the 
velocity could be considered constant in the Region II with the 
reaction or detonation and shock fronts traveling together. 

Actually confinement, especially in the earlier stages, will 

build up pressure and temperature due to the reaction taking place 

and will further increase the front velocities. Assuming the degree 

of confinement in terms of the masses surrounding the reaction, the 

pressures and temperatures can be estimated, and from them the -.rave 

7 

velocities. Studies ha\e shown however thrt the velocities in 
liquids do not increase very fast with increase in pressure so that 
the original assumption of essentially constant velocity in this 
region seems justified. 

As the traveling waves approach th^ boundary of the liquids, 
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the confinement decreases, and so the velocity differential between 
unconfined, or sonic, and the actual velocity decreases. Therefore 
a maximum velocity will be reached in this region which is a function 
of confinement , and which is also dependent upon the missile config- 
uration and yield for specific propellants. 

Calculations have been made for confinement but seem to be of 
lesser importance than other considerations. 

c igure 16 represents schematically, for a particular condition, 
the velocity rise and following decrease in Region II . 

In Figure 16, detonation velocity is reached at point A. In 
an unconfined system the velocity may be considered constant until 
the edge of the mixing region is reached. Thus the line AB would 
represent this case. If confinement is considered, which is naturally 
highest on the left hand side of the Region II, the wave velocities 
will still increase until, due to continuous dec. ease in confinement, 
they reach a maximum, and then wi*l decrease to a value above or 
equal to the unconfined system, depending upon what the confinement 
is at the liquid-air boundary. Thus the actual case will be more 
closely represented by the curve A-A'-B*. 

The calculations for these curves in Region II are relatively 
simple. They can be carried out provided the equation of state, 
heat generation and losses from a certain region, and constant volume 
processes for complete confinement, or appropriate volume changes 
for the various degrees of confinement, are considered. 

Verification of this paxt of the hypothesis is presented in the 


next section 
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Region III - Separation of the Shock Wave and the Reaction Front and 
Travel of Bach Through the Inert Liquid 
When the shock wave, driven by the reaction front, arrives at 
the edge of the mixing region, the shock wave will enter this region 
being attenuated through the new medium. The reaction front not 
being driven as it was in the mixing region, will be attenuated much 
more rapidly. The shock wave moving ahead of the reaction front, 
will arrive at the container wall giving an inpulse to this wall, and 
upon being reflected it will reverse its direction and move toward 
the reaction front. When the two fronts meet, the shock front will 
again be reflected in the original direction and thus will bounce 
back and forth between the walls and the reaction front. While the 
shock wave gives continued pressure impulses to the container walls, 
the wall may actually start expanding and moving out. This phenomemon 
has been observed in high speed photography of exploding missiles. 

Figure 17 shows more detail of Region III and a schematic indi- 
cation of the phenomenon occurring. The actual velocities of both 
the she k front and the reaction front can be calculated with the 
shock front taking up energy from the reaction front every time it 
is reflected from it. Thus the velocity of the shock wave will 
increase and the distance traveled between reflections will rapidly 
decrease. These phenomena continue until both the shock wave and 
reaction front arrive at the already moving boundary together. 

Region IV - The Liquid Propellant -Air Interface, Actually A Region 
Where the Liquid Boundary Begins to Move and Where the 
Reaction Front Forming the Fireball and the Shock Wave 
Emerge. 
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When the reaction front and shock front reach the liquid 
propellant-air interface which most probably has begun to move 
slightly, two distinct phenomena occur. 

A. The shock wave in the liquid is transformed into a shock 
wave in air. This is a wave-to-wave phenomena transition, 
essentially not involving mass transport. 

B. The reaction front, or detonation front, is transformed into 
a moving mass front. This is a transition from a wave phe- 
nomenon to a particle phenomenon. 

It seems that this may well be the reason why at this point the 
two phenomena, as for instance represented by the velocities, follow 
different paths. 

The shock wave experiences an increase in velocity as it crosses 
from the liquid into the air through possibly a small vapor layer. 

Very little information could be found about the transition of a shock 
wave from one medium into another and especially for the case where 
the compressibility is so different. 

The reaction front, when reaching the interface also as a wave 
phenomenon will now have to change into a particle phenomenon where 
the propellant particles have to move out forming the fireball boundary, 
thus involving tremendous amounts of mass transport. This transition 
must produce a rati;*'- abrupt change in the velocity in a rather small 
region. 

These phenomena are represented schematically in Figure 18 for 
Region IV. 


i 
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Region V - The Region in Which the Shock Wave Travels Through the 
Atmosphere as an Air Shock and Where the Fireball Grows 
and Develops Separately Far Behind the Shock Wave. 

After separation of the shock wave and the reaction front, each 
of these phenomena follow their own physical laws and relationships. 
See Figure 19. 

Fortunately some experimental information is available in this 
region on both the shock and the fireball. All this information, 
however, is for considerable distances from the liquid-air interface, 
and theory again will have to bridge this gap. 

Air Shock 

The attenuation of the air shock can be approximated by the well- 
known equations of compressible fluid flow and can be expressed in 
terms of the pressure ratio across the shock. 


where 


✓ 


c < 


Y + 1 

2y 








(49) 


v s = shock velocity 

c Q = velocity of sound at P Q 


Y * 1.4 for air 

P^/P. * pressure ratio across shock 


The pressure ratio can also be expressed in terms of the pro- 
pellant weight, the yield, the distance from the point of ignition 
and sane constants. 
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Equation (50) can now be substituted into Equation (49) to obtain 
the desired relationship, giving the air-shock velocity in terms of 
the distance from the point of ignition. Spherical geometry was 
assumed throughout in these derivations and whenever the fireball 
moved, r was taken as the radius of the fireball. This was necessary 
since otherwise information like wind velocity, etc. would have to be 
added into the above equations. Equation (51) is the desired result 
giving the air- shoe-, relationship for the fireball hypothesis in 
Region V. 


\T s f* * 

V 1 

fio b 

m 

[®) 1/3 \ . T - 1 tl 

V * V 

S 0 

2y 


[ r / y+T 


1/2 

> (51) 




The symbols in Equation (51) have the meanings as defined earlier 
with b and m constants . 

Equation (51) now allows the calculation of the velocity of the 
air shock at any distance from the missile. 


Fireball Boundary 

The fireball boundary can be calculated by utilizing the perfect 
gas relationships and considering either the spherical or hemis- 
pherical configuration. Then considering the heat generated through 
the chemical reaction processes and the heat losses, the necessary 
constants describing the process can be evaluated. Information such 
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as that found in the literature^ > 49 > 50 , 51 , 52 ver y helpful in 

this treatment. 

Further considering that the heat or energy released minus the 
amount used in raising the temperature and minus the amount lost pro- 
duces the kinetic energy which is observed in terms of the velocity 
at any time t or at any distance r . This is then the fireball 
velocity as expressed in Equation (52) . 


In Equation (52) , Q is the energy released and available for 
acceleration of the mass m . As can be seen the fireball velocity, 
that is the boundary velocity, is inversely proportional to the 3/2 
power of the radius. Q can be expressed in terms of the total weight 
of the propellants and the thermal yield. 

Figure 19 represents schematically the shock velocity and the 
fireball velocity in Region V indicating the separate paths which 
they follow. 

It might be well to mention again that in the presentation of 
the fireball velocity the actual center of the fireball was used as 
reference rather than the original point of ignition. 

If that had not been done, then the rise of the fireball and 
draft due to winds and atmospheric conditions would have complicated 
the presentation. 


2Q 1/2 1 C 

mrnii] ~ -rn 
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Complete Fireball Hypothesis 

Having discussed the five parts or regions of the complete fire- 
ball hypothesis, they can be combined to give the complete picture, 
which is presented in graphical form in Figure 20. The different 
regions are not plotted to the same scale but rather to a scale which 
allows the presentation of the variations in each 'on. This was 

done by enlarging Region I, and especially Region IV, and by shrinking 
Region III, and especially Region V. 

Figure 20 presents graphically the complete fireball hypothesis 
(not all detail shown) for a specific yield. Different yields would 
change it slightly, specifically the relative sizes of Regions I, II, 
and III. The higher the yield the smaller is Region I. 

The various regions of Figure 20 can be compared with the de- 
tailed regions and their discussion under the specific headings, 

Region I, Region II, Region III, Region IV, and Region V. 

Comparison of Results from the Hypothesis and Some Actual Data 

To show how thi hypothesis agrees wxth the sparse experimental 
information available, films of the S-IV Test were analyzed and 
analyses from the various reports were used. The hypothesis is plotted 
in Figure 20 and the available experimental points are superimposed. 

The equations used for plotting the calculated curves in Figure 20 
were Equations (51) for the air shcck velocity and (52) for the 
fireball velocity. 

These equations for the S-IV, where W ® 91,200 lb., y » 0.045 avg., 



Fireball and Shock Velocities, 



Figure 20 fireball Hypothesis Applied to the 
S-IV Experiment with Fxpe.-imental 
Points Superimposed upon the 
Calculated Curves 
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with their constants evaluated, have the following fora: 


Shock Velocity: 


v * 1100 
s 

Fireball Velocity : 


r 

j 14 r (91,200 x 0.045) 



1.4 


♦ 


1.17 


1/2 

l (51a) 


V FB " 


157.500 

372 “ 


80 


(52a) 


The results from these equations are the solid curves plotted 
in Figure 20. In this figure which is plotted to scale. Regions I, 
II, III, and IV are so small that not much detail can be shown. 


SIM4ARY 


In this section a hypothesis was discussed which seems capable 
of giving insight into the actual processes taking place from the 
time of ignition in a liquid propellant rocket explosion until the 
shock and the fireball have separated and dissipated. 

This hypothesis makes it possible to calculate curves for par- 
ticular missile configurations and yield estimates giving a package 
for complete analysis. 

The assumptions going into the analysis have been verified 
through actual experimental work the results of which are presented 
in the next section. 

The thermocouple grid analysis discussed earlier was the method 
used to supply the experimentation and verification. 
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E. Fireball and Post-Fireball Analysis 

Introduction 

After the detonation has occurred the behavior of the fireball 
from the explosion which is formed and then gradually changes into a 
combustion products cloud is of importance. How large a fireball is 
formed, what is its temperature and what are the pressures inside? 

To be able to obtain this information the knowledge of how the fire- 
ball comes about, how it cools and then changes into a combustion 
products cloud is essential. Thus, its behavior is really one of the 
last groups of phenomena or processes in a series. 

So for this phase of the work the knowledge of the fireball and 
combustion products .loud, volume -time, pressure- time and temperature- 
time histories have been assumed known and then the composition of the 
combustion phenomena has been determined. The composition of the 
fireball and of the combustion products cloud are important as well as 
their interaction with the atmosphere, expecially when toxic materials 
such as Flourine are used in the propellants. 

The volume -time, pressure- time, and temperature -time histories 

of the explosion from liquid propellants were chosen as input since 

20 21 

they may be determined theoretically * or may be measured in ex- 
periments thus givxiig a check on the theoretically detexmined infor- 
mation, with its statistical variations, etc. 

For this investigation the best infoimation available at this 
time has been used as input and then rather elaborate computer programs 
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have been used in obtaining the desired results . Homogeneity of 
the fireball and of the combustion products cloud have been assumed 
in all calculations and this seems to be a reasonably good basis 
since the turbulence of the reaction processes is great enough to 
traid to mix the different products well. 

With the volume-time, pressure-time, and temperature- time history 
available, the mathematical equations controlling these processes 
were set up and then computer solutions worked out to give the desired 
results. 

Only a portion of all the information which was generated by 
this analysis is reported here but it seems ample to show the method 
of approach and the kind of results which can be obtained. 

The fuel -oxidizer combination chosen here as examples for this 
investigation were ones which are used and some which may become 
important in the future development of liquid propellant rockets. 

They are: 

W 2 /W 2 

rp-i/lo 2 

^ 2 /LF 2 

RP-1/LF 2 

LH 2 /RP-1/L0 2 

LH,/LO, + 1% F 

L L + 51 F 

+ 10% F 

RP-1/L0, + 1% F 

z + 5% F 

+ 101 F 
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LH-/RP-1/L0- + 1% F 

c L + 54 F 

+ 10% F 

For the above combinations of fuel and oxidizer, assuming a 
quantity of propellants of about 100,000 lbs., the results which are 
all normalized are presented as follows: 

1. Fuel consumption versus time 

2. Volume of entrained air versus time 

3. Partial Pressures (of combustion products) versus time 
Partial Volumes (of combustion products) versus time 

4. Partial Weights (of combustion products) versus time 

Both the input information as well as the results are given in 
graphical form since it is believed that this method of presentation 
will give the maximum amount of information in the minimum amount of 
space. 

It might be mentioned that the method and computer program de- 
veloped are rather general and by no means restricted to the above 
fuel -oxidizer combination. 

Theory of Approach 

Equilibrium Composition of Chemical Reactions of Liquid Propellants 

Taking Place in the Atmosphere . 

The purpose of this phase of the research program is to theoret- 
ically determine the amounts of product gases foimed, as a function 
of time, as the result of a reaction involving liquid propellants and 
entrained air. This type of reaction is continuous since all of the 
available fuel does not react immediately and furthermore the resulting 
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fireball (which grows with time as more fuel reacts) continually 
entrains air. Given the initial amounts of fuel and oxidant as well 
as the volume - time history of the fireball (theoretically determined 
or as observed from high speed films), equilibrium compositions can be 
detemined. 

The equilibrium composition for the system of n products of 
reaction is determined by the simultaneous solution of n+1. equations 
consisting of the equations of mass balance, pressure balance and 
the dissociation equations involving equilibrium constants. 

Assuming a constant pressure process as well as an instantaneous 
reaction time and making use of either theoretically obtained or 
experimentally determined pressure- time and temperature- time histories 
of the fireball, a solution is found such that the total theoretical 
volume of the products of reaction is made identically equal to the 
total experimental volume by adjusting the fuel burning rate and/or 
adjusting the amount of entrained air. As a first approximation, it 
is further assumed that no air entrainment exists until all of the 
available fuel is burned. 

To efficiently meet these demands, a computer program has been 
developed. The program is general but limited here to the following 
reactants: liquid hydrogen, kerosene (RP-1), liquid oxygen, liquid 

fluorine, and air. Fifteen products of reaction were considered. 

The equations and method of solution follows. 
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(A 1 + A 2X)H 2 + (Aj + 0.209A 5 )0 2 + + 0.792A 5 N 2 + A^C 

Consider the following products of reaction; such that the right 
hand side of the equilibrium equation is 

F ( p l H 2° + P 2 C0 2 + P 3 H 2 + P 4 N 2 + P 5 F 2 + p 6°2 + ^7^ + 

+ pgNO + p l 0 °H ♦ P n H ♦ P 12 N ♦ P 13 F + p 14^ ^ p i5®) (53) 

where, 

N = total number of moles of products of reaction 
P ■ total pressure 

p^ = partial pressure of ith product 

The unknowns are p.(i ■ 1, . . .15), N. Hence 16 equations are 
necessary for a solution. The balance equations are: 
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(a) Pressure Balance 

P ’ P A * 0 

(b) Hydrogen Balance 

2(Aj + ■ p (2pj + 2pj + Py + Pjq + Pjj) ~ 0 


(54) 


(55) 


(c) Oxygen Balance 

2(Aj + 0.209A S ) - p (Pi + 2p 2 + 2p 6 + P 8 + P 9 + P 10 + p lg ) = 0 

(56) 

(d) Nitrogen Balance 

1.584 A s - p (2p 4 + p Q + p 12 ) - 0 (57) 


(e) Carbon Balance 

V " V (?2 + p 8 + p 14 5 * 0 (58) 

(f) Flourine Balance 

2A 4 - £ (2p s + p ? + p 13 ) = 0 (59) 

The above 6 equations can be reduced to 5 equations by eliminating 
p . Since we will be dealing with hydrogen (either LH 2 and/or RP-1) , 
p can be eliminated by dividing the 2nd equation into the last 4 
equations. The 5 equations are: 


• Pis ■ 0 


p - Pi - p 2 - 


(60) 
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<»1 - Y 2^ Pi - **2 + - *V« + Y lP 7 * V 8 " V 9 

+ <*1 " Y 2^ PlO + Y lPll • VlS * 0 W 

2 y 3 Pi + 2Y 3 p 3 - 2Y 2 p 4 + Y 3 p ? - Y 2 p 9 + Y 3 p 1Q + Y 3 p n - Y 2 p 12 = 0 (62) 

2Y 4 P 1 ' Y 2 p 2 + 2^*3 + Y 4 P 7 ‘ Y 2 p 8 + Y 4 P 10 + Y 4 P 11 " Y 2 P 14 * 0 
2Y s Pi + 2Y g p 3 - 2Y 2 p 5 + (Yg - Y 2 )p ? + Y 5 P 1Q + YjP^ - Y 2 p 13 * 0 (64) 
where 

Y x = 2(Aj + 0.209A 5 ) 

Y 2 = 2(A^ + A^) 

Y 3 = 1.584A S 

Y «’V 


The remaining 10 equations required for a solution are the dis- 
sociation equations. Tables of equilibrium coefficients are available 
(23) in terns of partial pressures rather than concentrations of the 
products of reaction. The appropriate equations are given below 


2H + 0 -*> H 2 0 
C + 20 -*• C0 2 
2H -► H 2 
2N -»■ N 2 


2 

such that p^ - KjPjjPjj. ° 0 

2 

such that p 2 - K 2Pl4Pl5 * ® 
2 

such that p 3 - KjP^j = 0 
2 

such that p^ - K 4 p^ 2 * 0 


(65) 

( 66 ) 

(67) 


( 68 ) 
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2F -*• such that pg - KgPjj * 0 (69) 
20 -*• 0 ^ such that Pg - K^p^g = 0 (70) 
H + F HF such that p^ - K 7 PjjP 23 = 0 (71) 
C + 0 -*■ CO such that pg - KgPj^js = 0 (72) 
N + 0 -*■ NO such that Pg - KgPj^PlS = 0 (73) 
0 + H -*• OH such that p 10 - KjoPijPj^ = ® (74) 


The equilibrium coefficients vary systematically with the 
temperature of the reaction. It is assumed that the products of 
reaction calculated at a particular temperature and pressure, are 
formed instantaneously. Hence, one need only solve the above 15 
equations for given values of P, T, and the amounts of reactants, to 
arrive at the equilibrium composition. 

Solution of These Equations 

Since the dissociation equations are non-linear, there exists no 
direct solution. The Newton-Raphson method is used to obtain a 
"trial and error" solution. 

(1) Initially, estimates of p^(i = 1, . . .15) are taken and 
each of the 15 equations is expanded in a Taylor's series 
about the estimated point, p^ . 

(2) Corrections to p^ are then found (Ap.) and the new esti- 
mates of p^ , given by p^ + Ap^ , are used in place of the 
initial estimates in (1) . 
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(3) The procedure is repeated until Ap^ becomes negligible. 
Consider a solution to two non-l^ear equations; f(x,y) * 0, 

g(x,y) « 0. 

Let the initial estimate of the required solution (x,y) by the 
point (x^yp. Expanding f , g in a Taylor's series about the point 
(Xpyp, then 


f(x,y) - o « f(x r yp + ^ (x^yp + ^ (x r yp . (y - yp + .... (75) 

g(x,y) - o « gfXj.yj) + ffcx^yp • (x-x x ) + I^Cxj.yp . (y-yp +....(76) 
or 


"-M* *♦$ • *♦••• 

Ag * |£ • x + |& . y + .. . (78) 

where 

Af » f(x,y) - f(x p yp » -f(x lf yp (79) 

Ag * g(x,y) - f(x 1 ,y 1 ) - -g(x 1 ,y 1 ) (80) 

and 

Ax * x - Xj (81) 

Ay * y - y l (82) 


Hence the non-linear equations have been transfomed into linear 


correction equations of the form 
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+ f x . Ax + f y . Ay = 0 (83) 

gCXj.yj) + g x .Ax + g y .Ay«0 ;j 

3£ 

where f » , f » f- , etc., and the derivatives of second order 

x dx y d/ 

and higher have been neglected. 

The equations are solved for Ax, Ay and new estimates are given 
by x^ + Ax, y^ + Ay. The procedure is repeated until the desired 
accuracy is obtained. 

The full procedure can best be demonstrated with an exanple . 
Example 

Consider a solution to 

f(x,y) « x 2 y + y 2 + 3 - 0 
3 2 

g(x,y) »x + 2xy + 4y 8 0 

then 

f x - 2xy 

f y • x 2 . 2y 

^ - 3x 2 - 2y 2 

gy = -4xy + 4 * 4(1 - xy) 

Let x^ * 1, y^ ■ -1, be the initial estimates. 

Then substituting into the linear correction equations 
3 - 2Ax - Ay - 0 
-5 + Ax ~ 8 Ay * 0 
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the solution to the above equations is 

x • 38/30 a 1.3 
y - 7/15 a 0.5 

The new estimates of x, y are given by x 2 , y 2 where 

x 2 = Xj + Ax * 2.3 
y 2 - y x + Ay * -0.5 

Substituting into the linear correction equations, then 

0.6 - 2.5 Ax + 4.3 Ay ■ 0 
9 + 15.4 Ax + 10 Ay - 0 

The solution to the above equations is 

Ax $ -0.4 
Ay 2 -0.3 


The new estimates for x, y become 
x 2 + Ax * 1.9 

y 2 + Ay ■ -0.8 etc. 

The correct solution is (2, -1). 


For more than 2 unknowns, the linear correction equations take 
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the form 


f(x,y,z,...) + f . Ax + f . Ay + f . Az + .... 

A / Z 

■ 0 

(85a) 

§(x>y >£>•••) + * Ax + gy * Ay + g^ . Az + •••• 

■ 0 

(85b) 

h(x,y,z,...) + . Ax + hy. . Ay + h z . Az + .... 

■ 0 

(85c) 

i(x,y,z,...) + 

* 0 

(85d) 

j(x,y,z,...) + 

* 0 

(85e) 


j 


etc. 

where the subscripted variable, f for example, represents the partial 

Jv 

derivative of f(x,y,z,...) with respect to x. 

Denoting equations (60) to (74) by (i • i, ...» 15); the 
correction equations are given by 

B 1 + ^,1 * Ap l + ^,2 • ^2 * ••• + Aus • A Pi5 * 0 


B 


IS + ^15,1 * Ap l J ^5,2 


Ap 2 + 

m 


* *15, IS 


Ap 


15 


0 ( 86 ) 


where A. . is the partial derivative of B. with respect to p. . For 
IfJ 1 J 

example, Ay ^ is the partial derivative of equation (7), i.e. By , 
with respect to p^. The equations are solved for Ap^ (i ■ 1 to 15) 
by first assuming initial estimates of p^ . Subsequent estimates of 
p. are given by p. + Ap. and the procedure is repeated until Ap ; 

approaches zero. 
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The coefficients of the correction equations (A. .) are denoted 
by matrix A and the constants are denoted by the vector, -B. 
Hence, in matrix notation, the set of linear correction equations is 
given by 


B = A . Ap (87) 

and its solution is given by 

Ap = A -1 B (88) 

where A * is the inverse matrix. 

Outline for the Fortran IV Computer Program 

The program is presently designed to handle nine sets of values 
of pressure, temperature and volume for a given propellant mixture. 
That is, equilibrium coefficients are incorporated into the program 
for values of temperature between 3000 K and 1400 K in 200 degree 
increments . 

Input Data 

The following information is required: 

(a) Weights of reactants, i.e., the total amount of fuel 
and oxidizer available. 

(b) Yield 

(c) Temperature of reaction 

(d) Pressure at which reaction occurs 

(e) Volume of products of reaction 



The following assumptions are implied: 

(a) Constant pressure process 

(b) Instantaneous reaction time 
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(c) No air entrainment until all of the available pro- 
pellants are used up. 

Procedure 

For each data point (i.e. for a given value of P, T, and V) 
the program determines the partial pressures of the products of 
reaction such that the theoretical volume of the product gases is 
identical to the given input volune. 

For the first data point, however, since no value of volume is 
available , the yield is used to determine the initial amounts of 
propellant burned and the partial pressures are then determined. 

For subsequent data points, the fuel burning rate is continually 
adjusted and partial pressures are calculated in turn so that finally 
the resultant theoretical volume becomes identical to the given 
(theoretically determined or experimental ly evaluated) volume. 

This latter procedure is repeated for subsequent data points 
until all of the available fuel is used up. From then on, air is 
added as a reactant combined with all of the available fuel in order 
to satisfy the "identical volume" condition. 

The program also converts the resultant partial pressures into 
the following: 

1 . Pound Moles 
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2. Pressure-Ratios, Mole-Ratios, Volume-Ratios 

3 . Pound Weights 

4. Weight -Ratios 

The fuel burning rate, the amount of entrained air, and the theo- 
retical volume for each data point are also determined . 

Subroutine Invert 

The subroutine solves the set of linear equations 

Ap = A _1 B (89) 

The input data card is 

CALL INVERT (A, NA, NAD, B, NB, NBD, DETERM, IERROR) 

where 

A = matrix of order NA 
B = vector having NB = 1 constant vector 
NAD = row dimension of A in main program 
NBD = row dimension of B in main program 
DETERM = dummy 
IERROR = dummy 

The output consists of A * placed in A, p placed in B, and 
the determinant of A placed in DETERM. IERROR is an error signal 
equal to 0 for successful inversion; equal to -1 for overflow, equal 
to +1 if no inverse is obtainable. 

The maximum size of A can be 100 x 100. 
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A - coefficients appearing in the correction equations 
B - constant appearing in the correction equations 

C - equilibrium constants 

P - partial pressure 

PR - partial pressure-ratio, mole-ratio, volume-ratio 

PT - total pressure 

T - temperature 

TNT - partial moles 

V - volume 

WMOL - molecular weight 

WT - partial weight 

WTR - partial weight-ratio 



RNAIR - weight-ratio of entrained air 

TN - total theoretical moles of products of reaction 

TNE - total experimental moles of products of reaction 
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TVOL - total theoretical volume of products of reaction 
WAIR - weight of entrained air 

X - number of carbon atoms in the RP-1 molecule, 

YIELD - percentage of fuel burned at time "zero" 

The fixed point variable, MA, is the number of experimental 
runs wit’-! combinations of U^, RP-1, LC^ and LF^. 

The Fortran IV program follows with an example of the output 

2 and 

entrained air. 

Input Information 

Many different quantities could have been chosen for the input 
information based upon which the desired fireball composition and 
atmospheric chemistry could be calculated. 

For this investigation the 
Volume-time history 
Pressure-time history 
Temperature-time history 
were taken as the principal input information. 

The reason for this choice was that another phase of this over-all 
program deals with the theoretical determination of these functions 
and most of all that it is possible to measure the above quantities 
and thus verify any theoretical results by actual field experimentation. 
This latter fact seems to be of extreme importance if theories are 
developed since without experimental verification they are of little 
use and certainly not much credence can be given to them. 


data for one of the nine data points using Ii^/RP-l/LC^/LF 
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Other f; ctors such as fuel burning rates , etc . were selected by 
20 21 

others * but the investigators of this project do not see how 
such quantities could be verified experimentally and therefore would 
remain assumptions throughout the work. 

As mentioned above, much work is being done on the determination 
of the volume-time, pressure- time, and temperature -time histories of 
the explosion phenomena from a theoretical point of view. Rather 
than wait for the results from this separate investigation and be- 
cause of contract conmitments it was decided to present the methods 
of obtaining the fireball and combustion cloud composition from such 
input data as mentioned above and for the present combine both theory 
and experimental information to obtain the most plausible functions 
at this time. 

A brief description of how the volume-time, pressure -time, and 
temperature- time functions have been determined for this report follows: 

Volume-Time History of Fireball and Combustion Products Cloud 
from Liquid Propellant Explosions 

The volume of combustion products produced by liquid propellant 
explosions transgresses a number of stages with time, changing in 
shape from one typical configuration into another. These stages can 
be observed in the high speed film records of such explosions and 
can be, in part at least, analyzed mathematically or theoretically. 
These major stages are: 

1 . Hemisphere 
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2 . Truncated Sphere 

3. Sphere 

4 . Pinched Sphere 

5 . Toroid 

The above 5 stages are distinct and can be observed in at least 
the larger explosions. 

Stage 1 . Hemisphere 

This stage is the earliest one which can be observed and is of 
relatively short duration. It involves a very rapid growth of the 
combustion products both along the ground and up into the atmosphere 
so that the shape can best be approximated by a hemisphere. The 
size of this initial hemisphere depends upon the yield of the liquid 
propellant explosion, the very rapid combination of the fuel and 
oxidizer so as to form detonation and shock waves. The larger the 
yield the larger the initial hemispherical fireball. 

Stage 2. Truncated Sphere 

Following the very rapid formation of the hemispherical fireball 
from liquid propellant explosions the hot combustion products begin 
to rise. This upward motion and the convection currents due to the 
bouyant forces undercut the rising mass thus forming a truncated 
sphere, in contact with the ground at the flat base. 

As the center of the mass rises the fireball changes more and 
more from the original hemisphere into a sphere, the shape which is 
attained when the combustion products become essentially tangent to 
the ground. 
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This stage in the development is referred to as "Lift Off", at 

20 21 

which most of the fuel seems to have been consumed * . 

Stage 3. Sphere 

Having attained essentially a spherical configuration at "Lift 
Off" the combustion products continue to rise as a rather turbulent, 
well mixing sphere which however gradually changes shape from the 
almost perfect sphere into the first slightly pinched and then rather 
pronounced pinched sphere. 

Stage 4. Pinched Sphere 

The change from the spherical configuration to the pinched sphere 
is rather gradual and then as the indentations become larger and 
larger, the appearance of the sphere is lost. A cross-section by a 
vertical plane through the center would give the appearance of a 
"Bar Bell". 

As this process continues the indentations will eventually touch, 
forming a toroid. 

Stage 5. Toroid 

From the time the toroid is formed the initial contact point of 
the indentation becomes a hole with the general configuration of a 
ring or doughnut . 

As this toroid grows in diameter the size of the hole increases 
but the volume now at this stage of development increases relatively 
slowly. 

Finally this well defined configuration diffuses into the 
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Figure 21 Typical Development Configuration Stages 

of Liquid Propellant Explosions 



Stage 1-. Hemisphere 



Stage 2-. Truncated Sphere 
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Figure 22 TYPICAL VOLUME-TIME FUNCTION FOR LIQUID PROPELLANT EXPLOSION 

PRODUCTS ( Yield =4.5 Percent ) 
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atmosphere losing its resemblance to any characteristic shape and 
being controlled to a great extent by the prevailing atomspheric 
conditions. 

Each of these stages as described above and schematically 
shown in Figure 21, takes a longer and larger part on the time scale. 
Stage 1 may occur in fractions of a second while the last stage will 
be a matter of minutes. 

Utilizing this 5 stage concept for the purpose of analysis a 
volume versus time curve can be obtained, either theoretically by 
the use of restricting assumptions or by the actual analysis of 
high speed film records of liquid propellant explosions. 

The variation is greatest in Stage 1 which is controlled by 
the yield while the statistical differences are rather small (but 
somewhat dependent upon atmospheric conditions) as long as the same 
quantities of propellants are involved and it is assumed that es- 
sentially all the propellants take part in the formation of the 
fireball and cloud^*^. 

Figure 22 shows the volume versus time curve for the S-IV PYRO 
experiment. The yield as reported was about 4 1/2% which is in 
agreement with the predictions of reference (1) . 

Similar volume versus time curves have been developed for the 
various fuel oxidizer combinations considered and reported upon here. 

Pressure-Time History of Fireball and Combustion Products Cloud 
from Liquid Propellant Explosions 

The pressure -time history as presented here and as used as input 
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data for the determination of the composition of both the fireball 

and the combustion products cloud was determined partially from 

preliminary theoretical considerations^^ and partially from the 

analysis of field data obtained by the liquid propellant explosion 

22 

program of project PYRO . 

The theoretical analysis was necessary for the early time 
processes since no experimental data is available ard the results 
were then checked and agreed with experimental re . in the later 
stages. 

In general it might be said that the pressure immediately after 
ignition rises very rapidly to very high values inside the missile 
due to the confinement of the propellants and the tanks, reaching a 
maximum somewhere as the reaction front progresses toward the boundary 
of the missile configuration. After this maximum is reached the 
pressure falls very rapidly to almost atmospheric conditions. 

From the time of "lift-Off" of the fireball which occurs at 
essentially atmospheric pressure 4 ^ the pressure drops very slowly due 
to the rise of the explosion products and the effect on atmospheric 
pressure due to altitude. 

The pressure -time history presented here for approximately 
100,000 lbs. of propellants was used for all the propellants reported 
upon here. 

Analysis of the sparse experimental information of liquid pro- 
pellant explosion experiments seems to support this general pressure- 
time history. The yield produced by the explosion will change the 
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early values of the pressure. Again for the analysis here a yield of 
4 1/2% was taken based upon the most likely value as given in reference 
( 2 ). 

The actual curve used here is presented in Figure 2 . If better 
information is to be used, an experimental program could be instituted 
to actually measure these pressures, an important reason for this 
choice of input information is because it allows theoretical deter- 
mination and experimental verification. 

Temperature -Time History of Fireball and Combustion Products Cloud 
from Liquid Propellant Explosions 

The third and last principal input information needed for the 
determination of the composition of the fireball and combustion pro- 
ducts cloud including air entrainment and atmospheric interaction is 
the temperature- time relationship. 

Again theoretical considerations und the available rough experi- 
mental observations of fireball temperature and variations with time 
indicate that the intial temperature is close to the maximum obtainable 
by the t -rticular propellants involved. Then, at low yields at least, 
since only a small part of the propellants take part in the initial 
stages of the fireball formation the reaction of the remaining fuel 
and oxidizer, both in the propellants as well as the atmosphere, make 

the temperature drop with temperature in an almost linear manner. 

20 21 

This is observed in theoretical work * and seems to be closely 

22 

approximated by the available experimental information . 
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EXPLOSION PRODUCTS ( YIELD « 4.5 PERCENT ) 
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This linear decrease of the temperature with time continues 

until the incandescence of the fireball ends often referred to as 

24 

the "duration" which can be approximated as shown in . 

For the purpose of analysis here it was assumed that the actual 
variation can be closely approximated by further linear decreases 
changing the slope to a value 1/2 the previous one for each subse- 
quent "duration" time interval. 

By this method a curve representing the temperature-time history 
of the liquid propellant explosion is obtained which from both theo- 
retical and the available sparse experimental observations seems to 
approximate the actual conditions. This again is taken here for low 
yield (4 1/2% in this case) liquid propellant explosions. 

Again it is believed that an experimental program can be designed, 
if desired, to obtain this temperature- time history for various cases 
and verify or modify the presently used information, which is presented 
in Figure 2 4 and again was used for all the fuel -oxidizer combinations 
analyzed and reported upon here. 

It should be mentioned again that the volume- time, pressure -time, 
and temperature- time histories were selected as the principal input 
data because it is felt by these investigators that this information 
which can be generated with appropriate assumptions theoretically, can 
be verified experimentally. In addition these volume-time, pressure- 
time, and temperature-time histories are of great interest to other 
investigators for various reasons. A number of groups are presently 
engaged in trying to measure pressures and temperatures within fire- 
balls and of combustion products clouds produced by liquid piopellant 
explosions. 
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FIGURE 24 . TYPICAL TEMPERATURE-TIME FUNCTION FOR LIQUID PROPELLANT 

EXPLOSION PRODUCTS ( YIELD » 4.5 PERCENT ) 
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It should also again be emphasized that in Uiis investigation 
homogeneity of the fireball as well as of the combustion products 
cloud was assumed. These assumptions seem to be reasonably well 
satisfied because of the tremendous turbulence observed within the 
fireball which tends to produce thorough mixing with a relatively 
short time. 

In addition to the principal input information, the volume-time, 
pressure- time, and temperature- time functions it is necessary to 
know 

4. The Type of Propellants 

5. Propellant Composition 

6. Propellant Quantities 

7. Yield 

4, 5, and 6 are easily obtainable as original data, while 7 is 
selected to obtain the results for this particular value of yield, a 
value which may again be dictated by theoretical considerations . 

Type of Propellants 

The type of propellants selected for this presentation are com- 
binations of fuel and oxidizers which are presently used in liquid 
propelled rocket systems or combinations which may become important 
in the future development of these rockets. 

The method however used is perfectly general and any propellant 
type and combinations could be analyzed in the same manner. 

The types selected for this presentation are: Ref. 75 
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W 2/ W 2 

rp-i/lo 2 

RP-1/LF 2 

LH 2 /RP-1/L0 2 

LH-/IXL + II F 

c 1 + 5% F 

+ 101 F 

RP-l/LO- + 1% F 

L + 51 F 

+ 10% F 

LH-/RF-1/L0, + II F 
L 1 + 5% F 

+ 10% F 

Propellant Composition 

The propellant type was outlined above with the composition of 
fuel to oxidizer chosen as follows: 


IU 2 /U> 2 

1 : 5 

by weight 

RP-1/L0 2 

1 : 2.25 

by weight 

w 2 /if 2 



RP-1/LF 2 



LH 2 /RP-1/L0 2 

1 : 2.6 : 5.86 

by weight 


In the combinations with Fluorine traces the weight ratios of 


the main constituents were the same as given above. 
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The chemical composition of the RP-1 was taken as C,, ,1L, ^ 

ll.o Z3.Z 

which was obtained from reference (25) . 

Propellant Quantities 

The propellant quantities were taken as 100,000 lbs. in all cases. 
This allowed the standardization of the pressure-time and temperature - 
time histories for the present analysis, since it seems that the quan- 
tity of propellants used has the major effect on the time axis of 
pressure and temperature. 

Yield 

The yield, the energy release as a fraction of the theoretical 

maximum, for these calculations and analyses was taken as 4 1/21 

2 

which from previous theoretical investigations and from experimental 
2 22 

observations ’ seems to come close to the statistically probable 
value. 

Again it might be mentioned that other values could be taken 
just as well without changing the method of analysis. The resulting 
compositions of the fireball and explosion products cloud would, how- 
ever, be different. 

With the input information as described above a number of cases 
were analyzed and many quantities calculated. Rather elaborate 
computer programs were developed for this purpose and the main program 
will be presented in the appendix. 

The results which seem to be most pertinent to this investigation 
are presented in the reference 75 , mostly in graphical form. 
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Results Obtainable 

Utilizing the data information as discussed above and the cal- 
culation and analyzing procedures outlined earlier, many important 
qunatities can be calculated. 

Because of the space limitations only one set of the more 
complicated combinations will be presented here. Many more have 

been worked out in connection with this work and have been reported 

75 

by Dr. Farber and his group in the references . 

As mentioned in earlier work the same pressure-time, and 
temperature -time history was used for all the propellant combinations 
but individual volume-time histories had to be worked out for each 
different propellant combination. 

Through a rather large iterative computer program such quantities 
as partial moles, partial pressures, partial volumes, partial weights, 
volume of air entrained, weight or air entrained, unburned fuel 
present, etc., were calculated. 

Some of these quantities were then normalized and one set is 
presented here graphically as a function of time. They are 


a. 

Fuel and oxidizer consumption wgt. 

(normalized) 

Fig. 

25 

b. 

Volume of entrained air 

It 

Fig. 

26 

c. 

Partial Pressures 

• t 

Fig. 

27 

d. 

Total Volumes 


Fig. 

28 

e. 

Partial Weights 

If 

Fig. 

29 
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These results and the manner of presentation, it is believed 
give a good picture of the composition and its time variation of the 
fireball and combustion products cloud. 

It is believed that the graphs are self-explanatory and the 
results can easily V* compared when different fuel-oxidizer 
combinations are used. 

The Fluorine tracer quantities added seem to have a hypergolic 
effect upon the cryogenic propellants to render the predictions of 
the most probable ignition times and delay times, and thus the yield, 
more certain. 
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FIGURE .75 . FUEL CONSUMPTION FOB LH 2 /RP-I/LO 2 + 10% F LIQUID PROPELLANT 

EXPLOSION ( YIELD - 4.5 PERCENT ) 
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FIGURE 26 . VOLUME OF ENTRAINED AIR FOR LRgVRP-l/LOg + 10% F LIQUID 

PROPELLANT EXPLOSION ( YIELD - 4.5 PERCENT ) 
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FI GURE 27 PARTIAL PRESSURES AND PARTIAL VOLUMES FOR LHg/RP-l/LOg + 10%F 

LIQUID PROPELLANT EXPLOSION PRODUCTS ( YIELD - 4.5 PERCENT ) 
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PART II 

EXPERIMEMAL VERIFICATION 

In this section the experimental work which was carried out in 
support and for the purpose of verifying the theoretical work will be 
discussed. The experiments verified the theoretical predictions of 
explosive yield and other characteristics of liquid rocket propellants. 

All the experiments can be divided into two large classes, labora- 
tory experiments and explosive field experiments. 

The laboratory experiments were carried out with inert components 
to reduce the hazards and only with small quantities of explosive 
materials so that explosions were not expected. 

The laboratory experiments can be further subdivided into Mixing 
Function development and electrostatic charge generation. 

The explosive field experiments were carried out to substantiate 
the information obtained by theory and then verified in the laboratory. 
The fie 1 d experiments were carried out with actual liquid rocket pro- 
pellants, and substantiated the laboratory data, and in addition pro- 
vided actual explosion information. They can be grouped as mixing and 
explosion experiments, electrostatic generation and auto ignition ex- 
periments. 

The present section will discuss in detail the experimental methods 
used and the results obtained. 

Laboratory Experiments 

The first large program carried out in the laboratory was the 
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development of methods by which the Mixing Function could be deter- 
mined, either for verification of the theoretical one or as a direct 
input into the Seven Chart Approach. 

Mixing Function Determination (i) 

It is always desirable to check any theoretical work with cual 
experimentation so in this section it will be shown that four methods 
have been developed which can give the Mixing Function experimentally. 
The Mixing Function is defined as the proportion or fraction of the 
propellants mixed at any time t . This fraction can be expressed in 
a number of ways as will be seen. It can be defined in terms of con- 
tact area or surfaces proportional to this area, in terms of liquid 
interfaces, of mixing surface profiles, etc. 

It will be shown how this was done in various manners by four 
experimental methods which are independent but complement each other. 
In many experiments two of these methods were used simultaneously so 
that the results could be compared for one and the same experiment. 

These four methods are: 

A. The Film Analysis 

B. The Wax Cast Analysis 

C. The Vibration Mixing Analysis 

D. The Thermocouple Grid Analysis 

The film analysis is a method which depends on high-speed pho- 
tographic coverage of the mixing phenomena, followed by an analysis 
and proper interpretation of these records. 

The Wax Cast Analysis "freezes" the mixing process at different 
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stages of the mixing so that the solidified image can be analyzed at 
a later date and at leisure. 

The Vibration Mixing Analysis simulates small fluid elements by 
different density, different color, different shaped particles such 
as marbles. By shaking them on a vibration table for a certain length 
of time, and by removing particles periodically, such things as evapor- 
ation or spill losses may be simulated. 

Thus the fixing phenomena can be studied in small incremental 
steps, using marbles or particles as fine as powder. 

The Thermocouple Grid Analysis, the best and most powerful of the 
four methods, consists of a three dimensional grid of fine thermo- 
couples, which produces a time record of temperature and its variation 
at many points in a region in which fluids at different temperatures 
are being mixed. These readings can be interpreted as will be explained 
to give insight into the phenomena leading up to and producing the 
explosive yield. 

In the following pages these four methods will be discussed in 
detail, and some of the results which have been obtained will be pre- 
sented. These results will give an indication of the value of these 
study methods in investigating the mixing phenomena of liquid rocket 
propellants. The results presented here are not intended to simulate 
a particular missile or missile failure (This will be done in the last 
two sections of this report) but rather to demonstrate how these methods 
can be employed in the study of the mixing produced by many different 
missile configurations and types of failures. 
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THE METHODS 
A. Film Analysis 

For the application of this method for analysis, it is necessary 
to obtain a high-speed pictorial record of the mixing phenomena, in- 
volving high-speed photographic equipment. For the present investigation 
transparent configurations were chosen and then, by placing mirrors 
at various angles, three-dimensional views were obtained on the same 
frame of the high-speed recording film. For non-transparent arrange- 
ments, X-ray, Gamma-ray, or tracer techniques could be used in a 
similar manner. Seme of the latter techniques have been used by one of 
the authors in concentration studies and profile interface studies of 
two- and three-phase mixtures*’ *5, 49, 75, 104. 

Figure 30 presents an overall view of the experimental apparatus 
used in the study of the mixing processes between hot or cold oil and 
water, hot wax and water, and kerosene, etc. It consists essen- 
tially of a Pyrex glass tube filled to a desired level with one liquid, 
and another Pyrex tube above filled with the desired amount of the 
other liquid. A diaphragm between the two glass tubes holds the upper 
liquid in place. At time zero the diaphragm is pulled out by a quick 
motion so as to remove it almost instantaneously. Slower removal of 
the diaphragm according to a programmed input can, if desired, simulate 
a progressive failure. If a complete diaphragm or bulkhead failure is 
not desired, a second fixed diaphragm can be inserted in addition to 
the removable one so that when the latter is removed, a desired size 
opening at a desired location remains, simulating a particular type of 
failure. 
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Figure 31 shows a number of frames for oil above and water 
(below) simulating bulkhead failure in a 1 1/2 in. I.D. tube and a 
ullage space of 4 1/4 in. Successive frames taken at 64 frames per 
second show the progression of the mixing process from three views, 

120 degrees apart. 

The analysis can be made directly from the frames as projected 
on a screen or from a more permanent record by either enlarging them 
on photographic paper or tracing the outlines of the mixing fronts as 
shown in Figure 32. 

Figure 32 further indicates the reproducibility of the mixing 
experiments by showing the traces of three identical runs and how 
close the total profile areas match. 

From the three-dimensional views or traces, the total surface 
areas were determined by first cutting the total volume into thin 
irregularly shaped discs and then graphically obtaining the peripheral 
surface of the discs. The resulting outside surface area, or A q , 
can then be plotted versus film frame or versus time (see Fig. 33). 

Since there are, however, droplets or particles of one liquid 
(vapor or solid) surrounded by the other within the space circumscribed 
by the profile, to get the true contact area between the liquids, 
these additional surfaces must be taken into account. This was ac- 
complished by a tedious method of counting these droplets, determining 
their size and surface areas, and adding these new areas to the profile 
areas. The ratio of the total area to the profile or outside area 
of any particular frame was defined as the "Turbulence Factor" for 
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this frame. This slow process of determining the total areas or 
;....ct areas was shortened considerably later by the very satis- 
factory approximate method of comparing frames for different runs 
with standard frames for which the "Turbulence Factors" had been de- 
termined carefully, and ascribing the same "TVirbulence Factor" to 
otter similarly appearing frames. 

Figure 34 presents the "TVirbulence Factor" for the runs presented 
in Figure 33 . 

Using Figure 33 and Figure 34, and combining than, gives the 
total area or an area directly proportional to the true contact area 
between the liquids. The result is seen in Figure 35. 

The "TVirbulence Factor" has been further substantiated by the 
Wax Cast Analysis which will be described below. 

Figure 35 , in addition to giving the total area, shows the re- 
markable reproducibility of these experiments. Furthermore, this 
figure demonstrates the effect of mixing energy. It shows three runs 
with a 2 1/4 in. ullage space and three runs with a 4 1/4 in. ullage 
space. Considerably more area is obtained for the 4 1/4 in. ullage 
space since the upper liquid obtained greater kinetic energy before 
mixing. In these, as well as all other experiments, it was observed 
that the reproducibility of mixing increases with the increas. in 
mixing energy. 

The Film Analysis as described above proved to be an invaluable 
aid in the study of the mixing phenomena of liquids. The results, 
however, were not taken as correct unless they checked with the results 
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Figure 35 Contact Area versus Time (Film Analysis) 
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of another of the methods described below, simultaneously applied to 
the same experiment. 

In the case of the mixing of hot and cold liquids with tempera- 
ture differences great enough to produce boiling of one of the liquids, 
another factor - the "Boiling Factor" - was introduced which accounted 
for the surface produced by the vapor bubbles of one of the consti- 
tuents. 

A number of different diameter Pyrex glass tubes were used in the 
mixing studies ranging from 3/4 in. I.D. to 6 3/4 in. I.D. to give a 
relationship for size. In these experiments the effect of surface 
tension could be observed since it altered some of the phenomena in the 
smallest sizes. 

B. Wax Cast Analysis 

The Film Analysis as described previously was a dynamic method 
of analysis and the results had to be obtained by interpreting recor- 
dings of the actual mixing phenomena through the profile area and 
utilizing the turbulence, boiling, and freezing factors to obtain the 
total area or contact area. 

The Wax Cast Analysis allowed the "freezing" of the actual mixing 
process at various stages of the mixing by mixing hot wax with cold 
water. By varying the temperatures slightly the wax would solidify 
earlier or later in the mixing process and the wax casts obtained in 
this manner could br' analyzed at leisure at any time thereafter. 

The experimental apparatus and setup for the Wax Cast Analysis 
is identical with the one for the Film Analysis, and so for many 
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experiments both methods of analysis were used simultaneously and 
checked against each other. 

Figure 3 6 shows a series of such wax casts representing a com- 
plete experiment. Excellent reproducibility indicated by the Film 
Analysis allowed the experiments to be carried out identically except 
for slight variations in temperature. The wax casts indicate the 
progress of the mixing process. The profile area was determined from 
these wax casts by two methods: A small grid was laid out on the 
surface and then the area determined by counting squares. This pro- 
cedure was checked by dipping these samples into paint and letting 
them dry. By determining the thickness of the paint film and its 
density and by weighing the sample before coating and after coating, 
the profile area could be detemined. Essentially the same answers 
were obtained by both methods. 

The serially sectioned samples or wax casts in Figure 37 allowed 
the total area to be determined by the same method used for the pro- 
file area again represented the "Turbulence Factor". 

Curves of the profile area versus time (Figure 38), the total or 
contact area versus time (Figure 39), and the "Turbulence Factor" 
versus time (Figure 40) were plotted. 

The Wax Cast Analysis was primarily used to check the results 
from the Film Analysis and, after it was found that the results from 
both analyses were in essential agreement, it was not used further. 

C. Vibration Mixing Analysis 

Another method by which the mixing processes and phenomena were 
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studied is the Vibration Mixing Analysis. In this method, the fluid 
particles are simulated by solid particles of various sizes, shapes, 
colors, densities, etc. These solid particles used in various pro- 
portions and configurations are mounted on a vibration table and 
shaken for certain lengths of time. Figure 4j shows the experimental 
arrangement simulating spill mixing configuration. Three liquids are 
represented by different color marbles and are arranged in a desired 
configuration. Shaking the configuration for predetermined times, the 
marbles will diffuse into each other. Thus in this manner the mixing 
process simulated here can be studied in steps. 

Figure 42 presorts a number of the views representing these steps . 
Any particular region can be studied in this manner by removing this 
volume and counting the particles of each of the constituents present. 
This gives the percentages of each one of these constituents, thus the 
degree of mixing in this region. 

Removal of some of the particles from certain regions as programned 
and governed by the fundamental relationships of heat transfer can 
further account for evaporation losses. This takes into account boil- 
off as part of the mixing process. 

In this manner a curve can be plotted for the fraction mixed as a 
function of time. The time scale is arbitrary since the amplitude 
and frequency of the shaking table have a pronounced effect on the 
speed with which the mixing progresses. To find the absolute mixing 
time either an experiment with liquids or theoretical calculations 
will give the time scale needed. This absolute time superimposed upon 
the curve gives the true mixing function - time relationship. 
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Figure 43 presents a mixing function curve representing the 
J-test series of the ' v “thur D. Little Spill Test Program 6 . 

It is interesting to note that all mixing experiments produced 
the same characteristic shape of the mixing function curve, only the 
actual values and speed with which the process occurred were different. 

D. Theiaocouple Grid Analysis 

The three methods for studying the mixing process, the Film 
Analysis, the Wax Cast Analysis, and the Vibration Mixing Analysis are 
excellent methods for obtaining insight into the phenomena taking place 
when different liquids mix. They are, however, simulation methods 
and are not readily applicable to actual explosive mixtures where 
detonation and explosion would, in most cases, destroy the records. 

The Thermocouple Grid Analysis overcomes this difficulty and can, 
therefore, be considered the most powerful and best of the methods 
discussed in this paper. It is the best but also the most elaborate, 
most expensive (and the data reduction the most tedious) of the four 
methods. For explosive tests, it is, however, the only method which 
will give information from time zero of the failure up to and beyond 
the time of ignition. 

The heart of this method is a grid of thermocouple junctions spaced 
throughout the region under study. These can be placed inside the 
tanks of a missile and extend around it if fireball data are desired. 

A continuous time record of the signals from the individual junctions 
is kept for the complete period of the processes under investigation. 
Thi. method, through its three-dimensional records of the complete 
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time history, can provide a great amount of information with respect to: 

1. The three-dimensional mixing front of n particular constituent 

2. The degree of mixing at a particular point 

3. The degree of turbulence at a particular point 

and, in explosive tests, in addition to the above: 

4. The location of the point or points of ignition 

5. The time delay from initiation to mixing to ignition 

6. The propagation of the reaction front 

7. The propagation of the shock front 

8. The separation of shock front and reaction front, etc. 

9. Other phenomena and events obtainable by detail data analysis 

The experimental arrangement and equipment for the Thermocouple 
Grid Analysis is the same as that for the Film Analysis and can readily 
be incorporated into static explosive test series. The only addition 
necessary is a three-dimensional grid of fine thermocouple junctions 
giving good response characteristics. Rill response times of less 
than 10 micro-seconds have been obtained in our laboratory. 

The signals from these thermocouple junctions can be fed to the 
recording equipment, which may be close by, in inert tests, or at some 
distance in explosive tests. 

The overall experimental arrangement for the work reported here 
is shown in Figure 44. The experimental apparatus and the ecording 
equipment are shown. All the control experiments using this arrangement , 
in addition, made use of high-speed camera recordings of the mixing 
phenomena so as to have a check and comparison between the Film Analysis 
and the Thermocouple Grid Analysis. 

Figure 45 pictures thermocouple grids which were used in rome of 
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these investigations. Our laboratory has the capability of monitoring 
over 40 individual junctions at present which can be extended to 65 
if needed. By high-speed periodic sampling through commutation this 
capacity can be increased manyfold, but the time continuity of the 
records has to be sacrificed. 

Figure 46 presents traces obtained from a mixing experiment. All 
twelve traces in this case correspond to junctions in a vertical plane 
at three different elevations with four junctions equally spaced in 
each one of these elevations and arranged in straight lines. 

Figure 46-A presents traces of the mixing of hot oil and water. 
Figure 46 B presents traces of the mixing of LN^ and kerosene. 

Figure 47a presents the resulting mixing profile as determined by 
interpolation between the junctions of the thermocouple grid at time x . 

Figure 47b also presents the resulting mixing profile of the same 
hot oil and water experiments as determined by the Film Analysis at the 
same time x . 

It can be seen that the results from both methods of analysis are 
essentially the same. More and closer spaced thermocouple junctions 
would give more points for drawing the mixing profile and would fix 
this surface with greater reliability. 

The degree of mixing at any time, x , around a particular junction 
can be determined from the time history at this j miction by writing 
the relationship for mass and energy balance at this junction incor- 
porating the laws of thermodynamics, fluid flow, and heat transfer. 

A Computer program to do this is of great help if a considerable number 
of junctions c.re involved. 
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a. Thermocouple Grid Analysis b. Film Analysis 

Figure 47 Comparison of Mixing Profiles Obtained 
Simultaneously by Thermocouple Grid 
Analysis and Film Analysis 



Time in Seconds 

Figure 48 Mixing Function or Spill Function for Simulated 
Two Component Liquid Propellant Mixing Test 
(Thermocouple Grid Analysis) 
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Averaging the conditions around each of the individual junctions 
can then give the mixing function, the fraction mixed at any time t , 
Figure 48. 


SUNMARY 

The preceding discussion of the four methods which can be employed 
investigating the mixing phenomena of liquids to establish a mixing 
function- time relationship shows the value of these methods, their 
relative merits, and singular advantages to do the job. 

The Film Analysis is relatively simple and easily carried out, but 
in its simplest form, the use of light requires transparent containers 
and rather transparent constituents. In its more sophisticated forms - 
using shorter wavelength radiation, such as X-ray, Gamma -ray or tracer 
methods - the equipment necessary becomes ouch more complex. 

The Wax Cast Analysis allows the checking of the results obtained 
by the Film Analysis by independently establishing the mixing profile, 
total or contact area, and turbulence factors. 

The Vibration Mixing Analysis provides a method by which, again as 
in the Wax Cast Analysis, the mixing process can be stopped at any point 
in its development and studied at leisure. 

- The Thermocouple Grid Analysis, the most powerful of the four 
methods, can provide all the information of the above methods (except 
for stopping the mixing process at any point in its development) and, 
in addition, can be used in actual missile configuration explosive ex- 
periments. 
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In the maimer described above these methods can provide information 
which is not available at the present time and, through this information, 
can give better insight into, and understanding of, the actual happenings 
during the mixing process and leading up to the explosion of liquid 
propellant mixtures - information which is needed if we ever hope to 
control and guide these processes. 


Electrostatic Charge and Voltage Generation (2) 

Mich of this method was described in the section on Critical Mass 
and some of the supporting results were given at that time. 

Basically it was theoretically predicted when fluids mix that 
electrostatic charges and voltages are generated and these v re then 
verified by actual measurements in the laboratory. 

To verify the results obtained by use of the Fluid Plug Model a 
simple experiment was set up. Figure 49 schematically describes it. 

A glass cylinder was used and filled to a desired level with RP-1. 
Above the RP-1 was a partition with a hole in the center which could 
be opened by sliding out two plastic sheets. The hole diameter was 
1/2 that of the glass cylinder. The space above the partition was 
filled with the desired quantity of LN 2 and when the hole was quickly 
opened the IJ^ dropped into the RP-1 through the hole, simulating a 
fluid plug. 

In the RP-1 space were mounted two screen electrodes, one just 
above the RP-1 liquid surface and the other in this case 1.25 inches 
below the upper one. Leads from the two screen electrodes were led 
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to an Electrometer and its output to a Brush strip chart recorder. 

In this manner, the charge and voltage build-up as a function of 
time could be determined for different combinations of liquids and for 
different quantities. 

The actual traces obtained look very much like the theoretical 
traces of Figure 11. Again it does not seem to make any difference 
■whether the whole quantity is used all at once or in steps, providing 
in the latter case leakage is prevented. 

First single probes were u* 2 d and only on occasion were charges 
and voltages observed, then double probes were used with similar 
results. The more probes were used the more often were readings ob- 
tained and finally with the use of screens which are really just multi- 
junction probles, readings were obtained every single time. It was 
observed that the screens had an averaging effect. It could be es- 
tablished at lerst that the charge and voltage generation is an integral 
part of the mixing process of liquids and occurs whenever there is 
mixing of dielectric fluids. 

When the large quantities were allowed to come together it seemed 
that the mixing takes place in very rapidly repeated jumps so that the 
same voltage was reached as if the same quantities were mixed in dis- 
crete steps. Since the mixing occurs in such small steps over and over 
again, the average value takes on much more meaning especially when 
applied to large quantities. 

The manner o ixing did not seem too critical since just pouring 
one component into she other produced essentially the same results. 
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The reason for this seems to be that the boiling action which provides 
considerable violence to the mixing process is the over-riding pheno- 
menon. If, however, considerably more mixing energy is used, such as 
the pouring from a height of 10 feet, the mixing and the charge gener- 
ation and voltage generation occur much more rapidly. 

It seemed, however, that the mixing occurred relatively more 
rapidly than the charge and voltage build-up then compared with the more 
gentle method of bringing the components together. Figure 5 0 shows an 
actual trace of voltage versus time for 200 ml incremental additions. 

Due to the small quantity of RP-1 and the rapid addition of 
the temperature difference between the components decreased during this 
test showing the effect of less vapor generation in the later additions, 
resulting in deeper penetration of the fluid plug and therefore more 
voltage generation. 

Hie high speed film analysis indicated that the bulk of the bubbles 
was very uniform in size and very close to 1/4 in. in diameter. There 
were larger and smaller ones but they were very few in number. 

Sometimes voltages of opposite polarity are recorded since with 
the turbu 1 ence in the mixing region the charged layers can come in 
contact with either electrode, and each one can became positive or 
negative with respect to the other. The actual amounts of voltage gen- 
erated and recorded did not in general vary too much but on occasion 
a large jump was observed. 

The average voltage generated was four volts with electrode spacing 
optimum based upon many observations for 200 ml of LJ^ • This voltage 
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Fig. 50 Typical VOltage-Time Trace for Successive Mixing of LN 
and RP-1 Using Small Quantities of LN, 
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was built up through the use of large quantities or small quantities 
in succession to several hundred volts during these experiments . Some- 
times, but not very often, a single jump of several hundred volts was 
observed with as little as 200 ml of I^. 

Hie actual values of charges and voltages generated increased with 
the quantity of propellants mixed. 

The theory and the laboratory verification resulted in the fonmi- 
lation and prediction of the existence of the critical mass which was 
several years later further substantiated in actual field tests. 

Another important observation in these experiments was that the 
amount of charge and voltage generated depended upon the energy which 
was used in mixing the components, the results of this were presented 
earlier. 

B. Explosive Field Experiments 

To determine how mixing of liquid rocket propellants compares with 
the laboratory studies which have substantiated the theories, a number 
of different experiments have been conducted bring out this comparison. 

1. FYRO Project 

25,000 lb. LCU/RP-l Explosion Experiments 

To show and demonstrate how the theory can predict the explosive 
yield for liquid rocket propellant explosions, two 25,000 lb. LO^/RP-l 
experiments were completely instrumented by the University of Florida 
Group with thermocouple grids. Hie thermocouples were installed inside 
the tank assemblies of the liauid propellant explosion test series, 
planned and conducted with project PYR0 at the Air Force Rocket 



Propulsion Laboratory at Edwards Axr Force Base, California. This 
instrumentation was to measure the phenomena following the initiation 
of failure. 

Because of the mode of failure selected by project PYR0, the 
region swept through by the star cutter had to be excluded from detailed 
analysis since the instrunentation in that region would have been des- 
troyed by the cutter (see Figure 51 ) before the events which were to 
be measured occurred. 

The overall purpose of this instrumentation was to: 

1. Correlate the mixing phenomena of the true propellants 
with laboratory mixing experiments using inert fluids for 
simulation. 

2. To substantiate experimental ly parts or all of the "Fire- 
ball Hypothesis"^ proposed earlier in these studies. 

Specifically it was hoped to be able to determine by this experi- 
mental procedure part or all of the following: 

After failure but before ignition: 

1. The three-dimensional mixing front, or boundary of the mixing 
region 

2. The degree of mixing at a particular point 

3. The degree of turbulence at a particular point 

4. The time delays from initiation of failure to start of mixing, 
and to ignition 

5. The location of the point or points of ignition 

6. The propagation of the reaction front 

7. The propagation of the shock front 
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8. The separation of the shock front and the reaction front 

9. Other phenomena and events obtainable by detailed data analysis 

Only in (2) and (3) above do the thermocouple response charac- 
teristics have to be considered since in all other cases only relative 
time differences arc needed. 

Excellent data were obtained in both experiments. Advantage was, 
however, taken of knowledge obtained from the data analysis of the first 
experiment. No. 278, to obtain the best results possible from the second 
experiment, No. 282. The main improvements were the moving up of sane 
of the thermocouples higher in the tank and closer to the star cutter, 
or into the mixing region, since it was found that the star cutter did 
not travel as far as was previously expected. The chart speed of the 
recording oscillographs was, in the second experiment, operated at four 
times the speed of the first experiment to increase the resolution 
capabilities. 

The instrumentation to accomplish the above consisted of: 

1. Very fine thermocouples inside the tank 

2. A reference junction box in an underground steel box near 
ground zero 

3. High speed recording oscillographs 

The sensing elements inside the tank (Figure 52 ) consisted of 
37 copper -constantan thermocouples , No. 36 gage, shellacked and Teflon- 
coated. These thermocouples were made so that no excess material was 
present at the junction and the whole thermocouple looked like a con- 
tinuous wire with no visible variation at the junction (Figure 53) . 
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A small plastic tube was slipped over the wires to about 1/4 inch 
behind the junction and filled with epoxy, both to support the junction 
preventing relative motion between the wires, and to provide a means 
for attaching the thermocouple to the support wires inside the tank. 

The thermocouples were located in the tank (Figure 54) so as to 
give an optimum pattern of the data. They indicated the arrival of 
events at their location and changes occurring at or in their region. 

The leads of the 37 thermocouples (74 individual wires) were 
guided along the support wires and loosely tied to them about every 8 
inches (Figure 55) . This gave the needed support to the fine wires 
and at the same time allowed than to give, in case they were hit by 
some small fragments from the shattered glass diaphragm. 

The leads were then fed by means of two Cannon plugs through the 
tank wall, and on the outside by copper -constantan thermocouple cables 
to the hot reference junctions, located in an underground steel boot. 
From there the signals produced by the thermocouples were carried by 
copper cables to three CEC recording oscillographs. These recorders 
were operating at a chart speed of 40 in. /sec for experiment 278 and 
at 160 in. /sec for experiment 282. 

The data obtained by the above instrumentation and methods are 
shown in Figure 56 for a few theimocouples ; Figure 56A giving traces 
from initiation of failure to some time after detonation; Figure 56B 
showing only the time increment during which mixing starts and indi- 
cating the passage of the reaction and shock fronts. 

Analysis of the data from these two liquid propellant explosion 
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experiments allow the following statements to be made: 

1. Experiments No. 278 and 282 were amazingly similar in propellant 
mixing, ignition and explosion characteristics. 

2. Practically all the mixing up to the time of ignition was confined 
to the volume swept through by the star cutter. This volume as 
determined from the data was about 12 percent of the total, based 
upon the RP, giving an ur r> - yield estimate of about 12 percent of 
the theoretical maximum or somewhat less than 15 percent TNT 
equivalent. 

3. The actual temperatures as recorded by the traces allow the calcu- 
lation of the degree of mixing by means of a number of simultaneous 
equations which can most conveniently be solved by computer . 

4. The degree of turbulence at a particular point (the points where the 
thermocouple junctions were located) can be determined from the 
fluctuations and fine structure of the traces in terms of frequency 
and amplitude. 

5. There was only one point of ignition in each of the experiments and 
this point was located in each case (Figure 57A and 57B) . 

6. The time delays between various events can be determined from the 
traces. If projected o t a screen these times can be read to a 
1/1000 of a millisecond. The time delay from failure or firing of 
of the ram to ignition in experiment 278 was 543 milliseconds and 
in experiment 282 it was 580 milliseconds (read to the closest 


millisecond. Greater accuracy is available if desired). 

The accurate reading of these delay times also allows the exact 
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determination of the film speeds of the various cameras used. The 
importance of this is brought out, since none of the cameras operate 
at their nominal speeds and vary from one time to the next. For 
instance, the two 4000 fps cameras during experiment 282 were oper- 
ating actually at 2620 fps and 3750 fps respectively. 

7. The propagation velocities of the reaction front and shock front 
could be determined and the analysis gives the results shown in 
Figure 58. 

8. By careful study of the fine structure of the traces the shock 
front propagation characteristics can be deteimined. 

9. From f7 , (8) above it can be determined where the two fronts 

travel together and where and when the shock front leaves the re- 
action front behind. 

10. Witu the delay times exactly determined and the film speeds known, 
then the volume-time or volume-distance characteristics of reac- 
and shock frontr. can be determined (Figure 59) . 

11. The fine structure of the traces reveals that some of the glass 
fragments from the shattered diaphragm hit some of the thermocouples 
in the uppermost of the four layers but did not damage than, so they 
recorded this and later events. rp he glass fragments did not pen- 
etrate during their high energy state to the lower layers since the 
attenuation in the liquid was too great. 

12. From the fine structure of the traces, at least three typical shapes 
are of .erved, Figure 60. The first is caused by the arrival of the 
LC "ol lowed by the reaction jnd shock fronts; the second is caused 
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by the reaction and shock fronts arriving simultaneously at the 
thermocouple , and the third is caused by the shock front arriving 
first closely followed by the reaction front. 

The above statements are only a few of the many which could be 
made in extracting information from these experiments . 

Simulating the configuration and the mode of failure of the 25,000 
lb. LD^RP experiments in the laboratory produced the curve shown in 
Figure 61. It again shows the typical quick rise and then the two 
typical limps for this type of failure. Naturally in the actual case 
the curve after the ignition point did not exist. 

The ignition point was determined by the critical mass technique 
which indicated that when 2900 lb. of the propellants are mixed, ignition 
is a certainty. Thus both the quantity and the ignition time were pin- 
pointed and then plotted on the mixing curve. Again perfect agreement 
is found between the experiments in the laboratory and the results ob- 
tained in the actual explosive field tests. 

Taking 2900 lb. and dividing it by 25,000 lb. or the total quantity 
of propellants present gives an indication or prediction of explosive 
yield of 11.6% . This value can be compared to the value which was 
measured, namely 121. 

200 lb. Cold Flow Experiment 

To actually compare the mixing characteristics of the inert fluids 
used in the laboratory for simulation and the explosive liquid rocket 
propellants a 200 lb. LO^/RP experiment was completely instrumented in 
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the same manner as the 25,000 lb. experiments. 

An explosion was not expected in this case and so the mixing 
process could be followed for a much longer period, as a matter of fact 
until the mixing actually stopped due to sludge formation. 

Again the cutter volume had to be excluded from the analysis but 
the results were in excellent agreement with the curves obtained for 
the simulation. Figure 62. 

The mixture was then set off by an explosive bolt and the yield 
actually measured was the same as predicted. 

The S-IV B Explosion Experiment 

To see how a prediction based on the simulation of the S-IV B 
explosion experiment would compare with the actual measurements, the 
mixing function was determined for the above configuration and the 
corresponding mode of failure. Figure 63. 

The ignition time was determined again as in the previous predic- 
tions by the critical mass method and thus was plotted on the mixing 
function curve which is the yield curve in the early stages. 

From the simulation it was shown that the mixing occurred with two 
distinct larger humps but which even at their peak values are relatively 
low and then stabilizes %x> an almost constant value. This is due to the 
evaporation occurring essentially at the same rate as new propellants 
are added to the mixture. 

The explosion occurred again very early thus not allowing large 
values of yield to be obtained. 

This self-limiting phenomenon is fortunate since large explosions 
with liquid rocket propellants are almost impossible to obtain. 




Time in Seconds 

Figure 62 Expected Yield as a Time Function 

(200 lb Cold Flow and Explosion Experiment) 
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Figure 63 Mixing Function, S-IV LO.^/LH^ (Based upon 3" Diameter 
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The ADL Spill Explosion Experiments 

The work done with liquid rocket propellants by Arthur D. Little, 
Inc.^ was to determine the explosive yield obtainable with liquid 
rocket propellants. The experiments were designed to give very good 
and rapid mixing so as to obtain large yields. 

When small quantities of propellants were used the ignition time 
could be controlled and was varied. It was observed tiiat early ignition 
produced small explosive yields, and late ignition produced small 
explosive yields. In between a maximum value was observed. 

When large quantities were used in the experiments early ignition 
always occurred due to self-ignition and the yields obtained could not 
be controlled. 

Again the mixing function was obtained in our laboratory and with 
the critical mass method, considering the higher mixing energy available, 
and the ignition time was determined. This time could also be checked 
against the measured value, Figure 64. 

Again excellent agreement between the predictions and the actually 
measured values could be observed. 

Electrostatic Charge Generation and Auto-Ignition Results of Liquid 
Rocket Propellant Experiments 

It has been observed that when large quantities of liquid fuels 
and oxidizers are brought together, either during experiments or 
accidentally, liquid propellant explosions result. 

Many different phenomena can trigger these explosions such as 
flames, sparks produced by striking or breaking metal, hot materials or 
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hot spots produced by slow chemical reactions of fuel and oxidizer, 
the breaking of crystals which are formed when one of the liquid com- 
ponents freezes the other and which are broken mechanically or by 
thermal stresses, or by static electricity which is the result of 
internal friction and which may produce a spark discharge. 

Many more possibilities could be cited but these additional sources 
do not change the basic picture. 

Thus if there is contact or mixing of a liquid fuel with a liquid 
oxidizer, ignition is possible if an ignition source is available. If 
not, the mixing process may proceed with more and more of the fuel and 
oxidizer mixing until an ignition source appears either through external 
or internal action. 

The Critical Mass concept was developed earlier^’^’^^’^ as 
also shown in detail in section I-C. Figure 65. 

2. KSC Auto-Ignition Experiments 

The following presents verification of the theory and predictions 
which were made ear ii er 1 »2»3»4,24,75,79,104,116,118 experin.-n 

tation and measurement of the critical parameters involved in the real 
explosive Liquid Rocket Propellant combinations. 

For the actual explosive experiments which w -• planned and then 
carried out at the Kennedy Space Center, the experimental procedure 
had to be modified as well as the equipment . Since the experiments 
had the potential of an explosion (although none or only a very few 
were expected with this quantity) precautions requiring remote oper- 
ation had to be taken. 



Explosion (Mixed) Mass, lb 
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Fig. 65 Explosion Mass - Mixing Energy Relationship 
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A test set-up was designed, consisting essentially of two cylin- 
drical containers, a larger one on the bottom and a smaller one on top. 
The exact dimensions were. Lower Tank: 12 3/4" O.D.X. 12 3/4" high, 

Upper Tank: 8 3/4" O.D.X. 8 3/4" high. 

The top tank was hinged so it could tip and pour its contents into 
the lower tank. See Figure 66. Both tanks had lids and were thermally 
insulated so as to reduce the evaporation losses of the cryogenic 
liquids used as fuels and oxidizers. 

Temperature probes, consisting of stainless steel sheathed thermo- 
couples, two in each of the tanks, were used; these tonitored the liquid 
levels in the tanks so that proper quantities were in each tank at the 
time an experiment was started. 

The lower tank was fitted with two copper screens, approximately 
one inch apart, spaced so that they would be surrounded by the mixing 
region during the experiment. Thes v screens formed the electrodes to 
measure the charges and voltages generated due to *he mixing processes. 
The signals from the probes arid screens were transmitted through am- 
plifiers to the instrumentation van where they were recorded on strip 
charts and tape, figure 67 shows the screens installed in the lower 
tarJc. 

The test set-up allowed experimentation with any one of the pro- 
pel la. ts in either tank. 

The experimental procedure, described in greater detail in the 
appropriate section, was to fill each tank to the desired level, then 
by remote operation: 




do unci: 



Lon in Lower Tank 

m 
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1. Lift the lid of the lower tank 

2. Poll the safety pin cm the upper tank 

3. Dump the upper tank 

Records of the probe temperatures and screen voltages were made 
Airing the above procedure. 

The 60 Pound Explosive Experiments 

The test set-up for the sixty pound experiments was essentially 
the same as for the six pound experiments except that the tanks needed 
were considerably larger. The dimensions were. Lower Tank: 27 1/4" 0.0. x 
27 1/4" high. Upper Tank: 18 1/2" O.D. x 18 1/2" high. 

Figure 68 gives the experimental arrangement for this series of 
experiments with the instrumentation part being the same as for the 
smaller tests except that the stainless steel probes which had proved 
to be erratic and slow in response were replaced with exposed small 
thermocouples and the screens were larger in diameter. 

The actual operation and running of the experiments was slower 
since larger quantities of cryogenic fuels and oxidizers were involved 
requiring longer chill-down and filling times as well as longer warm-up 
times after the experiments were conducted . 

The 240 Pound Explosive Experiments 

While no actual explosions were expected with the six pound ex- 
periments , maybe one or two with the 60 pound experiments (although 
none occurred for the number of experiments carried out) , about one 
out of ten experiments with the 240 pound quantities was expected to 
produce an explosion. 






175 


The experimental ret -up was essentially the same. See Figure 69 . 
The tank sizes were. Lower Tank: 43 1/4" O.D. x 43 1/4" high. Upper 
Tank: 30 1/4" C.D. x 30 1/4" high. The exposed thermocouple type probes 
used in the 60 lb. test had not always been able to distinguish be- 
tween the vapor and the liquid; therefore, these were replaced with 
heated thermocouples . These would remain cool in the liquid but would 
heat up rapidly in vapor where the heat gain was considerably larger 
than heat loss. Also the times for chilling -down, filling, stabilizing, 
dumping, and then warm-up after each experiment were considerably longer. 

Fluidic Oscillator Gas Sensors 

As an additional experiment not directly connected with the mea- 
surement of the electrostatic charge and voltage generation, fluidic 
oscillators, developed under a separate phase of this contract were 
used to monitor the cloud produced by the experiments . 

Figure 70 shows one of the earlier designs of this very simple 
device which can instantaneously detect a change in the density of a 
gas going through it, indicating the concentration of another gas. 

The device is simple, easy to make, and inexpensive . A transistor 
radio earphone can be used as the transducer to change the frequency 
impulses from the oscillator into an electrical signal which can be 
heard, displayed on an oscilloscope and/or recorded on tape. 

The oscillator frequency produced by the gas flowing through the 
device by applying a vacuum indicates the gas mixture and its concentra- 
tion. A separate report is being written on this part of the work. 
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Experimental Procedure 

General : The experimental procedure was designed to allow the 
filling of the two tanks of the experimental set-up to the required 
liquid levels and then allow the dunping of one of the components into 
the other with the resulting voltage and charge generation being measured. 

To accomplish this a test procedure had to be written since many 
people from different support groups were involved at Kennedy Space 
Center; this procedure had to be agreed upon and approved by everyone 
concerned. 

Evaporation or Boil -Off Experiments 

The first step in the procedure was to determine the heat losses 
from the two tanks so that the filling and holding procedure and timing 
could be worked out. 

Theoretical calculations were made on the evaporation rates of the 
cryogenic components and they were experimentally verified. The probes, 
two in each tank were set a certain distance apart and then the tanks 
were filled to the upper probe and the time determined for the liquid 
level to drop to the lower probe. 

This was done for IN-* IX^* and 11 ^ for each of the tanks in the 
6, 60, and 240 pound experimental configurations. 

Actually these boil -off experiments were carried out at the 
beginning of each of the series involving different quantities. 

Charge Generation and Auto-Ignition Experiments 

On the morning of each test day the calibration of the instrumen- 
tation was verified and the equipment was checked to see that it was 
recording properly. 
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1. to open the lower tank lid 

2. to pull the safety pin of the upper tank 

3. to dunp the upper tank 

The high speed camera and the documentary camera were started 
during the countdown. The recording instrumentation which was running 
all through the filling operation at low speed was switched to high 
speed during the countdown to produce maximum resolution for the dcta 
taken during the actual experiments. Figure 71. 

After completion of the test it was necessary to wait until the 
cryogenics vaporized before one could safely approach the test equip- 
ment, check it over, and ready it for the next experiment. 

This procedure was followed in each case for all the charge and 
voltage generation auto -ignition experiments. 

In all cases the LO^/U^ experiments were carried out first since 
both the LC >2 and the IH- vaporized without leaving traces of their 
presence and no cleaning operation was necessary. When RP-1 was used 
some of it froze and the waiting period until a-1 the LC^ vaporized 
and RP-1 gel melted was considerably longer. It also was necessary to 
remove any RP-1 which splashed into the upper tank which was to be 
filled with LC^ for the next test. 

Data and Results 

In this section e data and the results of the various experi- 
ments will be briefly reported and the significance and deeper meaning 
of these results will be discussed in the next section. 
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The documentary and high speed remote controlled cameras were 
loaded and readied for the experiments. 

With these tasks completed the chill -down of the tanks was 
started. was used first to conserve LC^ and 11^ . The tanks were 
both partially filled with LN^ and this then allowed to vaporize; 
this pre- cooled the equipment for the experiments. 

Regardless of whether the was to be in the upper or lower 
tank for the particular experiment the flow was started and the 
liquid level monitored. It took considerable time before the probes 
showed liquid levels. The lower probe was set at the level required 
for the experiment and the upper probe at the level to which the tank 
had to be filled to offset evaporation losses and have the correct 
quantity of liquid in the tank at the time of dumping. 

When the liquid level reached the upper probe the signal was given 
to terminate flow; then the LH, tanker truck was shut off, disconnected 
and driven away. This operation took from two to three minutes. 

When the tanker had cleared the hazard zone the signal was 
given by the safety director to start the LO 2 flow into the second 
tank. Again the liquid levels were monitored and this tank was filled 
to the upper probe and held at this level until the relative timing 
between the two tank levels was right so that with the proper count- 
down the liquid levels in both tanks were at the correct level at the 
moment of dump. 

At the command "dump", three control cables were pulled in se- 


quence: 
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Evaporation and Boil -off Experiments were necessary to predict 
the proper timing of the experiments . Evaporation losses occ»’i 1 in 
both tanks when filled with cryogenics; safety cansiderati reqiu <■'. 
these tanks to be filled sequentially. For this reason, theoretical 
calculations were made as to the rate at which the liquid leveJs 
dropped in the tanks and were then verified by experiment. 

The results obtained were as follows: 


Six Pound Experiments 


LN- 


L0- 


LH- 


60 Pound Experiments 
LN. 


10 - 




240 Pound Experiments 
LN_ 


10 - 


IH. 


Calculated Rates 


Upper Tank 
0.131 in/min 
0.08 
1.05 

0.113 

0.07 

.86 

0.084 

0.053 

0.67 


Lower Tank 
0.11 in/min 
0.067 
1.00 

0.089 

0.055 

0.72 

0.08 

0.05 

0.06 


Three different types of probes were used during the experiments. 
Each successive type was an improvement. Excellent agreement was 
shown between the experimentation and calculated values reported here. 
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Charge Generation and Auto- Ignition Experiment s 

The objective of this series of tests was to determine if electro- 
static charges are generated and if so to find their functional 
relationship to other parameters. An important parameter involved in 
these experiments is the weight or mass of the propellants involved. 

Six Pound Experiments 

Figure 72 shows the actual traces from a typical experiment. 
Recorded on the graph are the upper and lower probe readings in each 
of the tanks and the generated voltage. Table II lists all the 
LO 2 /IH 2 experiments. It should be noted that the actual voltage gen- 
erated exceeded the range of the instrumentation in the first five 
experiments. In the first test the range was + 100 volts, for the 
next four + 1000 volts and for the remaining experiments +_ 10,000 volts. 

For the LC^/RP-l experiments the range was set at + 10,000 
and after the first eight experiments changed to + 1,000 volts to 
obtain better resolution. Table III presents the LC^/RP-l results. 

60 Pound Experiments 

Figure 73 presents the actual traces of a typical experiment, again 
showing the temperature traces of the liquid level probes in both the 
upper and the lower tanks and the voltage generation trace indicating 
the "spike" voltage due to the mixing. Table IV presents the results 
for the LO 2 /IH 2 experiments in which the LO 2 was poured into the U^. 

The voltage range was set at + 20,000 volts. Table V presents the 
results for the LC^/RP - ! experiments with the voltage range set at 
♦ 1000 volts. The voltage range was then changed to + 10,000 volts 



Volts Lower Tank Upper Tank 
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Fig. 72 Chart Traces from the Six Pound Experiment 


Voltage Lower Tank Upper Tank 

Generated LH-> L0~ 
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TABLE II 

Voltage Data froa Six Pound LOj/U^ Experiments 
(LO-y is Poured into LHol 


Wo. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


Date Screen Diaaeter-In. Voltage - Volts 

1 / 12/72 U + 40 

>+100 

1 / 13/72 >+1000 

>+1000 
+ 350 
>+1000 

1 / 14/72 +4800 

+2800 

+1600 

+6400 


11 


+2000 
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TABUS III 

Voltage Data frog 6 lb, LO? — RP- 1 Experiments 
(IX ) 2 is poured into RP- 1) 


No. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 
17 


Date Screen Diameter- In . Voltage - Volts 

1 / 14/72 11 +400 

+- 80 
+600 

1 / 17/72 +- 80 

-350 

-140 

-350 

-300 

1 / 18/72 -240 

-500 
-720 
-760 
-120 
-400 
- 40 
-200 
-700 


18 


-680 
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No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 


TABLE IV 


Voltage Data from 60 Lb. LO^/LH? Experiments 
(LO 2 is poured into LH 2 ) 


Date 

Screen Dia. - in. 

Voita*e--volts 

1/27 

25 

+1200 

1/28 


+3200 



+ 800 

2/1 


+ 800 


+1600 
+1600 
+ 80 


2/2 + 400 

+’*00 

+1100 

+2800 

2/3 +1620 

+ 200 

2/4 +1300 

- 500 
+1000 


17 


600 
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No. 

1 

2 

3 

4 

5 

6 

7 

8 


TABLE V 

Voltage Data from 60 lb. LO^/RP-l Experiments 
(L <>2 is poured into RP- 1) 

Date Screen Di*» — in. Voltage- -volts 

2/11/72 25 -80 

-400 
-200 
- 80 
-160 
-390 
-380 
-255 


TABLE VI 

Voltage data from 60 lb. LH^/LOo Experiments 
(LH 2 is poured into LO 2 ) 


No. 

Date 

Screen Dia.--in. 

Voltage--volts 

1 

2/15/72 

25 

-1000 

2 



-1000 

3 



-1000 

4 



- 500 


5 


- 500 
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for the JH^LC^ experiments. Table VI. In this last set of experiments 
the LH 2 was poured into the LO^ in contrast to the first series of ex- 
periments with 60 pounds of propellants where the LO^ was poured into 
the LH 2 - 

The 240 Pound Experiments 

In these experiments , the ones with the largest quantities of 
propellants used, the I fl 2 was poured into the L0 2 in all but the last 
two tests. While in those two the L0 2 was poured into the Ui 2 . 

Figure 74 presents the traces of a typical experiment and Figure 
75 the traces of the last experiment in this series, carried out the 
morning of March 2, 1972, which resulted in auto -ignition and the 
predicted explosion. 

Table VII and Table VIII give all the data for the 240 pound ex- 
periments with the voltage range set at 10,000 volts for the first 
day and at + 5,000 volts for the remainder of the experiments. Since 
this series terminated with the explosion, no further runs could be 
made at this time and the LO^/RP-l experiments were postponed. 

Discussion 

The purpose of the L0 2 /IH 2 and LC^/RP-l experiments was to verify 
the Critical Mass Hypothesis which is based upon an ignition source 
inherent in the mixing process or, in other words, one which will 
always be present whenever liquid propellants are mixed. 

Among the many possible ignition sources, electrostatic charge 
generation with attendant measurable voltages is always present. This 
was shown in the laboratory with inert components such as LN 2 and RP-1 
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Fig. 74 Typical Chart Traces for a 240 Pound Test 
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Fig. 75 Chart Traces for the 240 Pound Test Explosion 


Voltage Lower Tank Upper Tank 

Generated LH~ L0~ 
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IAMJB VI I 

Data £ro« 240 lb. LH 7 /LO 7 Experiment 
CUI? is pecreed Into 1X»2^ 


Screen Dia. - in. 


Voltage- -volts 


+ 400 


+ 250 


+ 600 


- 200 


- 300 


- 300 
>-1700 


> -1300 


VA&jul VII I 

Data trout 240 16. LCH/LH^ Experiment 
CL02 is poeuretf Into LHj) 


Sc reen Dia. - In. 


Voltage-- volts 


+1000 


-3300 
+ 300 


+ 300 


+2500 (Spike) 
-5000 (Spike) 


+2250 


>+5000 (Spike) 
>-5000 (Spike) 
-3300 (Spike) 


- 500 


- 700 


-2100 
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and others. The question, however, remained whether this would be 

true with the actual propellants used in these tests. All theory 

indicated that there should be the same basic behavior. The voltage 

and charge generated were predicted from the physical constants of 

the constituents and the fonnulation of a physical model, the "Fluid 

Plug Model", which lends itself to mathematical treatment* *30,44,83,104,116. 

Basically this is the concept that one fluid falls into the other 
like a plug and the heat transfer between them will vaporize the other 
proportional to the contact surface area. The vapor formed makes the 
fluid plug bob or oscillate giving it a motion similar to that observed 
in the laboratory with inert systems. 

Then electrical analogy was used to transform the mechanical system 
into an electrical system. This allowed the determination of elec- 
trical properties such as voltages, charges, etc. The mathematical 
treatment involves basic heat transfer equations and vibration rela- 
tionships coupled with the mechanical to electrical analogy. 

The motion of one fluid layer across the other produced electric 
double layers of rather high voltages and since the constituents were 
good dielectrics these voltages built up until discharge occurred across 
one of the many bubbles formed in the mixing and boiling. 

High speed photography of glass contained inert systems showed the 
size of the great majority of bubbles to be about 1/4 inch in diameter; 
from this the necessary charge to produce discharge and consequent 
ignition was determined. 

The theoretical work and the inert laboratory experiments were 
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in excellent agreement and were reported in part in the "Prediction 
of Explosive Yield and Other Characteristics of Liquid Rocket Pro- 
pellants"^* 1 ® 4 by Farber , et. al. and in more detail in "Critical 
Mass (Hypothesis and Verification) of Liquid Rocket Propellants" 11 ® 
by E. A. Farber. 

The tests described in this report were planned, coordinated, and 
executed at the Kennedy Space Center to determine how closely the 
previous theoretical and laboratory work predicted the actual behavior 
of the real potentially explosive rocket propellants. 

Ihe most important relationship predicted both by theory and 
laboratory experiment was that there exists a limit to the amounts of 
liquid propellants which might be mixed without an explosion. Thus 
there is an upper limit to the size of explosion which can be produced 
regardless of the actual quantities of propellants present. The pro- 
pellants not taking part in the explosion may contribute to the sub- 
sequent fire. 

To show this relationship between charge generation and mass, 
tests were planned with quantities of 6, 60, and 240 pounds of propel- 
lants and if necessary with 5,000 and 25,000 pounds. 

The Six Pound Experiments 

These tests were the first ones in the series. Because of the 
small size they were easiest to handle, and the explosion probability 
was less. None or possibly only one explosion was expected for the 
number of experiments planned. 

The voltage instrumentation range was set at +. 100 volts since 
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enough sensitivity was desired to measure any voltage and charge 
generated even though it was expected that if all parameters were 
properly controlled the generated charge and voltage would by far 
exceed this value. 

The first two experiments showed that in both cases the pen of 
•the brush recorder rather violently reached the end of its travel 
and stayed there. 

Hie range was then changed by one order of magnitude and the 
following day the next three experiments again each produced enough 
voltage to "peg the pen", indicating that more than + 1,000 volts 
were generated in each case. 

The range was again increased by one order of magnitude and this 
time the actual readings were obtained with the maximum reaching 
6,400 volts. From the size of the screens and their spacings the 
charges which were generated could be calculated . 

After these experiments with LO^/U^ which produced rather high 
voltages in every test (thus showing that voltages and charges are 
produced in every case of mixing of these fluids) the experiments 
were repeated with L0-/RP-1. 

Again voltages and charges were generated in each case . Hie 
actual values were smaller than in the LO^/IH? case since the energy 
available for the mixing process was smaller . 

Tables II and III give the data for the six pound series of 
experiments . 

It was felt that the results from these experiments verified the 
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predictions and thus no more tests were required and the program 
moved on to the 60 pound experiments. 

The 60 Pound Experiments 

With the experience gained from the first series of experiments 
the 60 pound test equipment was set up and the test procedure repeated 
with L02 being poured into New values were worked «ut for evap- 

oration rates and the corresponding liquid leve 1 ’rop rates. New 
filling times and proper probe settings were c lined. 

The explosive likelihood in these experiments also was higher 
with the expectation of possibility of one or two of the experiments 
auto- igniting* 

The probes were set at levels which allowed proper timing to 
carry out the necessary change-overs in between the filling of Ilk 
and L0 2 tanks. The voltage range was left at + 20,000 volts to save 
time and not require recalibration, and the change to a smaller range 
postponed. 

A few preliminary experiments were necessary until the timing 
had been worked out satisfactorily and then the experiments were 
carried out just as were the smaller ones with the voltages recorded 
and the charges calculated for the larger screens. Since the screen 
areas in these experiments were 5.2 times as large, equal voltage 
generation indicates a charge 5.2 times as large. As shown by equa- 
tions (90), (91), (92), and (93). 


Given 


Q-CV, 


(90) 



Where, 


197 


Q * charge -coulombs 
C ° capacitance-farads 
V ■ voltage-volts 

then for the same voltage generated V, we have by division: 



C 60 V 



%0 


C 60 0 


(91) 


Where, 

Q b0 ■ Charge developed in 60 pound Experiment. 
Qg = Charge developed in 6 pound Experiment . 
C 60 * Screen capacitance in 60 pound Experiment 
Cg ■ Screen capacitance in 6 pound Experiment . 


Also c « 0.224 (n-1) (92) 

Where, 

A * screen area - in^ 

D s screen separation - in (constant) 
n * number of screens (constant) 

K * dielectric constant (3.3) 

C ■ capacitance - yy farads 


then substitute into (91) and simplify; (due to experimental similarity 
the K, d, and n values are equal and cancel) 
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5.2 


(93) 


therefore * 5.2 (for the sane generated voltage). 

A series of these experiments was carried out. Voltage and 
charges were again generated in each case with the larger screen having 
more of an averaging effect. 

When it was considered that enough of these tests led been con- 
ducted the next set of experiments involving LO^ being dialed into 
RP-1 were started. The probes were changed to respond to the proper 
propellant temperatures and the voltage instrunentation range was set 
at + 1,000 volts. The voltages generated and the corresponding charges 
were again smaller than those produced with IX^ and U^. This was in 
perfect agreement with the predictions and after it had been established 
that the voltage and charge generation was present in each of the cases 
the last set of experiments was initiated. 

These experiments involved again IX^ and but this time the 
Ulg was poured into the L£> 2 . The voltage range was set at ♦ 10,000 
volts. No significant difference was found in the behavior of the first 
and last series except that the voltages and charges generated with 
w, dumped into the Lf^ might have generated slightly higher voltages 
on the average. This can be explained by the greater energy supplied 
to the mixing process by dumping the heavier LC^ into the lighter IJ^. 

No explosion occurred in this series of experiments, although 
some of the experiments, based upon the theory, came close to auto- 
ignition. 

The actual results are reported for these experiments in 
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Tables IV, V, and VI . 

The 240 Pound Experiments 

Since the screen area of the 240 pound test is 13.2 tines the 
six pound screen area, equal voltage generation indicates a charge 
13.2 times as large. 

This nay be shorn as in the 60 Pound Experiment by substituting 
into equation (90) , and obtaining by division: 


Where, 


Q 240 . W 

5T «7“ 


or 


Q 240 


S4O 


% 


(94) 


g 240 - Charge developed in the 240 Pound fctperaent 
% » Charge developed in the 6 Pound bperinent 

* Screen Capacitance in the 240 Pound Experiment 

* Screen Capacitance in the 6 Pound Experiment 


Combining Equations (92) and (94) yields 



*240 

*T 



13.2 


(95) 


Therefore 


Q 240 * 13.2 Qg (for the same generated voltage) 


In these experiments it was expected that an explosion would 
occur; from experimental history it was predicted that one cut of ten 
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experiments with LO^/LH^ should auto- ignite. 

Because of the chill down times and then the filling times it 
was decided to start with the IH 7 in the upper tank and the LO., 
in the lower tank. 

Again the experiments were carried out in the same manner as 
the earlier ones. The timing improved with the number of experiments 
carried out and the last three experiments of the first day of this 
series turned out to be the best ones. It is believed that the first 
of these came very close to auto -ignition. 

The last experiment of that day was cairied out be reversing the 
propellants, having the LO^ in the top tank and the in the lower 

tank. 

The first experiment of the next day, carried out in exactly the 
same manner as the others, actually produced auto-ignition and the 
consequent explosion. These data are shown in Tables VII and VIII. 

The 240 Pound Auto-Ignition 

The tenth experiment of the 240 pound series was carried out the 
same way as all the others. Since it was the first experiment of the 
day, the tanks were first chilled down with to conserve L0 2 and 
Then the lower tank was filled with Li^ to the desired level and the 
tanker truck was disconnected and driven away. 

The upper tank was filled with to the required level. Then 
a short waiting period reduced both tank levels to the exact levels 
desired. At this time the L0 2 from the upper tank was dumped into the 
lower one containing the W 2 * 
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The experiment looked just like all the others with a vapor 
cloud appearing from the mixing process. It grew as the. others a::d 
up to this point the test proceeded just the same as before. 

Then Dr. Farber observed a small red glow looking like reflection 
appearing near the opening of the lower tank and changing to yell *w 
and then white as it grew. It then rapidly developed into the deto- 
nation and explosion with shock wave, fireball, etc. of & typical liquid 
rocket propellant explosion. 

Figure 76 shows six frames taken during this test by the high 
speed movie camera which was running at 1600 frames/sec. Frames 2 
through 6 were taken sequentially prior to and during the explosion. 

The development of the vapor cloud and the explosion can be seen clearly. 

The prediction of an auto- ignition with these quantities or pro- 
pellants was satisfied and the prediction that one of ten would produce 
an explosion (a rough estimate) was exactly met, probably more by 
coincidence than by design. This again indicates that much can be pre- 
dicted from theory. 

Figure 77 shows a plot of the generated charge ratio (to bring all 
measurements to the same basis) versus the quantity of propellants 
involved. 

In this plot the lower bound of all the charges generated can be 
extrapolated to produce auto- ignition with a mass of about 2300 pounds, 
a value predicted about four years earlier^*®^*^. When that value 
of charge producing auto-ignition is reached by a smaller mass it will 
explode. 
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MIXED WEIGHT OF LOj AMD LH { (5/ 1 RATIO) IN LB. 


Fig. 77 Charge K.ilio as a Function of Propellent Weight (Ll^/UI., Mixtures) 




204 


Again it should be emphasized that the theory and its verifi- 
cation by these experiments predicts that auto- ignition becomes a 
certainty whan the masses become large enough reaching the Critical 
Mass. Below the critical mass quantities the probability of an ex- 
plosion due to auto-ignition should decrease with decreasing mass, or 
in other words the smaller the mass the less likely is an auto-ignition. 
Since the probability of auto- ignition is not zero, however, for 
masses smaller than the critical mas% explosions with these quantities 
can and do actually occur. 

In these experiments there were no auto-ignitions witu the six 
or 60 pound masses but one out of ten with the 240 pound masses. 

This explosion in the 240 lb. series verified the hypothesis of 
auto- ignition from self generated charges and voltages in the mixing 
of liquid propellants. The voltage spike measured was more than ample 
to produce the electric discharge which triggered the explosion, and 
it exceeded the predicted value needed for auto- ignition. Since the 
explosion was considered the successful conclusion to this test series , 
these experiments were terminated . 

240 lb. LO^/RP Auto-Ignition Experiments 

Some time had to be taken to rebuild the experimental set-up 
and also ground zero was moved further from the control center where 
it was to be set up originally but was changed because of lack of 
cryogenic lines. 

Also because some of the shrapnel action during the 240 lb. 
LO 2 /IH 2 explosion more barricades were erected for protection of 
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equipment and personnel. TV monitors were provided so that the ex- 
perimental procedures and phenomena could be observed closely. 

Since more elaborate testing procedures and also more compli- 
cated data acquisition was used, involving more personnel requiring 
closer coordination, a more formal test procedure was established to 
the satisfaction of all concerned. 

After all this was done, the experimentation was started exactly 
like the one for the 240 lb. LO^Ii^ experiments. 

The results obtained from these experiments were again recorded 
and the important data is given in Table IX. 

charges and voltages produced by the mixing of the W 2 and 
the RP wore again characteristic of these experiments producing spikes, 
but none were large enough to actually produce ignition and conse- 
quent explosion. 

Figure 78 is a plot of the normalized charge versus the mass of 
propellants used and again gives the same curves as the one for LC^ 
and but the value predicted for ignition several years ago and 
indicated on the graph was never reached. It is again seen that the 
lower bound of the measured values merges with the predicted value 
needed. As a further check also the average value at each mass 
quantity produces a line which also merges at that value. The charge 
and voltage which is needed to set off the 240 lb. mixture is the 
minimum and also the average value. 

Non- Instrumented Experiments 

The question was raised as to what effect the instrumentation may 
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Fig. 78 Charge Ratio as a Function of Propellant Height 
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Table IX 

Voltage Data from 240 Pound LC^/RP-l Experiments 
(LOj is Poured into RP-1) 


No. 

Date 

Screen Diameter - In. 

Voltage - Volts x 13.2 

'it 

9/14/72 

40 

1,850 

2 



12,400 

$ 



1,850 

4 

9/15/72 


12,500 

S 



9,120 

6 

- y 


2,510 

7 



16,300 

8 

9/18/72 


7,900 

9 

9/25/72 


4,750 

10 



1,580 

11 

9/26/72 


1,580 

12 

9/27/72 


1,580 
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have on the auto- ignition phenomena. From analysis it was deduced 
that tiie screens have if anything mi averaging effect since they spread 
the produced charges over a larger region and do not allow the con- 
centration which otherwise may occur. 

When the equipment was rebuilt and set-up no internal instru- 
mentation was used in connection with the experiments , nothing was 
put into the tanks which may effect the charge and voltage generation 
or any other phenomena. 

Since no instrumentation was used in the tanks, the liquid 
levels in the tanks had to be estimated by timing. This procedure 
was checked out earlier and produced accurate predictions as to 
liquid levels . The chill time, the filling time and the evaporation 
time were closely monitored. 

With what was estimated to be the correct quantities of propel- 
lants in the tank, the second test produced the expected explosion 
indicating that the phenomena without instrumentation if anything 
allows the ignition conditions to occur more easily. This is exactly 
what was predicted and expected. 

Fully Instrumented Auto-Ignition Experiments 

Again the equipment was rebuilt and set-up at ground zero and 
this time the objective was to measure and determine as many parameters 
as possible . In addition to the probes to determine the liquid levels 
in the tanks and the screens which measured the charge and voltage 
generation, three high speed cameras were used to record the phenomena, 
a thermocouple grid was put into the lower tank consisting of 30 couples 
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in four layers, a photomultiplier tube to sense any light from the 
ignition and explosion before it was destroyed by the explosion, and 
two antennas, one looking immediately over the lower tank, the other 
a little farther away. 

The thermocouple grids were to measure the shock front and the 
reaction front velocities as well as mixing, etc. as in the tests 
described earlier which were carried out a* the Air Force Rocket Pro- 
pulsion Laboratory and Edwards Air Force Base. They were made of #36 
gauge teflon coated copper-constantan wires manufactured irto a single 
cable. The thermocouples were carefully welded under a unall 20 power 
microscope with a minimum of mass at the junction. 

As in the previous instrumented experiments, they were then fas- 
tened to thin support wires and passed from the tank through a conduit 
and to a junction box located at ground level. 

A series of blast gauges were set at ground level in three lines 
120 degrees apart radiating out from the lower tank. 

The two antennas were to pick up signals from the discharges as 
a result of the fields generated by the mixing processes. 

All the signals from all this instrumentation were transmitted 
back to several instrumentation vans where much of the information was 
put on tape and strip charts. 

Again the tests were conducted as previously and all the data 
recorded. The mixing was monitored with the thermocouples and the 
voltage generated measured with the screens. The eleventh auto- ignition 
experiment produced ignition and a subsequent explosion. 
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All the data was recorded from this explosion and the information 
which was on tape was later reporduced on strip charts for analysis. 

Again the expected number of experiments necessary to produce an 
explosion with a certain mess was verified and the yield as determined 
by comparing the blast gauge data to published literature war as expected. 

Some time had to be taken to pull ail the information off the 
tapes and to put it on strip cl .arts so that the various traces could be 
stabilized. Unfortunately the paper used by NASA was of the non-permanent 
type and the records darkened as they were used and a second set had to 
be obtained to complete the analysis. 

The analyses of the traces were simpler in this case than on the 
earlier experiments since the traces had more separation between them 
and it was easier co follow all the variations. 

A. Velocity Vector /Analysis 

This method was based on the average velocity components obtained 
from the ratio of the distance, As, between two thermocouples to their 
difference in response time AT as the reaction front passes over them. 

Case I The reaction front (assumed spher.) is oriented perpen- 
dicular to a line between two thermocouples and moves as shown in 
Figure 79A. This would give an average velocity of the reaction front 
directly equal to ^ and would be located as a good approximation on 
the midpoint of the line between the two thermocouples. If the reaction 
front moves in any other direction the ratio becomes a component of the 
front velocity only. 
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Case I I The reaction front (assigned planar) is oriented perpen- 
dicular to a plane containing three thermocouples and moves as shown 
in Figure 79B. This is essentially a two dimensional case with the 
edge of the front appearing as a straight line. 

We can see in Figure 79C the relation between the velocities of 
a point A on the reaction front and a coincident point B roovi , along 
the line connecting the Thermocouples 1 and 2 can be found from the 
relative velocity equation 


V A 


B + V A/B 


(96) 


The direction of v^ is perpendicular to the reaction front and 
v^/b is parallel to the front (Figure 79D) . Similarly, if a pair of 
coincident points C and D be taken with point C in the front and point 
D moving along line 1-3, we have the relation 


V D " 


V C + V D/c 


(97) 


Again the direction is perpendicular to the front and Vjy^ 
is parallel to it. Since the front is assumed to be planar and trans- 
lating Vp = v A then we can see from the figure that the line ter- 
minating the known vectors Vg and is also parallel to the reac- 
tion front and from the vector polygon Figure 79E, the velocity of the 
reaction front is seen to be a vector perpendicular to this line and 
is also terminated by it. 

Ak l-2 

Thus in this case if two average velocity components v p = ^ 

1 1-2 

- As l-3 

and Vp = ^ are known, v^ the average front velocity can be found 

from the velocity polygon. This average velocity is located as a good 
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approximation at the centroid of the triangle formed by the three 
thermocouples . 


Case III The last case is the most general in which reaction 
front (assuned planar) passes obliquely in space through 4 thermo- 
couples not all on the same plane. By an analogous development it can 
be shown that in space the reaction front velocity will be per- 
pendicular to a plane terminus of three component average velocity 
vectors - and will also be terminated by it. (Figure 79). 

As. . As. . A$. ^ . 

Thus if v n = — and v„ = -■ and v n - 


B At 


1-2 


'C At 


1-3 


’D At 


1-4 


are known v. can be found. 

A 

Again v^ is located as a best approximation at the centroid of the 
tetrahedron formed by the four thermocouples. This method could be 
carried out using three different techniques: 

1. Graphically using projections of the space vectors and 
descriptive geometry. 

2. By vector analysis using the equations 

a. Equation of the reaction front plane: 


$ * v B ) = 0 where $ = (v R x v c + i> c x v p + v D x v fi ) (98) 

b. Equation of the reaction front velocity: 


* 




V 


(99) 


3. 


By analytic geometry using the equations 
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a. Equation of the reaction front plane 


Ax + By 

♦ Cz + D = 

0 


where 

v By V. 1 

. 

v Bx v Bz * 

A * 

v Cy v Cz 1 

» B * - 

v Cx v Cz 1 


v Dy v Dz 

1 

v Dx V Dz 1 


v Bx v B ) r 1 


V BX v By v Bz 

C = 

v Cx V CY 1 

, and D = 

v Cx v Cy v Cz 


v Dx v Dy 1 


v Dx v Dy v Dz 


b. Equation of the reaction front velocity 

v. = £ — , vdth direction numbers (A,B,C) (101) 

A b 2 7 c 2 ' 

These methods gave a first approximation of the location of the 
ignition source in a region 6 to 8 inches above the top level of the 
thermocouple grid and near thermocouples 3 and 4. This allowed the 
use of a second method of analysis shown in part B. 

(B) Reaction from Displacement versus Time Analysis 

After the ignition source had been approximately located, several 
plots of front displacement versus time were made from points in this 
region. Theoretically the proper location should give a plot forming 
a smooth curve. Actually scattering was observed in the plots. This 
was attributed to the differences in reaction times of the individual 



thermocouples due to different bead sizes. The best plot showed a 
clearly defined group of thermocouples on a smooth curve (practically 
a straight line). This group would represent those thermocouples 
of the same size and the same reaction time. Method A was iterated 
with the select group of thermocouples as in the following example. 

Example: Using the analytic geometry equations Tart A Method B 

Given Thermocouple group 1,31,28,11 


Number Location in ft. (from #1) Time in sec in sec Velocity Comp. 








ft/sec 



X 

y 

1 



V 

V 

V 







X 

y 

4a 

1 

0 

0 

0 

3.65787 

0 

0 

0 

0 

31 

1.63 

1.53 

.582 

3.658"3 

.00016 

10,187 

9,562 

3,637 

28 

.822 

.083 

.582 

3.65797 

.00010 

8,220 

831 

5,820 

11 

.822 

.083 

.166 

3.65795 

.00007 

11,743 

1,187 

2,371 


Substituting 

in equation (100) 

wc have 






A «* 29.332 

B - .9036 

C “ 30.058 

D * -416.82 


then 


D 


A A^TTc? 


9,921.8 ft/sec 


with direction numbers (29.33, .903, 30.06). 
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In a similar manner other groups of thermocouples gave the following 
data: 


Group 

Velocity 

ft/sec 

Direction Numbers 

1,4, 9, 6 

10,601.1 

29.33, 

.903, 30.089 

9,10,19,16 

8,897.3 

21.23, 

54.597, 37.647 

11,28,10,19 

13,111.7 

53.91, 

-53.54, -61.53 

9,19,6,8 

7,563.9 

88.54, 

281.3, 182.4 


This allowed a refinement of the center location to a point 8” about 
the top level of the thermocouple grid and 2.b’' from 34 thermocouple 
horizontally and 8" from #12 horizontally (about 1.7" from tank wall). 
See Figure 80. 

Part B was iterated again from this new center and gave a smoother 
plot with less scatter (Figure 81) . Applying the least squares method 
to the data to obtain the slope of the plot gave a front velocity of 
10,600 ft/ sec as the best value of the reaction front velocity in the 
tank. 

Separation between the shock front and the reaction front was not 
observed as in previous experiments since due to the design of the 
experiment the whole container was filled with the mixing region. 

Also the fireball was not observed since no over-all fast cameras 
were set up and all of them were monitoring close in data. 

Time Relationship Between the Various Instrumentation Systems 

In the instrumentation set-up each individual part had its 
characteristic tune constant and also the transmission lines between 
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Figure 81 Reaction Front Displacement Versus Clock 
Time, Determined from the Thermocouple 
Response Times Shown 
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components effected the response to the signal input. 

The cameras were operating at different speeds and the timing 
lights were located at different points with respect to the shutter 
so that the physical event recorded was in a different spot on the 
film with regard to the universal timing code. Also since the physical 
event could have occurred while the shutter was closed, a certain 
amount of uncertainty as to absolute time for a specific event exists. 

The thermocouples had a certain amount of delay in response, but 
since they are made very similar the variation between them is not 
too great. Each will respond somewhat late to a signal but the dif- 
ferences in time between measurements are not great. They exist 
however, and therefore the leading edge of the plot was taken as the 
wave front since the data both bunched up and indicated the most 
accurate timing. 

Taking the data and allowing it to be shifted in time within the 
interval of uncertainty then allowed the lining up of the physical 
events as shown in Figure 82. Obviously all instrumentation recorded 
only one and the same physical event, the explosion, and so the event 
must have occurred only at one time and each instrument recorded it 
at the same time although taking the uncorrected records just as they 
were put on tape or film it seems that the same event was recorded 
at different times. 

Figure 82 shows how the different instrumentation systems line 
up. The thermocouples sensed the explosion first, naturally they were 
located closest to the point of ignition. Then the cameras recorded 
the first flash of light as the explosion emerged from the mixing tank. 






* Position of Shutter "Open time" increment 
may vary within total time interval. 

Note - No meaningful data vaa obtained from the 
photomultiplier, or the antennas. 

Fig. 82 Time relationship of instrumentation during 
last 240 lb. test explosion 
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Last in time the pressure sensor at ground zero picked up the explosion 
and recorded it. 

It must be realized that the time intervals between these recor- 
dings are very small but nevertheless existent and must be taken into 
account for exact analyses and interpretation of events. 

All the time constants and delay times and camera dark times , etc . 
evolved through numerous discussions with the KSC personnel responsible 
for the particular instrumentation and circuitry of each measurement. 
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PART III 

SATURN V DESTRUCT SYSTEM ANALYSIS* 

Having developed methods for ana’ zing and predicting, both 
theoretically and experimentally through modeling, the behavior of 
liquid rocket propellant mixtures and their explosive yield as described 
in the previous parts, this section will be concerned with the applica- 
tion of these methods to actual systems. 

Since the more certain prediction of explosive yield is a function 
of the mode of failure, which in many cases is difficult to predict, it 
was suggested by the principal investigator to analyze the Saturn V 
destruct system. Since, if it should give a low explosive yield, it could 
be used to control the mode of failure of the Saturn V in case a 
catastrophic failure is imminent. Rather than letting a failure occur at 
random and not knowing what to expect it is better to have the rocket 
fail in a certain predetermined manner and have the knowledge of what 
to expect under these conditions. 

With this in mind the analysis of the Saturn V destruct system was 
undertaken, and the analysis discovered a number of interesting facts. 

A. Theoretical Analysis 

For the purpose of the theoretical analysis several models were 
selected, and since the results were amazingly similar, the simplest 
ones were believed to be the best, or at least most useful. 

The first step in the analysis required a knowledge of the mode of 
failure and to this end it was necessary to obtain information as to 
the effect of the destruct system upon the Saturn V. This proved to be 
much more difficult than anticipated. Since this information was to be 

* 

Footnote at end of Part III, p. 34S. 
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furnished by NASA a number of knowledgeable individuals were consulted. 
Much disagreement was found among the groups in NASA. 

120 

Finally we received a letter from Marshall Space Flight Center 
which spelled out the effect of the destruct system upon the Saturn V 
vehicle, and since this was the only infoimation which we received in 
writing it was used as the basis of this analysis. 

If better infoimation should become available at a later date the 
analysis as presented here could be repeated using the identical relation- 
ships buv different numbers. 

Figure 83 shows the effect of the destruct system upon the Saturn V 
vehicle as given to us by NAS/.. 

The destruct system was to cut slots into the various tanks and let 
the fuel disperse out one side and the oxidizer out the other. This 
seems to be fulfilled except for the third stage and this problem will 
be treated later and in greater detail since it seemed to hold a poten- 
tially dangerous situation. 

Thus the analysis then, except for the above mentioned problem, 
consists i. the determination of the following: 

1. Outflow for each stage with respect to 
time for the slot encountered. 

2. Determination of the quantities of liquids 
on the ground at any time. 

3. Determination of the spill ellipses, including 
hydraulic jump, etc. which determines the puddle 
depths at different points or regions. 



Payload 


Instrument 

Unit 



S-IV B Stage 


»7" dla. hole in 
bottom of LOX tank 




'V-/ 


3‘ x 22' slot in LH^ tank 


S-II Stage 

3' x 18' slot 
in LOX tank 


2' x 41' slot 
in LOX tank ' 


S-I Stage 


=s=U)X 


II fi 


f. LOX 


^ " 


■2' x 31' slot in LH 9 tank 


'2' x 20' slot in RP tank 


Figure 83 Schematic Diagram of Saturn V, with Effects of 
Destruct Initiation Indicated— NASA Information. 
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4. Determination of the puddle size and 
the spreading velocity at any time. 

5. Determination of the liquid stream 
impact point at any time. 

6. The r^otting of the spill areas at 
sequential times. 

7. Determination of the amounts of liquid 
vaporized both due to ground action (small) 
and due to the mixing of the various 
components. 

8. Determination of the mixing characteristics 
in the overlapping puddles of both fuels 
and oxidizers. 

9. Determination of the mixed areas and 
volumes for LO 2 /RP-I, and LO 2 /LH 2 . 

10. Determination of ^as volumes - the gas cloud 
forming a hemisphere proportional to the 
areas (arsuming no wind). „ 

To determine the above the following procedures were used. 

Having accepted the mode of failure or the effect of the destruct 
system as given to us by NASA, computer programs were set up to determine 
the outflow of the liquids from the slots in the tanks. 

The next few pages show the theoretical treatment of t'r e tank 


outflow. 
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Liquid release rate from the tank 

Efflux velocity through the slot . 



"5 - mean head in jet 
- initj a! level in tank 

figure 84 Liquid Outflow (Slot) 

Figure 84 gives a schematic sketch used in the following 
ahalysis. 

At the instant w hen the liquid is at h^ , the velocity V, 
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at any point z in the jet is given by 


V 2 - V + | (Pi - p) + 2gh 


( 102 ) 


where. 


h - 

6 * 

P, * 


hi - z 

density of fluid 

pressure in tank at liquid level hi 


From continuity 


A 1 V 1 = A V 


(103) 


where. 


V g = average velocity in slot 
A = area of slot 


The average velocity in the slot can be found by integrating 
equation (102) for V between h = E - b/2 and h = E + b/2 

,E*b/2 


/ 


Vdh 


(104) 


E-b/2 


Substituting for V from (102) into (104) gives: 


„ , 2 ,, . i „ ,, v i .h-'.i, ii , v i , p r p , k . ii. . fl051 

Vs - j(2gt>) c v f + (- 2g£ F i> 1 l l05) 
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where c y is the velocity coefficient and is given by equation 
(103) . 

An iterative scheme, eg. Newton- Rhapson , can be used to solve for 
Vs, knoing pj. 

Equation (105) gives the average instantaneous velocity in the slot 
when the liquid level is exactly at h^. The equation is valid for 
K * b/2. 

When E < b/2, the slot becomes a weir. 

Efflux Velocity Through a Weir 

The velocity at any point z at the instant when the liquid 
level is at h^ , is given by 

V 2 = vj + 2gft (106) 

where h = h^ - z 



Figure 85 


Liquid Outflow (Weir) 
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The average velocity is given by 


f 



(107) 


Substituting for V into equation (107) gives 


1 

v w - fusfl) 7 c v 


[(1 


vf | 

- 


v 2 l 

r * 


(108) 


Equation (108) gives the average instantaneous velocity in the 
weir when the liquid level is exactly at h^. The equation is 
valid for IT < b, or for K ^ b/2. 

The areas of the slot and weir are ab and aFf respectively 
where a » width of slot. Hence, the average flow rates are given by 


JP 

N 

ab V s 

(109) 

% m 

aR V 
w 

(HO) 


Tank Pressure, 

Assuming a polytropic process, then 

Pi ■ Pi<*V<u> k 


cm) 
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where 

k 

k 

*1 


initial pressure at level 
1, isothermal expansion 
Cp/c^ reversible isentropic expansion 
initial volume above level h. 
volume above level h^ 


For uniform area tank with flat top, equation (111) becomes 


J-' h ,k 

Pl U~ 5 


(U2) 


Equation (111) may be written 


p i f l 


where f ^ = ( 


( 113 ) 

( 114 ) 


Time to Lower Liquid Level in Tank , t 


t* / E / h i 

/ at - - / $ah . J idh 

'a r V 


for the slot then 


hi v * p i f h . K . 1,7 


1 3 

+ 5 ♦ I)’ 


- (- 


J. „ f -l 3 

1 - i h + 5 - i)^] _1 dh 


— X- / IC-3L 
|(2gb)7 c y h 2gb p / Jgb b 2 2gb ? /° gh b 2 


(115) 
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where « function of K 
■ function of K 


for the weir then 


Sr 



b _ 1 y2 3 

f U &) 1 [(1 + ™-) 7 

IT 


v 2 I 

(-MY 1 dH 
2gR 


(116) 


where Vj « function of TT 


Jet Velocity at Datum and Jet Trajectory 

The velocity at the ground is given by the vector sum of the 
horizontal velocity leaving the slot, equations (105) and (108), and 
the fall velocity of a particle of water falling a distance d ♦ b/2 
for the Slot and d + H/2 for the weir . Calling this distance D, then 

1 

V g - (V 2 + 2gD) 7 (117) 


and the cross-sectional area of the jet at the ground is 

A - - AV/V (118) 

5 V 8 

g 

Letting D be the mean Vertical distance from the datum to the slot 

* 

or weir, the time t£, for a fluid particle to fall distance D, 
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is given by 

t f - (~) ? (119) 

The horizontal distance x, that the fluid particle travels before 
reaching datum is 

x » Vt f (120) 

Procedure 

1. Fquations (105) and (108) can be solved numerically for 

i 

V and V with given values E and IT and the corresponding 
values of f^. 

2. From continuity ™d (1) above, Vj as a function of E 
and H can be determined. 

3. Equations (115) and (116) can be integrated numerically 
to determine t $ and t w the time required to lower the liquid 
level to the specified values of E and R. 

4. From (1) or (2) and (3) above, the flow rate Q , as a 
function of time can be determined. 

5. The curve of the flowrate Q, as a function of time can then 
be displaced by time t£ so that the velocity at the ground, Vg and 
the flow rate, Q, at the ground can be determined as a function of time. 

6. The horizontal displacement of the jet, x, as a function 
of time is readily determined from equation (120) . 
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7. The area of the jet at the ground can be determined as a 
function of time or distance. See figure 86. 

< , 

Time' required for jet stabilization 

After the destruct system is initiated and the slots opened 
in the tanks it will take a short but finite time to establish the 
stabilized flow. The following will give a short analysis of this 


Using the momentum equation 



from the above 


t - - In ( 1 - 



( 121 ) 


( 12*0 


and 
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(1 



V 2 - — - 2gh ( 1 - e -t ) 

c l 

1 -t 7 

V - - [ 2gh( 1 - e x )] z 


(123) 


The above equation is the outflow jet velocity at any time 
t after the opening of the slot. A quick analysis gives the following 
results. 


Table X 


e- 1 

8 

1 - 0.368 

8 

0.632 

e' 2 

- 

1 - 0.135 

- 

0.865 

e’ 3 

- 

1 - 0.05 

• 

0.95 

e' 4 

8 

1 - C.018 

8 

0.982 

e' 5 

. 

1 - 0.007 

B 

0.993 
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Table XI 


0 


0.795 

0.925 

0.975 

0.991 

0.997 


V 

max 

V 

max 

V 

max 

V 

max 

V 

max 


Shape of the Liquid Pool 

Next it was necessary to establish or estimate the shape of 
the liquid splash pool with the liquid jets impinging upon the surface. 

The derivation on the next pages indicates that the pool is 
essentially circular and the eccentricity of the splash ellipse 
is too small to warrant its consideration in the further treatment of 
this problem. 

Below, the method for estimating the shape of the liquid pool is 
outlined and its variation with time given. 

Basic Assumptions 

1. The pool is elliptical with the stream impact center 
coinciding with the focal point nearest the rocket. 

2. The pool has constant depth as long as fluid 
continues to be added. 

3. The fluid stream is circular with uniform velocity 
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as it hits the ground. 

Figure 87 is a schematic sketch of the physical configuration used 
in this analysis. 

Figure 88 gives some of the nomenclature used in the analysis 
of the splash puddle configuration. 

Derivation 


Vrt. ■ at * ** 

t - C-if-)* 

1 

V wrt = ( 2gz ) 7 (124) 


z - distance from the bottom of 
the slot to ground 


Letting V A and Vg be relative velocities with respect to 


the velocity 


of the focal point of the ellipse 



At tine t 

. _V 

V 


(125) 


and 



V, 


SP 


H 


V SP + V H 




(126) 
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origin 


z 



Figure 87 


Schematic Sketch of Slot, Fluid 
Stream and Splash Puddle 
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Figure 88 Nomenclature 


where in equation (126) Vgp is the spreading velocity due to 
the vertical component of the mass flow. is the velocity caused by 
the horozontal momentum. 

To determine the distances 5?^ and , from the definition of 
the ellipse 


and 


1 

Of a (b 2 - a 2 ) 7 



b - 
b + 


1 

(b 2 - a 2 ) 7 * 
(b 2 - a 2 ) 7 


(127) 


(128) 
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Pram (125) and (128) 


V, 


B 


b - (b 2 - a 2 ) 1 
b + (b 2 - a 2 ) 7 


(129) 


multiplying out 


2 _2,T- 


2 _ 2 >?, 


V, fl> ♦ (b - a ) ] - V B (b - (b‘ - a 4 ) 4 ] 


and collecting tetns 


b (? A - V B ) + (V A + V B ) (b 2 - a 2 ) 7 = 0 


(130) 


also from (126) 


<VV ’ ' 2v h 


- 2 ?, 


SP 


(ia) 


Substituting (131) into (130) 


1 

b V H - Vgp (b 2 - a 2 ) 7 - 0 


and 



(132) 


(133) 
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If Q is the total voltrae on the ground and 
d is the puddle depth 


Q « T a b d 


and 


a 


T b d 


Substituting (135) into (133) 


and 


and 


b2 ^5 " *sp 0,2 ■ 3 > 


V 2 

V H 


V 2 

V SP 


VsP ( 1 ' FFi? 
* ‘ 

t fef 

T* 0& - %) 


1 

,T 


(134) 

(135) 


(136) 

(137) 


Determination of Vgp and 

For the first step assume that the puddle is circular 
with Q cubic feet on the ground and Q being the volume flow rate 
into the puddle. Figure 89. 


* r 2 d. 


2 d r d 


dr 

at 


but since dr/dt 


(Q/ f d) 7 


r ■ 
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1 

Q » 2 d ( <V d ) 7 

1 

Vg, » Q / 2 ( Q d ) 7 

.Q 


C9 



(138) 

(139) 


Figure 89 Splash Puddle 


Let 


V| K)r Horizontal component of velocity as 
the stream leaves the rocket 


A o 

*1 

V 

vert 


(vertical component = 0 ) 

Cross-sectional area of the stream leaving 
the rocket 

Cross-sectional area of stream as it hits 
the ground 

Vertical component of velocity jas stream 
hits ground [ 7 yen - (2gz) 7 ] 


Fran Continuity 


and 


V hor " \ V vert 


A, - -2-^2E , ri . ( V I,J 

^vert 


(140) 


(141) 
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Since 


Q “ A V. 


o hor 


a 'd using (141) 


and from (142) 


* 2 Tj d 


'SP 


1 

7 


1 2d ( TA o W 7 vert> 


H < A o V V vart >’ 


Assuming that 


’SPj SP 

7 V 

v hor V H 


V. 


H 


V V 

V SP hor 


SP, 


using (139) and (143) and substituting it o (145) 


V < C *0 * V vert’ 


1 

7 


Solving for a and b 


(142) 


(143) 


(144) 


(14 S) 


(146) 


from (135) 
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from (137) 


since from 


a 


T b d 


'4> Q 2 


i ’sp-* , 


1 

,T 


[ (S-) z C- 

II d 


1 - (-J 1 ) 2 

V SP 


■)] 


1 

? 


(144) , (146) and (139) 


(147) 



V SP 


V 


hor 


^har A o d ^ ( 4I> Qd ^ 
^ ^vert ^ 


with (142) 


A)2 

V SP 


( t\^ ,d 3 

’ 7 vert ^hor A o 

4 * d Vr 

v 

vert A q 


(148) 


(149) 


and substituting (149) into (147) 
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b 


T d 


1 

1 - ( 1 n d Vhor ^ 

V . A 
vert o 


1 

J 


(ISO) 


From the compute" programs ^hor' A o’ d, and fj 

can be determined at any time t, and with this information a and b 
can be calculated giving the configuration of the splash puddle 
ellipse. 

Taking the first stage of the Saturn V as an example, since it 
gives the largest impact angle and therefore has the greatest effect 
upon the results, they are as follows: 


IX> 2 t - 100 sec V^ or « 6.82 

Q « 37,877 

^ - 3.626 

V 

vert » 74.5 


b - 158.8 ft 
a - 152.2 ft 


t ■ 160 sec 



4.45 

38,700 

1.726 

74.5 


b - 161.5 ft 
a - 152.7 ft 
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RP-1 t ■ 100 sec 



\ * 
V 

vert 


6.30 

18,958 

3.142 

51.2 


b - 113.8 ft 
a ■ 106.2 ft 


t * 160 sec 


b - 

117 

ft 

a * 

107 

ft 


^hor 

Q 

A o * 

7 

vert 


4.16 

19,642 

1.542 

51.2 


From the above it can be seen that the splash puddles are 
essentially circular and using circles makes the analysis much simpler. 

Having shown that the flow pattern on the ground is essentially that 
of circular splash puddles, another phenomenon of the splash puddles , 
namely the hydraulic jump was investigated. This last phenomenon 
produces a variation in thickness of the puddle and this may become 
important in estimating the quantities mixed by having the puddles flow 
together aiid overlap. 

The following pages treat this phenomenon and characteristics. 
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Radial Flow Patterns 

The purpose of this part of the study was to provide a method 
for predicting the depth of flow within a radial flow pattern containing 
a circular hydraulic jump. The results can be used to aid in predicting 
mixing rates of the puddles formed by spilling propellant components 
from a rocket. 

Initially a purely analytical approach was tried. 

Using the notation as shown in Figure 90, expressions were 
developed by applying the equations of conservation of mass, momentum, 
and energy. 

I. Conservation of Mas s (Continuity), (steady flow and 
inconpressibility are assumed) 


a. within impinging jet 


Q 





(151) 


b. at toe of jump 


Q 



*6 V o 


(152) 


c. at heel of junp 

Q - 2lr 1 y 1 V 1 (153) 

d. downstream of junp 


Q - 2 T” r y V 


(154) 



FIGURE 90 NOTATION FOR CIRCULAR 
HYDRAULIC JUMP 
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Equations (151), (152), (153), and (154) may be equated 

assuming no splash before the jump, and no losses due to evaporation 
or ground seepage. 

II. Momentum Across the Jump 

It is assumed that the slope of the jump is constant. 
The forces on an elemental control volume within the jump are as 
shown in Figure 91, and can be expressed as follows: 

f o ■ Po A o ■ \ I'd r o * 

£ 1 ’ Pi *1 * \ r : * 

f s " P s A s * 

The sum of the forces in the horizontal direction, assuming that 

ctm d<|> d<J> . 

SIN 1 if 15 

I 

£ . y 2 0 V y\ rj . 1 (r r r 0 )(y^ * ♦ y*) 

The change in momentum in the horizontal direction is 

a m - ^ qar - I (Vvrd* )(V.- v ) 
g g v o'o o ■ ,v 1 o' 

The required relationship is (making use of equations (15 2) and (15 3) 





2S1 


Vo ' r i + i (r r V^o + Vo + $ 



/Vo 

tvi 


i 



(15S) 


Equation (155) can be written as 



R * 

IK 


1 

+ T 


R + 6F* + 2 6F* 

t — + seamy “ 0 


where 




and F 0 



(156) 


III. Energy Balance Between Jet and Flow Field Up to Toe of 
Jimp 

The following analysis is made assuming that there are 
no friction losses and jet impact losses. 

The amount of energy passing station 1 (Figure 92 during time 
interval AT, is equal to the amount of energy contained within the 
region "S”, where AT is the time interval required to form the 
region "5" at a given flowrate. 





Figure 92 Energy Balance 
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The energy passing station 1 during time interval AT is given by 
m( AT)V? 

E i- -v 

The energy contained within the region "S" is 




i»l 


2g 


dm. * 


ri I* 

J n 


Assuming that the depth of flow is constant from r ^ to r Q , then 
dmf - in r .py 0 dr 


From continuity. 


v B ■ dr 
V i 2* Pgy 0 r 


Substituting Equations (159) and (160) into Equation (1S8) with 
r.< r i < r 0 , gives 


h 



m^ dr 
4wpgy 0 r 


• 2 
m 


4wpgy r 


In 



Evaluating the expression within the limits shown, the element c> 
with radius r^ and depth y Q has been onmited, and the energy c. 
element must be added to Equation (161). This energy is esse;. - is i; 
potential energy, and is given by 

E • m-? y 0 p(H - y 0 ) 

where (H - y Q ) is the height of fall. Since y Q <« H , H - y o * 


(157) 

058) 

(159) 

(160) 


(162) 

H 
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and Equation (162) can be written 


E ® irrT y 0 pH 


Adding Equations (161) and (162) gives 


e s =i p* e ? 4^ i "| * n } v* 

Evaluating AT in Equation (157) 


(163) 


AT 


2 

trr y p 
o 7 o 


(164) 


Substituting Equation (164) into Equation (157) gives 

n l V v i 


w 

Equating Equations (165) and (163) by conservation of energy gives 


(165) 




i - In + mrty pH 

2g 4vpgy o r^ j 7 o 

V 1 m 2 2 A 2. 7 

Since H • and m ■ tt , Equation (166) c<ai be written 

1/2 


(166) 




r i 


4 


2 5-1 


«- r 


J 


(167) 


The description of the flow given by Equations (151) , (152) , (153) , 


(156), and (167) is incomplete, because the location cl the jump, r Q , 
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i 3 unknown, as is the length of the jump, r^- r Q . Since analytical 
methods have been exhausted in the derivation of Equations (151) , (152) , 
(153) , (156) , and (167) , an attempt was made to derive empirical ex- 
pressions containing the unknowns, and r Q . 

Length of Jump 

Experimental evidence indicates that for rectangular jumps, the 
length of the jump varies between 4.3 and 5.2 times the thickness at 
the heel of the jump. An empirical expression can be given by 

r.- r * ny, ; 3.5 < n < 4.5 

1 o # i 

Letting n - 4, and rearranging terms gives 

h - -li 

r o T l_ r o 


•] 


(168) 


Location of Jump 

Local experiments were conducted to establish a relationship be- 
tween the impinging jet's flow characteristics and the radius, r Q , 
at which the junp begins. It was assumed that the relationship is of 
the form 

^ - c!\? (169) 

r j 1 


F i 


UgTj) 


17 ? 


where 
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The results give values for c and b of 0.245 and -0.672 respec- 
tively. Substituting these values into Equation (169) , gives the 
following: 

^ = 0.245 fT 0 * 672 (170) 

Experiments Leading to Equation (170) 


Procedure 


These experiments were conducted on the roof of the Mechanical 
Engineering Building . Hie apparatus was set up as shown in Figure 93 
and 94. The flow rate was varied giving circular hydraulic junps 
ranging in diameter from 2.5 to 11.75 inches. Experiment No. 1 was 
performed without a nozzle, and Experiment No. 2 was performed with a 
3/16 inch I.D. nozzle so as to produce higher jet Froude numbers. 

For each flowrate, the Diameter, D Q , was measured, and the flow- 
rate was determined by measuring volumes and elapsed times. 

The Froude number of the jet, F. , is 


given by 


F. 

J 



TTL 


where 


v i 


Velocity of the jet 


and Dj 


Diameter of the jet. 


Results 

A plot of r./t Q vs Fj is shown in Figure 95. A curve of the 
form r^r^* was fitted to the data, where c » 0.245 and b * -0.672. 
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Equations (156) , (167) , (168) , and (170) are presented graphically 
in Figures 96 thru 99 respectively. 

A plot of the simultaneous solution of Equations (151) , (152) , 
(153), (156), (167), (168), and (170) is shown in Figure 100. After 
making several sample calculations using Figure 100, it was found that 
predicted depths at the heel of the jump were approximately twice those 
observed. Reviewing the assumptions made in the derivation of Equation 
(167) it was felt that the error involved by neglecting the losses 
incurred in the abrupt transition from vertical to horizontal flow, by 
assuming the depth to be constant, and by neglecting frictional losses, 
accounts for the sizeable departure from observed results expt rieneed 
when using Figure 100. Unable to account for losses neglected in 
Equation (167) analytically. Equation (171), using continuity only. 


where a 



r j 


1 

2a 



(171) 


was used to replace Equation (167) . A plot of Equation (171) , for 
various values of a , is shown in Figure 101. A value for a of 
0.14 to 0.10 gives results which best agree with observed results. A 
simultaneous solution of Equations (151) , (152) , (153) , (156) , (168) , 
(170) , and (171) , is given by 


RF b + SF? b - TF^ - UFj b - VF? +b + WF? (1+b ^ + X = 0 


(172) 
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where 

R * 2QA 2 a ; S - 32A 2 ac - 2Ac - 3c/4Aa 
T * 8A 2 c 2 + 4c 2 /a ; U = 4Ac 3 /a ; V = 6a 2 /A 

W = 3ac/A 2 - 24a 2 c ; X * 3Aa/c 

A * y^r. ; a = V 0 /Vj ; c = 0.245 ; b = • 0.672 

A plot of Equation (172) for vari . values of a is shown in Figure (102) , 
Computer programs for the solution of Equations (156) and (172) 
are presented in the Appendix. 


Conclusions 

An effort was made to apply Equation (156) to the regions down- 
stream of the jump, with the limitation that the Froude number at the 
heel of the jump be less than one. Computer programs for the solution 
of the equation, both considering the side forces on the elemental 
volume and neglecting them, are presented in the Appendix . T » was 
found that a very rapid decrease in thickness was predicted immediately 
foil -swing the jump, and b, applying continuity to the predicted di- 
mensions an increase in velocity resulted for an increase in radius. 
This not only contradicts intuition, but also observed results. There- 
fore, Equation (156) is only applicable to the situation in which a 
jump occurs. 

From observation, it is felt that the depth downstream of the jump 
can be considered approximately constant for all practical purposes, 
its value taken as the depth at the heel of the jump. 
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Tb above analysis showed that the splash puddles can be considered 
circular for all practical purposes since the eccentricity is very small. 
This allows for mich less laborious treatment of the problem and a such 
shorter method for determining the nixing function. 

Fron the slot size and basic fluid dynamics the out-flow velocities 
and flow rates as a function of tine can be determined giving the amount 
of liquid on the ground at any tine t . This was done by computer 
program and the results are shown in Figure 103, normalized for the 
Saturn V. 

The impact point also can be calculated by the computer program 
giving the time variation of the iiqpact point on the ground surface. 

The results of these calculations are shown in Figures 104 through 108. 

The depth at any one time t of the puddles was assumed to be 
constant as supported by the theoretical investigation, and the variation 
with time was calculated using the momentum spreading relationships 
during the early times and the hydraulic head spreading at the later 
times. Both these curves were plotted and where they crossed the change 
from one result to the other was made. The complete curves representing 
the depth of the puddles are also given in Figures 104 through 108. 

From the quantity of propellant on the ground and the depth at 
any one time the radius of the puddle can be calculated and the vari- 
ations with tine are again presented in the same Figures 104 through 108. 

Having the impact point for each fuel or oxidizer on the ground 
and the amount of each component on the ground at any tine t , and 
the depth, the resulting radii can be plotted to give the splash patterns 
which allow the determination of the mixing regions and the quantities 
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Figure 103 Volume on Ground from Each Tank Based On NASA Slots 
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Figure 105 Theretical Calculations of Inpact Distance, Radius, and Depth 

of Puddle Versus Prototype Tine for Stage I (RP Tank) . 
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Figure 106 Theoretical Calculations of Impact Distance, Radius, and Depth 
of Puddle Versus Time for Stage II (IXD 2 ) 
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Figure 107 Theoretical Calculations o£ Impact Distance, Radius, and Depth 
of Puddle Versus Time for Stage II (LH 2 Tank) . 



Figure 108 Theoretical Calculations of Impact Distance, Radius 
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which axe actually mixed at any time t . 

Figure 109 gives the splash pattern and the mixing regions around 
the Saturn V rocket due to the activation of the destruct system and 
shows the progressive mixing as time goes on. The actual conditions 
are plotted for 5, 10, 15, 20, 30, 40, and 50 seconds after the ini- 
tiation of the destruct system . 

A scale change was necessitated at this point to be able to pre- 
sent the continuing splash patterns on the same size of paper. So 
Figure 110 presents the continuing splash patterns for 50, 60, 80, and 
120 seconds. The first one was selected to give the overlap with the 
previous Figure and a comparison can be made. 

Another scale change was made for the last part of the splash in- 
vestigation to again shrink the pattern so as to fit it onto the page. 
Figure 111 presents the patterns for times of 120, 180, and 240 seconds 
after the initiation of the destruct system. 

From the splash patterns the average depth at any time x was 
obtained by dividing the total splash area into the known volume on the 
ground. The areas of the mixed constituents LO^IJ^ , LO^RP-l , and 
were then multiplied by the average depth at that time to 
obtain the mixed volumes. The volume of the three way mix was most 
frequently assumed to be one half LC^/U^ and one half LX^/RP-l. Other 
combinations were also investigated such as having all the LH^ react 
before any of the RP-1 was allowed to, or the evaporation was such as 
to produce a certain detinite mixing depth. These simplifying assump- 
tions were made to get sane approximate answers early in the investi- 
gation and after much more detailed work with computers, proved to be 



Splash Patterns for time intervals of 5, 10, 15, 20, 30, 40, and 50 
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Figure 109 Splash Pattern for Saturn V Destruct System 
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Splash Patterns for 50, 60, 80, 120 sec. 



: igure 110 Splash Pattern for Saturn V Destnict System 
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Splash Patterns for 120, 180, and 240 sec. 



Figure 111 Splash Pattern for Saturn V Destruct System 
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rather good and in basic agreement with the later much more elaborate 
methods of analysis. 

For the more elaborate work the evaporation losses were theoreti- 
cally calculated using the methods of Part I. This allowed the deter- 
mination of the mixture on the ground at any time T , considering 
the heat exchange among the propellant components and the propellant 
components and the ground. 

The evaporated material was assumed to form a hemisphere , (no wind) 
above the puddles. The clouds came together and mixed before the 
liquids could on the ground and so they show the first potential for 
producing an explosion and a yield, but the liquids started mixing 
soon after and they had much more energy which could be released in an 
explosion. Using the methods as discussed in detail in Part I, for 
the analysis of the Saturn V destruct system the resulting explosi\ 
yield curve can be plotted. Figure 112. 

It is well to mention again this curve of Figure 112 is the ex- 
plosive yield which can be expected at any time x if ignition should 
occur at that time. From the previous discussions in Part I and Part II 
of this report and the many papers and other reports discussing this 
subject, it has been shown that this curve would actually be terminated 
rather early in its development either by external or internal ignition 
sources so that the actual yield values are held to rather small numbers. 

The third and second stage, S-IVB/S-II Interstage failure is 
discussed in detail later in this section. 
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a. 1/373 Scale Sand Model Analysis 

b. 1/50 Scale Water Model Analysis 

c. 1/28 Scale Water Model Analysis 


Basically the procedure was to build models and then use them to 
determine the mixing function and from it the yield function. This 
was done by having different colored materials (sand or water) flow 
from the model tanks and mix. The analysis of these experiments re- 
sulted in these functions. 

Rather than going through the complete process for each method and 
thus doing repetitive work, it was decided to obtain the mixing curves 
foi each of the models and methods without consideration of evaporation 
and then compare the results at that stage. With good agreement of 
these three experimental procedures, the best resulting curve was pro- 
cessed to finally obtain the modeled yield function. 

From the modeling as will be described for each syst v the percen- 
tage mixed as a function of model time is obtained schematically as 
shown in Figure 113. 
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Figure 113 Mixing Ratio Versus Model Time 


Two curves must be found, one for LO^/RP and another for LO^/UL 
since they result in different yield values for respective unit volumes. 
They will have to be added to obtain the total possible yield at any 
model time. 

The mixing areas under analysis will form basically three regions, 
one where LC^ and RP mix, the second where LC^ and Ii^ mix, and the 
third where ail three components mix. A fourth region where Lll^ and 
RP mix is only of interest for evaporation loss determination or 
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possibly atmosphere interaction. 

In the third region a decision uist be made as to the proportion 
of LO^ mixing with and that mixing with RP. 

The next step in the analysis is to expand the model time to the 
corresponding prototype time. This relationship will be different for 
each of the systems used, but the method of determination is the same. 

Since the outflow for the prototype is calculated by the computer 
program, while that for the models is determined experimentally, the 
relationships of prototype time to model times are known while the 
tanks are emptying. Then with the assunption that the ratios of vol- 
umes mixed to volumes on the ground are equal for both model and pro- 
totype, step by step construction of the V^^xotal versus prototype 
time function is possible during this time. Figure 114. The evapor- 
ation losses due to the mixing and the minor effect due to the ground 
are ignored. 

The step by step construction of the prototype mixing curve beyond 
the emptying time shown in Figure 114 can also be expressed mathematically 
as follows. 


V . 
mix 


Ax 


M i-1 



\ 

/ 



(173) 


The above equation is taken where both the model and the proto- 
type relationships are known so that the next interval can be calculated. 
By selecting desirable intervals successively on an overlapping 
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Computer Program 
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figure 114 


Construction Method for Determining the Prototype 
Mixing Function from the Model Mixing Function 
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basis, the complete prototype mixing function can be constructed from 
the model mixing function. 
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ETC. 


Following the above procedure all the curves for each of the models 
can be constructed, stage by stage or for the total, keeping LO^/RP 
separate from U^/UL. 

These curves can then be compared as shown in Figure 115. 

Now either any one of the curves or an average from the three 
models can be selected and an evaporation factor selected, or the actual 
evaporation calculated at each time from the contact areas of the 
propellant constituents and the heat transfer characteristics. The 
characteristic curves are shown in Figure 115 c,d. 

Since one unit volume of LO^/RP'l has a different amount of energy 
than one unit volume of LO^UL* the ordinates of these curves are 



Figure 115 c 


Figure 115 d 



Figure 115 Explosive Yield Analysis 
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multiplied by their respective energy factors and result in the yield 
which can be produced at any time t from each of the constituent com* 
binations. The sum of these yield components is then the total yield 
which can be expected at any time t . Figure 115 e. 

If the cloud is considered in addition to the liquid propellants 
on the ground, then assumptions have to be made as to the meteorological 
conditions. These are rather difficult to predict in general and so 
such conditions as zero wind velocity, etc. can be taken. 

With the complete procedure laid out, the detailed modeling 
methods are now discussed with the 

1/373 Scale Sand Model 
1/50 Scale Model 
1/28 Scale Model 

Sand Model Analysis 

To determine the mixing function the first modeling method used 
was to take one of the Saturn V models which were on sale at the KSC 
Visitor's Center and to subdivide this plastic model into compartments 
simulating the tanks in the actual Saturn V. The scale of this model 
is 1/373 of the prototype. 

Tnen slots were cut to scale to simulate the effect the destruct 
system would have on the prototype as specified by MSFC. These slots 
were pre-cut after it had veen tried to actually make them by hot wire 
explosions. None of these worked out as satisfactorily as the pre- 
cutting. 
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Hie various tanks were filled with colored sand, red for RP-1, 
yellow for , and blue for LC^. A clam shell was put around the 
model to prevent the sand from flowing from the tanks. Then the sand 
in the tanks was fluidized by air at 3 inches of Hg to make it behave 
like a liquid. 

The model was set upon a platform which rested on a vibration 
table which could be controlled in both frequency and amplitude. 

This complete system is shown in Figure 116. The sand behaved 
like a liquid by-being fluidized in the tank and then when hitting the 
vibrating platform kept on flowing just as a liquid would. 

With this physical system it was only necessary to start at time 
zero by opening the clam shell for let us say 1/2 second. With the 
clam shell closed and the shaking table stopped at that time, the 
splash pattern on the ground and the mixing could be studied at leisure. 
The areas and the depth of the various regions could be traced or 
photographed for later analysis. 

This process was repeated as many times as was desirable and 
thus each record represented a point on the mixing-model time curve. 
Figure ll 7 shows photographic records of the mixing process. 

By the above method each stage can be analyzed individually or 
all of them combined giving the mixing function for the desired con- 
dition . 

By the method discussed in the earlier parts of this section, the 
prototype mixing function can be determined step by step from the 
prototype computer program and the complete model mixing curve. The 
results are shown in Figure 118. 
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Figure 116 Saturn V Destruct System Vibration Analysis Apparatus 
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If this were the only modeling method used here the evaporation 
could be calculated point by point and the yield versus time curve 
obtained, completing the prediction modeling. 

Since in this work, however, other modeling methods were used, 
they will be discussed and the results compared before the mathematical 
reductions are introduced. In this manner, three modeling methods 
were compared and used for the predictions. 


1/50 Scale Model Analysis 

Having carried out modeling with a rather small model and fluidized 
sand, it was thought to be desirable and probably more accurate to to 
to a 1/50 scale model of the Saturn V. 

The possibility of having this model constructed was investigated 
by our group and MSFC as well as private industry, and it was found 
that we could construct this model at a small fraction of the cost 
it would have taken to have it made. 

So for this reason the model was constructed here following the 
Saturn V configuration. Again the slots were cut into the tankage as 
specified by MSFC, and after clamping arrangements and other methods 
of opening the slots quickly were tried, the rip- tape technique developed 
by our group was found most satisfactory. Duct tape was put over the 
slots and then fastened to ropes with pullies so it can be ripped off 
almost instantaneously. Since 5 of these tapes had to be pulled at 
the same time, it was necessary to have a coordinated system available. 

The 1/50 scale model was filled with colored wateT simulating the 
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liquid propellants, again yellow for red for RP, and bit e for LOg* 

High speed cameras were set up which recorded the mixing of the pro* 

* • * 

pellants as they spilled out onto a marked surface. Figure 119. 

' ’• i 



‘ Figure 119 1/50 Scale Model Splash Pattern ' 

Using high-speed films and utilizing the flow out calibration Tuns, 
the areas of the spill and mix could be determined frame by frame and 
the mixing curves determined. These presented the quantity mixed as 
a percentage of the total propellants present in terms of model time. 

The model time was then transformed into the prototype time just 
as was shown earlier in this report and therefore, is not repeated here. 
The end result is shown in Figure 120. 



1/50 Scale Model 



Prototype 
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1/28 Scale Model Analysis 

To have a range of sizes to also get a feeling for the scaling of 
these methods, it was decide to also carry on the experimental ion with 
a model 1/28 the scale of the actual Saturn V. 

Again this model was constructed by our group for a small fraction 
of the costs involved in having it done elsewhere and also at a tre- 
mendous saving in time. Figure 121. 

The model was then used exactly as the 1/50 scale sndel, except 
that the experimental ion had to be carried out at the Solar Energy 
Laboratory about 12 miles from the campus, since the roof upon which 
the smaller experiments were carried out was not large enough to give 
the splash areas room to develop before running off the edges. 

Otherwise the experiments were basically the same. Again high- 
speed cameras recorded the mixing processes for later analysis frame 
by frame. 

The analysis again resulted in a mixing function versus model 
time, and with the method described earlier could be tr.msformed into 
the prototype mixing function, again without evaporation . Figure 122. 

Having finished the analyses for each of the three models, the 
Sand Model (1/373), the 1/50 Scale ..fodel, and the 1/28 Scale Model, the 
results from the modeling can be compared. 

Figure 123 shows that mixing curves predicted for the prototype 
without evaporation are about the same. So the rest of the analysis 
needs to be carried out only once and will give the same result 
regardless of which model is taken as the reference. 










Figure 123 


Experimental Volume Mixed - Time Function 
As Predicted for °.atum V 
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Evaporation Losses 

Since the mixing process of the constituents produces evaporation 
of the cryogenic components, less and less total quantity of propel- 
lants is present and the lc^s occurs in the mixing region. 

If the mixing could be carried out without ignition occurring, 
which has been shown not to be possible, then eventually all the vola- 
tile components would be gone and nothing would be left for any reaction 
to occur. With this, the explosive yield, which starts at zero when 
no mixing has taken place, will return to zero. In between these limits 
in time the explosive yield would reach a maximum when the greatest 
amount of propellants is mixed and available for reaction. 

Since little information is available with regard to evaporation 
rates, etc., our group carried on some experiments with spilling 
an a sand bed which was resting on a scale. 

Weighing the combination of sand and at short time intervals, 
and subtracting out the ice formed from the atmosphere, the evaporation 
rates and the heat transfer characteristics could be estimated. Figure 
124. 

Using the heat transfer characteristics obtained and the informa- 
tion that the contact area is essentially proportional to the profile 
area, the amount evaporated can be determined and the actual amount of 
mixed propellants present at any time t is found. 

Translated into explosive yield. Figure 125, as described earlier 
and expressed in percent of the total, it is shown that the explosive 
yield for the Saturn V through destruct system activation starts to 
rise in value and reaches a maximum of about 9 percent at about 150 
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Figure 125 Experimental Expected Yield Function for the 
Saturn V Destruct System 
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seconds, then slowly falls off. 

It is very doubtful that the actual phenomena can follow this 
curve very long since ignition, either through an external source or 
through internal action such as charge and voltage generation, will 
occur leaving us with a rather low explosive yield. 

This last fact has been supported a number of times with smaller 
rockets and liquid propellant quantities which produced rather small 
explosive yields. 

From the above analysis it is seen that the destruct system is 
essentially doing what it is supposed to do except in the third stage. 

There a potentially dangerous condition seems to exist. Hie 
activation of the destruct system could lead to an explosion in the 
second stage which in turn could produce a considerably larger secon- 
dary explosion in the third stage. The total explosive yield value 
is still very low as far as the total Saturn V rocket is concerned, 
but it occurs near the payload and may through its proximity do con- 
siderable damage to vital systems and possibly personnel. 

A detailed analysis of tiiis case will be made in the following 


pages. 
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S- IVB/S-II Interstage Failure Analysis 

During the Saturn V Destruct System Analysis it was found that a 
serious problem may exist due to the activation of the Saturn V Pro- 
pellant Dispersal System. 

Hie destruct system is to cut slots into all the fuel and oxidizer 
tanks except for the S-IVB LQX tank^'^* where the ring charge will 
cut a 47 inch diameter hole in the bottom of the tank. The objective 
of this destruct method was to disperse the propellants and oxidizers, 
with the propellants pouring out on one side of the rocket configuration 
and the oxidizers on essentially the opposite side. 

The above objective seems to be met in all cases except in the 
case of the S-IVB stage and its interaction with the S-II stage. 

As our analyses, both theoretical and experimental (modeling), 
indicate , the LOX pouring out the 47 inch hole of the S-IVB LQX tank 
will, through its weight and momentum, tear loose the thrust cone 
carrying the rocket engine for this stage . LQX will pour over the 
engine and through port holes in the thrust cone following some of the 
supply pipes which go from the tanks to the engine. 

The S-IVB engine is only a few inches above the very thin Itfg tank 
dome of the S-II stage. The tearing loose of the thrust cone and the 
resulting downward movement of this engine will cut the cone of the 
S-II stage. This sequence of events will allow LQX to pour directly 
into the Ii^ - a phenomenon or mode of failure investigated both under 

22 11Q 

Project PYRO* and the University of Florida program on Liquid 
Rocket Propellant Explosion Qmracteristics*’^*^ , ^*^*^* 5 ®*^ ,7 ®* 
122,123,124,125 . resu iting in an explosion of as few percent of the 
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total possible. 

This primary explosion will occur essentially under the engine 
and thrust cone driving than upward or forward with explosive velocities. 
These velocities may be of the order of several thousand feet per 
second. Such velocities have been measured for large pieces of metal 
in liquid propellant explosions through film analyses of the ADL^ 
and FYRO experiments. 

The forward motion of the thrust cone and engine is very likely 
to break up the common bulkhead between the LOX and in the S-IVB 
stage, thoroughly mixing the LOX with Li^ in a very short time incre- 
ment. This EXPLOSIVE MIXING MODE has not been covered in any of the 
experiments and test programs carried out such as Project FYRO, the 
ADL experiments, etc., but was simulated by the modeling techniques at 
the University of Florida. 

The indication is that the EXPLOSIVE MIXING MODE provides much 
more thorough mixing before ignition than is accomplished by the much 
more gentle mixing methods of the above cited experimental programs. 

But even in those programs it was shown that the greater the energy for 
mixing the greater were the yields which could be obtained. Mr. 

Fletcher’s work® 0 also indicated that more thorough mixing before ig- 
nition produced larger explosive yields. 

Thus our analysis, the only one to the best of our knowledge 
which covers the EXPLOSIVE MIXING MODE, indicates that this secondary 
explosion produced by the primary explosion, of the type covered in 
the ADL and FYRO experiments, can be of much larger magnitude than 
that produced by the non-explosive mixing modes covered in previous works. 
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This effect of the destruct system was not originally expected to 
occur, and even though it does not appreciably effect the over-all 
explosive yield of the Saturn V, still predicted as low, because of 
its location near the payload, it may be more critical in nature. It 
may even trigger further explosions in the LEM, Service Module, etc. 
with the possible destruction of a radioisotope generator. 

This S-IVB/S-II interstage failure or interaction, is man-made, 
forced upon the Saturn V configuration by the design of the destruct 
system. It seems that the removal of the ring destruct charge from 
the S-IVB LQX tank would eliminate the possibility of this particular 
interstage failure. 

Details of the Analysis 

Below will be given some of the details and specific results ob- 
tained from both theoretical and experimental modeling analyses. 

S- I I /S-IVB Design Features : 

Figure 83 shows schematically the Saturn V configuration 
with the location of the destruct charges of the present Saturn V 
destruct system. Figure 126 gives the detail of the important sections 
of the S-IVB and S-II stages involved in the interstage failure analysis, 
the subject under discussion here. 

Effect of the Destruct Charges Upon the Saturn V Tank 

Configuration : 

Figures 127 and 128 give information as to the amount of 
LH 2 in the S-II and S-IVB stages respectively as a function of time 
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Figure 128A Volume-Time Function of LH 2 In the S-IVB Stage. 
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as based upon the failure iiuonnatian supplied by NASA. Hie liquid 
levels are also given. 

LQX Escape Prom the S-IVB Stage : 

The action of the destruct system according to NASA infor- 
mation will cut a 47 inch diameter hole into the bottom of the LQX tank. 
This will allow the LQX to drain out into the thrust cone supporting 
tiie engine. Figure 129 gives information as to the volume and liquid 
level of LQX in the tank as a function of time. 

Forces on the Thrust Cone : 

With LOX pouring out under the influence of pressure and 
gravity through the 47 inch diameter hole in the bottom of the S-IVB 
LQX tank after initiation of the Saturn V destruct system, tremendous 
forces are produced on the thrust cone supporting the engine. The 
LQX jet, the pressure behind it due to evaporization of LQX in the 
LOK tank, not even considering the impact load of the jet, will produce 
a loading on the thrust cone in tension of about 700,000 lbs. The 
impact will considerably increase this value. 

Figure 129 shows maximum amount of LQX in the S-IVB LQX tank during 
the expansion process and while filling the thrust cone. 

Strength of the Thrust Cone and Engine and Puncturing of the 

S-II 1H 2 Dome : 

The thrust cone carrying the S-IVB rocket engine is attached 
to the bottom of the LQX tank by 96 bolts which have to carry the 
entire load in tension, lb determine the strength of this configur- 
ation we took a number of the bolts and tested them in our laboratory 
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both at room aid temperatures. The results of these tests are 
shown in Figure 130. 

The combined strength of these bolts and thus the system under 
uniformly distributed loading (which is very unlikely because of the 
violence of the LOK jet) is less than the force produced upon this 
system even when only static loading is considered. Thus it must be 
concluded with certainty that the thrust cone assembly will tear loose 
from the bottom of the LOK tank and drop, partially guided by the fuel 
and oxidizer lines. 

The sequence of events likely to follow the failure of the bolts 
is shown and described in Figure 131, both schematically and pictor- 
ially based upon model studies. 

Dropping of the Thrust Cone and Engine and Penetration of the 

S-II LHg Tank Dome : 

As discussed above, it can be concluded with mathematical 
certainty that the engine and thrust ccne will drop. After movement 
of only a few inches the rim of the expansion nozzle of the engine 
will contact the doom of the S-II LH*, tank causing tremendous stress 
concentrations , which will cut a hole into this dome letting the 
engine penetrate into the S-II IU 2 tank. LOK Tunning through the holes 
in the thrust cone, down the engine, and spilling over the edge of the 
thrust cone will thus enter the S-II LH 2 tank mixing with Ui 2 below 
the engine and thrust cone assembly. 

Primary Mixing of LQX and LH, 

All the events leading up to this point happened in a very 



Figure 130 A Strength of S-IVB Thrust Cone Anchor Bolts. (At room temperature 
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(a) 



Schematic configuration prior 
to Initiation of destruct system. 


(b) 



Destruct system Initiated; slots 
blown In LHo tanks and circular 
hole blown Tn bottom of LOX tank. 


(c) 



LOX pours Into thrust cone. Some 
LOX escapes through openings (for 
pipes) from thrust cone. 


Figure 131A 

Possible Sequence of Events 
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(d) 



Height and liqpact of LOX In thrust 
cone breaks mounting bolts and entire 
engine assembly falls. Engine 
penetrates upper bulkhead of S-II 
LHg tank* LOX flows Into LH 2 tank to 
start primary mixing. 



(f) 



Engine assembly Is driven explosively 
through lower bulkhead of LOX tank- 
primary fireball expands. 


Figure 131ACont1nued 
Possible Sequence of Events 
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( 9 ) 



Thrust cone-engine assembly continues 
forward, penetrating upper bulkhead 
of LOX tank. Introducing LOX Into 
LH 2 tank of S-IVB. 


(h) 



Engine continues through LH 2 tank 
producing further mixing. 


( 1 ) 



Extreme turbulence and large amounts 
of LOX Introduced Into LH 2 tank 
produces extensive mixing. 


Figure 131A Continued 
Possible Sequence of Events 
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Violent secondary explosion 
occurs In LH 2 tank of S-IVB. 


Secondary explosion fireball 
expands destroying S-IVB Stage. 


Figure 131A Continued 
Possible Sequence of Events 
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Figure 131BPoss1b1e Sequence of Events. 
Letter notations correspond to sketches In Fig. 
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Figure 131B Continued 
Possible Sequence of Events 
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FI gure 131B Conti nund 
Possible Sequence of Events 
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short time. LOK from the S-IVB and from the S-II will continue 
to nix until ignition occurs producing the primary explosion. A 
typical curve for the amounts nixed as a function of tine is presented 
in Figure 132 for both the models which were used in simulating the 
events and for the prototype based upon correlation analysis. 

These curves represent the amount mixed as a function of time if 
ignition does not occur. When ignition occurs, it will terminate the 
curve at the time of ignition. 

Ignition and Predicted Yield for the Primary Explosion : 

Based upon the theoretical investigations of ignition phe~ 
^^^^75,124, 125,126^ ^ can be concluded that ignition of the mixture 
would occur rather early during the mixing phenomenon and thus result 
in a relatively low yield. Experience also supports this. 

The type of failure involved in this case has been covered by the 
University of Florida investigation 79 * 15 * 75 * 124 * 125 ^ ^ by 

Project PYFD 22,1J ^. The S-IV test 22 * 119 with an explosive yield of 
about four percent seems to be applicable in this case. Also, the 
thoroughly instrumented 25,000 lb. explosion experiments 79 * 75 * 124 * 125 
correlate well with the low yield prediction. 

Effect of Primary Explosion Upon Engine and Thrust Cone Assembly : 

The primary explosion, occurring essentially below the thrust 
cone engine assembly, will, based upon vector analysis, be propelled 
upward or forward carried by the shock and reaction front of the explo- 
sion at velocities of up to several thousand feet per second. Velocities 
of this magnitude were verified through analysis of high-speed movie 




Time, sec. 





327 


£ilm coverage of liquid rocket propellant explosions both in the ADL 
experiments^ and the PYRD experiment where trajectories of large 
pieces of metil can be traced, giving their motion, velocities, etc. 

The engine thus moving forward penetrates the lower bulkhead of 
the S-IVB LQX tank then the upper one which is cannon with the LH 

2 

in the S-IVB. 

Failure of the LQX Tank (Hole or Shattering) : 

The impact of the thrust cone and engine assembly upon the 
LQX tank lower and upper bulkheads causes penetration of the thin 
metal walls. The minimum effect will be penetration by the thrust cane 
and engine, which will become mangled in the process. The maximum 
effect will be the shattering of the bulkheads. Both opinions have 
been advanced with, however, no experimental support of how the ma- 
terial of the S-IVB stage will react to such severe impact conditions. 

The amount of mixing of LOX and ll^ will depend on the size of 
the hole produced and the subsequent amount of motor penetration into 
the Uk tank. Again experimental determination of the material behavior 
could narrow the range of results predicted. 

Secondary Mixing in the S-IVB Stage : 

Depending upon the size of the primary explosion and the 
time at which it occurs, partial oi complete penetration of the thrust 
cone and engine into the Ili. tank of the S-IVB can be expected. This 
phenomenon, producing violent mixing, is referred to in our work as 
the EXPLOSIVE MIXING MODE. This mixing mode produced by an explosion, 
ramming the propellants and oxidizers together with extreme velocities 
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and violence* was simulated in our modeling studies and showed that due 
to the extreme mixing rates very high explosive yields can be expected. 

None of the other experiments with liquid rocket propellants such 
as the ADL experiments, Project FYRO, and the work of Mr. Fletcher 
covered EXPLOSIVE MIXING MODES and therefore cannot be used to make 
predictions about the explosive yields obtainable from such a mode 
of failure. Our work of simulation through model studies covered this 
mode as well as the others and indicated the much higher explosive yields 
obtainable. 

Even though the other works cited above did not cover the EXPLO- 
SIVE MIXING MODE, they indicated that the higher the energy supplied 
for the bringing together of the propellants such as dropping, firing 
into concrete walls, etc. , the higher the yield obtainable . The 
EXPLOSIVE MIXING MODE is one of the extreme cases of energy supply for 
the mixing process 

Figure 133 presents the mixing function for the secondary ex- 
plosion produced by a one tenth of one percent primary explosion. 

Ignition of Secondary Mix and Resulting Explosive Yield : 

Previous work has shown that the mixing energy affects the 
delay time before ignition less than the amount of fuel and oxidizer 
mixed. Thus, the greater the amount of energy supplied, the more 
mixing which will be accomplished before ignition occurs, and thus, 
the higher the explosive yield obtained . 

In the EXPLOSIVE MIXING MODE the size of the primary explosion, 
up to a value where maximum secondary mixing can be obtained, is 
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ng Function for the S-IVB LH 2 Tank Using a 1/50 Scale Mod?!. 
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important. If the primary explosion is larger than this value it 
consraes part of the LQX which otherwise would be available for the 
secondary explosion, resulting in a decrease in the total yield. 

Figure 134 presents the theoretical maximum explosive yield of 
the secondary explosion based upon the primary explosive yield. It 
can be seal that for very low primary explosive yields the secondary 
yields produced are also low. With a yield increase of the primary 
explosion the secondary maximum yield increases and then decreases 
again depending upon the delay time involved before the primary ex- 
plosion. Finally, if the primary yield is rather large the secondary 
explosion is smaller since not much oxidizer was left for this reaction. 

Figure 135 presents the experimentally predicted primary yield 
as a function of ignition time after destruct initiation. This figure 
also presents the secondary explosive yield resulting from an approxi- 
mate 0.11 primary explosion. 

Had the primary explosive yield increased (in the modeling studies) 
the resulting secondary explosion yield curve would have been shifted 
with the max imm value moving up and to the left. 

Results of the above analysis and phenomena following a destruct 
initiation, which can produce the S-U/S-IVB interaction, are sunmarized 
in Figure 136. It presents theoretical primary and secondary explo- 
sive yields and experimental primary yield predictions. It also pre- 
sents experimental secondary explosive yield predictions based upon an 
approximate 0.1% primary explosive yield. A higher primary explosive 
yield shifts the secondary explosive yield curve up and to the left. 
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Figure 134 Theoretical Prediction of Maximum Yield as a Function of Time After Destruct. 
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Figure 135A Experimental Prediction of Maximum Yield as a Function of Time 





Figure 13SB Expert aental Prediction of Maxi mm Yield as a Function of Tiae 




33S 



(A *1VS) *>f ‘PISW SAtS0ldX3 


Figure 136 Comparison of Theoretical and Experimental Maximum Yields 
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There may be same desire to establish the exact "likelihood** of 
the events described above. It does not seem too important to estab- 
lish this. Theoretical studies and the modeling, both two and three 
dimensional, simulating the primary explosions, show that these events 
will happen. The initiation of the destruct s)-stem may well be necessary, 
causing a sequence of events such as described above. The reported 
explosion on the pad of the Russian 10-14 million pound thrust super 
rocket during the summer of 1969 is a case in point. 

In the case at hand the problem of the S-II/S-IVB interstage 
failure is self * tposed through the design of the destruct system. The 
above problem can be eliminated by removing the ring charge from the 
S-IVB LOX tank resulting in a simpler, better, and safer Saturn V 
destruct system. 


Recommended Action for Preventing the S-IVB/S-II Interstage Failure 
The detailed analysis of the destruct system reported else- 
where, is believed to be able to function more or less as intended 
except for the S-II/S-IVB interstage failure. 

The sequence of events as analyzed and discussed in this paper are 
triggered by the setting off of the destruct charge on the S-IVB LOX 
tank. For this reason it seems that, from a technical point of view, 
the removal of the ring charge would prevent the occurrence of this 
sequence of events. 

Everyone with whom this problem was discussed, at Huntsville and 
Cape Kennedy, agrees that we would have a better destruct system with 
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that charge removed. 

If, however, positive destruct of the LOX tank is required, pan- 
cake charges could possibly be used to accomplish this vd chout let* ng 
the S-IVB and S-II interact as discussed. 


Conclusions 

It is believed that upon the foregoing work the following 

conclusions can be made: 

1. The methods previously developed by Dr. Farber's group regarding 
liquid rocket propellant explosions seemed well -suited for the 
detailed analysis of the Saturn V destruct system. 

2. A considerable amount of uncertainty still exists regarding the 
actual effect of the destruct system upon the Saturn V tankage 
configuration. Our reconmendation of several years a™o concerning 
experimental verification of the various opinions is still valid. 

3. An S-IVB/S-II interaction caused by the initiation of the destruct 
system can happen, and in our opinion, based upon theoretical 

and experimental modeling studies, is likely to happen. 

4. This S-IVB/S-II interaction can produce a primary explosion of the 
size predicted by other works, and in turn a much larger secondary 
explosion produced by the explosive •"ixing mode not previously 
investigated . 

5. The sequence of events described in this paper can be prevented by 
removing the ring charge on the S-IVB LOX tank, preventing the 
occurrence of a high yield explosion in this manner. 
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6. If positive destruct of the S-IVB LQX tank is required this can 
possibly be done by using pancake charges, thus preventing the 
above sequence of events which lead to a high secondary yield 
explosion. 


S-IVB Pancake Charge Destruct Analysis 

Because of the possible difficulties with the destruct 
system as it was designed and intended to be used, opening a hole in 
the bottom of the LC^ tank of the third stage, which in turn can 
produce a primary explosion leading to a larger secondary explosion 
in a critical region, alternative destruct methods were investigated. 

This seemed to be necessary if positive destruct of the tanks 
was to remain a definite requirement. 

For this reason the possibility of rupturing the tmks of the 
third stage, especially the IX^ tank, was studied. This tank 'oes not 
have actual contact with the outside skin and thus cannot be ripped 
open like the other tanks for propellant dispersal. This could be 
done by the proper design and placement of pancake charges. 

The pancake charges were simulated in both two and three dimen- 
sional models, see Figures 137, 138, by several methods. 

1. Pressurized air, blowing one propellant into the other. 

2. An air rifle which produces a blast opening the tanks and 
allowing the components to mix. 

3. A 22 caliber blank which acts similarly to the air gun but 
has a sharper and more powerful blast. 
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4. The pressurization of the tanks so when the opening are 
made manually one liquid is propelled into the other by the 
pressure in the tank. 

5. Pressurized, unpressurized and vented systems were analyzed 
and compared. 

Again depending upon the design of the charge different results 
can be obtained. If the shaped pancake charge is to basically blow 
a hole into the walls of both the IHg tank and the LC^ tank and both 
of them are pressurized, relatively little mixing is produced and with 
the liquids streaming out of vhe hole an external Hj fire may be 
produced. 

If, however, the shaped charge blows one constituent into the other 
backed by the pressure in the tank, then considerable mixing can occur. 

It might be observed that the size of the pancake charge in the 
modeling was relatively larger than would be used in the comparable 
prototype since true scaling would not have produced much effect on 
the models. Hie models used here weTe again 1/50 scale. 

The actual explosive yield even for th Tst case is only a 
small fraction of the maximum of the total theoretical maximum for the 
Saturn V. 

Hie position, however, may be critical and it may be better to 
produce a small explosion, or not to destruct the tank at all, if 
this is admissible for actual operation. 

Figure 139 presents the progression of the mixing process by show- 
ing a number of frames of the high-speed movie film. 
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An analysis of the records, frame by frame, results in a mixing 
function and expected yield as shown in Figure 140, and it can be 
seen that by this method of failure the mixing and thus the explosive 
yield will always be a small fraction of the theoretical possible for 
the Saturn V. Thus in the over-all picture the pancake charge destruct 
does a nice job in producing low yield but the location near the pay- 
load, nuclear packages, or the personnel capsule may be cause for 
some concern. 

On the last Figure, 141, in this section several typical yield 
functions are presented mostly to indicate the various possibilities 
of modeling. 

The fastest rising curve is the result of a glass model which was 
partly shattered by a 22 caliber blank which was used to simulate the 
pancake charge. Because of the shattering and the large contact area 
the mixing progressed rather rapidly in this case. 

The next curve over is essentially the same experiment repeated 
but with a plastic model which had a predrilled hole which was covered 
by tape and the tape was blown in with a 22 blank. Because of the 
smaller hole and therefore less contact surface the mixing progressed 
more slowly. 

The next curve presented gives the results when a plug is removed 
in the plastic model to open up the holes which otherwise would be 
produced by the pancake charge. Since the action is much less violent, 
the mixing which occurs is much slower. 

When an air gun is used instead of the 22 blank or the plug, the 




Prototype Tine 





Figure 141 Comparison of Various Pancake Charge Produced Failures 
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tape over the hole absorbed much of the energy and the mixing is 
still slower. It does not seem as good a simulation to the pancake 
charge effect. 

The last curve presented is that of the LC^ tank pressurized and 
not vented or opened up by the pancake charge. This makes it difficult 

to flow into the LO^ tank except for an ini* 
tial surge and the mixing never reached the completion it did in the 
other simulations . 

In all cases the mixing and thus the resulting explosive yield 
was a small fraction of the total potential of the Saturn V. 

Since this type of destruct system was not seriously considered, 
the work was not carried out as completely as in the others. It is 
of interest, however, since the type of failure and the subsequent 
mixing phenomena are different. 


to have the simulated 


This work is presented late in the report since it represents an 
application of the theories and methods developed under this 
contract. It was actually carried out in 1968. 
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Part IV Space Shuttle Analyses (Early Configurations) 


Introduction 

Since methods had been developed earlier that proved to be 
reliable for predicting the explosive yield, it was decided to apply 
the same procedure of analysis to the various configurations of the 
Space Shuttle and its most probable modes of failure. This information 
then became available parallel to the development of the shuttle rather 
than after the designs and configurations had become hardened, as proved 
to be the case with the Saturn V. 

It is much easier to use infoimation,if available, for guidance and 
decision making than to use that same information for making changes 
later. 

Hie early configurations were of the all liquid type with a booster 
and two different types of orbiters, the Delta Wing Orbiter and the 
200 Mile Orbiter. 

Again, 50th scale models were constructed from wire frames and 
plastic to simulate the tankage configurations and then they were mounted 
on the pad and in flight configurations so that with high speed pho- 
tography various modes of failure can be studies. See Figs. 142-145. 

With these models available the following modes of failui_ ere 
studied and the mixing functions obtained: 
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Figure 143 Space Shuttle Delta Wing Orbiter 
Bulkhead Type Failure 
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A. Bulkhead Type Failures 
Booster 

Delta Wing Orbiter 
200 Mile Orbiter 

B. Toppling on Pad 
Booster 

Delta Wing Orbiter and 2C0 Mi Orbiter 
Drop Off $ Fall Over 

Booster Hus Delta Wing Orbiter 

Booster Plus 200 Mi Orbiter 

c. C rash landing 
Booster 

Delta Wing Orbiter 
200 Mi Orbiter 


Other configurations and a number of different modes of failure 
could have been used, but the above were chosen as the most probable 
ones both iy our group and an independent study made by the Aero-Space 
Corporation . 

The results of the studios are shown below aid the discussion of 
the analyses, which are essentially the same as those used for tiie 
Sr'arrn V, are w>t repeated. 

The crash landing dynamic tests were simulated by hanging the 
r r 'ce shuttle models on chains and letting them swing down irum dif- 
ferent height'- r simulate different landing speeds, and letting tner 
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fly into the ground at different angles. A dashpot constructed from a 
55 gallon oil dnra nude it possible to stop the models in distances and 
tines realistic for the modeling of actual cradles and crash landings. 

Again, High-speed photography was used in recording the events, 
and the analysis of the mixing frame by frame produced the results for 
these different configurations and inodes of failure. 

Since .he actual configuration of the Space Shuttle now has taken 
on a Solid-Liquid combination, the model ii^ results discussed here 
have value mostly in showing how the modeling and analyses procedures 
can be applied to almost any real system once the configurations are 
known and the most likelv modes of failure selected. 

A. Bulkhead Type Failures 

These types of failures are probably the most likely 
ones and are investigated in this section. 

Booster 

On Pad Failure - He or RP Tank Explosion Rupturing 

Main LO^ and ill. Tanks : 

The booster model was in vertical or launch position and holes 
were precut into the tanks and then covered by duct tape. 

At time zero the duct tape was ripped of remotely and the simulated 
liquid propellants were allowed to mix while high speed cameras recorded 
the processes. 

Tn deciding upon the damage done to the propellant tanks and the 


hole sizes, it was assumed that the He tank explosion progresses 
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spherically outward and the intersection of it with the tanks produces 
holes in the respective tanks allowing the fuel and oxidizer to flow 
out and nix. 

The frame by frame analysis, translated into prototype tiro as 
explained and executed earlier, gives the mixing function for the 
Space Shuttle booster of this configuration. 

Five different failure modes are presented here with their re- 
spective results: 


1. He Tank Explosion, Snail Goncentric Holes: 

A He tank explosion cuts holes into the LO^ and the LH^ 
tanks one concentrically above the other and both of than 2.6 inches 
in diameter. It was assumed that the He tank explosion produces a 
spherical shock wave which cuts the holes into the main propellant tanks. 

The 1X>2 exiting from the top hole falls into the hole in the 
tank and mixing progresses. Translated into prototype time the mixing 
function is obtained for this case and this mode of failure. Fig. 146 



2. He Tank Explosion, Small Off-Set Holes: 

This case is similar to the first discussed above except 
that the two holes in the main tanks are off-set by one radius. 

When the LO 2 stream falls into the hole in the IJ^ tank off center, 
i xlows same of the 11^ to splash out, giving more room o v « n ullage 
space for the LO^ to enter. More mixing is obtained than was possible 
in the first case. Again with the same frame by frame analysis the 
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mixing function is obtained as presented in Fig. 147 


3. He Tank Explosion, Large Holes 

Case 3 is similar to case 1 except that it is assumed 
that the explosion of the He tank is more powerful and produces con- 
siderably larger holes in the main propellant tanks than were produced 
before. 

The simulated hole size in this case was 4.5 inches in diameter. 

While in the first two cases most of the mixing occurs inside the 
tanks, in this case a considerable portion of the total mixing occurs 
external to the tanks where the two liquids from both the L0 2 and IH 2 
tanks came together. 

The results from this mode of failure and configuration are shown 
in Fig. 148 


4. Kidney Shaped J P Tank Explosion (No Son) 

The next case analyzed was that of the kidney shaped J P 
tank exploding and again cutting holes into the main propellant tanks 
of the Booster. 

For this particular case it was assumed that the skin of the Booster 
was damaged or removed allowing the propellants to flow to the ground. 

The propellants were mixing as the streams were flowing together before 
they hit the ground and further mixing occurred after impact and puddling 
on the ground. 

Again making the frame by frame analysis and determining the mixing 
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in the impacting streams and then on the ground, the mixing function 
for the prototype Space Shuttle Booster is obtained. Fig. 149 

5. Kidney Shaped J P Tank Explosion (With Skin) 

The last case simulated and discussed here on the Space 
Shuttle Booster is similar to case 4 except that the assumption was 
made that the skin of the Booster was left intact by the explosion and 
the L 02 and are confined by the vehicle. 

Due to confinement and less chance for the propellants to spread 
out and disperse , the mixing produced is considerably higher and the 
mixing function again for the prototype is shown in Fig. 150 


Delta Wing Orbiter 

Two cases are taken up here. They are both on-pad failures due 
to He tank explosions. In the first case the He tank explosion blows 
holes into the main propellant tanks allowing the fuel and oxidizer to 
mix. In tne second case the He tank explosion produces a seam failure 
which opens a rough larger hole and therefore the mixing is much greater. 

The results from the analysis are again presented after frame by 
frame anal is and prototype time correlation. 

Case i Holes produced in the tanks give the results shown in 
Fig. 151 . 

Case 2. A seam failure gives the results shown in Fig. 152 
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Fig. 153 -» Tank Explosion in 200 Mile Orbiter 
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200 Mile Orbiter 

The case analyzed for the 200 mile orbiter is a He tank explosion 
which opens holes in both the LC^ and the LH 2 tanks with subsequent 
mixing. 

Again this failure occurs on the pad with the orbiter ii: vertical 
position. 

With the mixing occurring both inside and external to the tanks , 
rather high values are obtained at the end of t. i^ing process. 
Naturally, as was pointed out earlier, the "critical mass" will be 
reached quickly, not letting the mixing process ^rog'ess very far, re- 
sulting in a low explosive yield. The mixing fun: tion, however, since 
not terminated in the modeling, can proceed undisturbed right to the end. 

The results from this analysis are shown in Fig. 153. 


B. Toppling on the Pad 

In this section cases are taken which may occur if the 
booster or the orbiter or the combined Space Shuttle configuration 
should topple over on the pad. 

Again, these modes of failure have been considered a ,n ong the more 
likely ones and for this reason they have been modeled. 

Boo ster 

The toppling of the boos alone has been divided into two cases, 
one where, due to the toppling, the seams fail in ‘an. ’id the tanks 
split open, and two where the bulkheads fail or thi 'ndi c' th t; 


come open. 



365 


1. Seaa Failure 

In this simulation the booster was allowed to tip over 
when at the same tine the seams were opened by rip tape procedure. All 
this was done on a crossha tched platform so that the high speed records 
of th >nenomena could be analyzed accurately and carefully. 

The results from this series of experiments are shown in Fig. 154 

2. Tank Bulkhead End Failure 

The second case was essentially the same as the first but 
with different damage resulting to the tanks, namely the heads were 
allowed to come off the tanks. 

This failure allowed the liquids to splash forward and a certain 
amount of mixing occurred. Fig. 155 gives the results from the de- 
tailed film analysis. 

Delta Wing Orbiter 

The Delta Wing Orbiter in this series of experiments was allowed 
to topple over and at the same time the seams were split open simulated 
by rip- tapes. 

The resulting splash patterns as recorded by high speed photography 
were analyzed and resulted in the mixing function presented in Fig. 156 

200 Mile Orbiter 

The experiments carried out with the 200 mile orbiter were iden- 
tical to those of the delta wing orbiter, but due to the different 
internal configuration of the tankage, the mixing patterns and the 
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resulting mixing function look ! different. Fig. 157. 

Booster and Delta Wing Orbiter 

In these series of experiments the combination of the booster 
and the delta wing orbiter were allowed to topple over on the pad 
and while they fell the tanks were opened with rip-tapes simulating 
The failure of the tankage. The booster and orbiter separated upon 
impact, leaving an ever increasing splash and mixing area which was 
faithfully recorded by high speed photography and later analyzed frame 
by frame to obtain the mixing function. Fig. 158. 

Booster and 200 Mile Orbiter 

These experiments were carried out very much like the ones for 
the booster and delta wing orbiter and the results were obtained in 
exactly the same manner. 

Again the mixing function is obtained which represents the results 
of this mode of failure and this combined configuration. Fig. 159. 


C. Crash Landing 

In this, the last, series of experiments on Space 
Shuttle Configurations under consideration at the time, the experiments 
carried our included dynamic systems. Investigated were primarily the 
crash landings of the booster, and each of the orbiters, separately. 

To simulate these conditions, a physical system was set up to 
allow the booster or orbiters to be suspended by thin chains. In this 
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maimer the models could swing from various heights simulating different 
landing velocities and could allow, by varying the length of the chains, 
different landing or crash angles. 

To prevent bouncing of the models and to stop them, simulating 
the actual cases, a large dashpot made from a 55 gallon drum was at- 
tached to the models. This stopped them at different and predetermined 
rates of deceleration. Rip- tapes were also attached to the models and 
they were automatically removed as the model simulated the eras! landing. 

Again, the complete procedures and series of events were documented 
on high-speed films which were later used for analysis and the deter- 
mination of the mixing functions. 

More trials were needed for this procedure, until the right velocity 
simulation and dash-pot deceleration values were obtained. After this 
method was perfected the results obtained were reproducible. 

Booster 

The booster model with rip- tapes and the braking or dash pot 
system attached was filled with the simulated propellants. Then the 
model was allowed to swing down from a predetermined height and to 
impact at a shallow angle upon the landing surface. 

As it landed, the rip tapes were pulled, opening the tankage, and 
the braking system stopped the forward motion rather sharply. 

The resulting splash patterns were recorded with high-speed pho- 
tography and then analyzed frame by frame by tracing and then planime- 
tering the various areas, (total, splash, and mixing). 

From the analysis the following results are obtained. Fig. 160 
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Delta Wing Orbiter 

The crash landing simulation was repeated for the delta wing 
orbiter, and except for the configuration and the braking system and 
rip- tape attachment, the experiments were carried out the same as for 
the booster. 

The analysis of the high-speed photographic records again allowed 
the determination of the mixing function for this type of experimental 
simulation. 

The results are shown in Fig. 161. 

200 Mile Orbiter 

The last in the series of experiments of the early Space Shuttle 
Configurations was the crash landings of the 200 mile orbiter. 

The experiments were carried out in the same manner as those with 
the booster and the delta wing orbiter. The 200 mile orbiter after 
being filled with simulated fuel and oxidizer was allowed to swing 
down and crash land on the simulated landing platform. The resulting 
splash patterns, again recorded with high-speed photography, were then 
analyzed giving the results reported in Fig. 162. 
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Final Closure 

The foregoing discussions and presentation of results cover 
a large amount of both theoretical and experimental work. 

A number of new concepts are presented, and verified experimentally , 
to give several new approaches for the evaluation of liquid rocket 
propellant phenomena, and to allow their analyses with regard to mixing 
characteristics, basic behavior and explosive yield predictions. 

For the first time, data were taken inside exploding missiles, 
measuring mixing characteristics, location of ignition points, shock 
front and reaction front behavior and velocities, etc. 

Three methods for determining explosive yield were developed and 
verified in the laboratory and field and then used for analyses and 
prediction. They are the 

1 . Mathematical Model 

2. Seven Chart Approach 

3. Critical Mass Analysis 

These methods, although quite different in approach and independent 
from each other, give the same results when used for analysis of a 
particular problem. 

It is believed that this work as presented here is a valuable 
contribution to Space Technology, and Liquid Fuel Technology in general , 
and a significant step in making Space travel and Industrial processes 
safer. 
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Computer Programs for the Solution of Equations ( 156 ) end ( 172 ) 



EQUATION (2.<b) With Fences . 


,_./J 0 '» c 1 lAeia,3»6).»**._Sf. STEWART. _ . ... 394 

C DCSTWUCT SPILL PATTERN (CIRCULAR HYDRAULIC JUMP) 

C . RO » YO.f 0 .=_ RADI US. DEPTH »FROUDE_.NO_A.T_BEGI NNI NC ..OF„JUMP 

C Rl.YI.FI = RAD I US. DEPTH .FROUDE NO AT END OF JUMP 

C R » Y » F _ = RADIUS. DEPTH. FROUDE NO. _DQWNST.BE AM _0f _J J MP 

C RR1=R/R 1 ,RRO=R 1/R0.DR0=Y1 /YO.DRI *Y/Y 1 


1 

2 

3 

— 

D IMENSION F0( 10) .RRO( 1 5) , RRI (10) 
REAC( 5. 10 ><F0( J) . J=1 .10) 

READ(£, 1 1 >(RRO< J) • J=1 .15) 

\ 

4 


READ! S. 10) (RRI C J) • J*1 .10) 


5 


DO 6 1*1.10 


6 


WR 1TE( 6.12) 


7 


WR ITE( 6, 13) 


a 


WRITE(6.14 ) 


9 


WR I TE( 6. IS) 


10 


rN=6o*F0( I )**2 


i 1 


DO 6 J=1 , 1 5 


12 


A2=-(RR0( J )- lo )/( 2o*RR0( J) + l» ) 


13 


A1=-(RR0( J ) +FN+2© )/(2 c *RRO( J) + i© ) 


14 


AO=FN/(RRO( J)*( 2o*RR0( J) + lo ) ) 


15 


Xl=FO( I ) 


16 

1 

X- / 1 


17 


EQN=X**3*A2*X#*2<-Al*X+A0 


ia 


DER IV*3«*X**2+2o*A2*X+Al 


19 


2= EON 


20 


ODER IV 


21 


X 1* X— ( 2 /D ) 


22 


I F ( ABS( X-Xl )— 0*0005)2 *2 .1 


23 

2 

ROOT 1 = XI 


}24 


fcl*A 2+ROOT 1 


F 25 


AC=- AO/ROO T 1 


26 


D*B**2-Ao*AC 


27 


IF(D)4. 3. 3 


28 

3 

ROOT2=(-B+D**Oo5)/2o 


29, 


ROOT 3*(— B— D##0o5) /2© 


30 


GO TO 5 


31 

A 

ROO T 2=0o 


32 


ROO T3=0o 


23 


WRITE( 6, 16 )F0C I ) »RR0( J) .R00T1 .ROOT2.ROOT3 


34 


GO TO 6 


35 

c 

WR I TE( C.17)F0( I) .RR0( J) .ROOT1 .ROOT 2 • R00T3 


36 ~ 

6 

CONTINUE 


37 


DO 112 J*l»10 

• . 

38 


Fl=l ©05 


39 


WRITE16. 18) 


40 


WRITE! 6. 19) 


41 


WRITE! 6.14) 


4 2 


WR ITE( 6 . 20 ) 


4 3_ 


DO 112 1=1.20 


44 


Fl = Fl-0o0S 


45 


FN=6o#Fl**2 



46 A~2*-TkRl( J )-1 o )/( 2© *RR 1 ( JKl© ) 

4 ? A 1=-(RR 1 ( J ) +FN +2 • ) /( 2c »RR 1 ( JHlt ) 

40 ' A0=FN/(RR1( J) *(2o *RRI( J)-»-U )•) 


A <i . _ X 1-F 1 . 

50 , 7 X*X1 

51 EON = X**3 +A24X *_+.2+A_l * X+A O_ _ i 

V UTk I V-3©*X**2+2, *A2«X*A1 , 

• • \ 
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56 Xl~X-(//L>) 

56 JF(AKS( X-Xl )- 0o00ii5) 8*8.7 

£ 7 g UUO T 1= XI 

58 ~ B= A2+WOOT 1 

50 AC=— A 0 /ROO T 1 

I 6 C 0=U**2-4«*AC 

61 1F!D)11C.9.9 

62 9 M00T2=!-D+D4*0o5)/2» 

6 3 ROD T3= ( -B- D **0o_5 » /2o 

64 GO TO 111 

6 5 110 ROOT2=Co 

66 ROOT 3~0o 

67 WR ITt ( 6*21 » FI *RR1 C J) .ROOT1 . ROOT2.ROOT3 

68 GO TO 112 

69 1 11 fcR ITS! 6.22 1 F 1 »RR1 ( J) . R OO T! . KOOT2 . R OOT3 

70 112 CONTINUE 

7 1 .10 FORMAT! ICF 60 I ) 

72 11 FORMAT { 15F£o2> 

73 12 FORM/ { ✓// * 1 ‘'•X.23HCI RCULAR HYDRAULIC JUMP ) 

74 13 FORMAT { ✓ » 6 X * 41 H~ ROUDE~ RADIUS DEPTH RATIO CY 1/YO 1 > 

75 14 FORMAT! 6X.36HNUM3CR RATIO THREE ROOTS ) 

76 1 5 FORMAT !7X , 45HIF0 ) ( R l /RO) ROOf l ' ROOT 2 ROOT 3 V 

77 17 FORMAT ( SX* F 60 1 »3X»F 5© 2*4X*F9o4»lX *F9» 4 » 1 X »F 9 c 4 ) 

78 16 FORMAT! 5X.F6& 1 »3X.F5©2 »4X.F9o'4%9H1 MIGINARY 12 .9HIMAGINARY 12 1 

79 18 FORMAT! /// . 1 2X .28 HP0 WNSTRE AM OF H YDRAULIC JUMP > 

60 19 FORMAT! ✓,6X,4CHFROUDERAOi US DEPT HR AT IO !Y/Y1) ) 

e 1 20 FOR MAT! 7X *45H! F I) ! R/ Rl ) ROOT 1 RO O T 2 ROOT 3 ) 

82 21 FORMAT! 5X.F6o2.3X.F6ol . 4X »F 7 o 4 ,1 X.9HI MAGI NARY 12 .9 H IMAGINARY - 1*2 1 

83 22 FORMAT! 5X.F6p2 .3X .F6o 1 »4X «F7 q 4»1X» F9i>4.1X»F9o4 > 

84' " ’ STOP ~ 

|85 _ . END 




C If'CULAR 

HYDRAULIC. 

JUMP _ 

- * ----- 

froude 

NUMOLR 

crc ) 

K \ r> » L s 
RAT 10 
Clil/RO) 

DEPTH RATI O < Y 1 /YC) • 

THREE ROOTS 

ROOT 1 ROOT 2 ROOT 3 

\ 

20 C oO 
2 Of c 0 
2 0 v e 0 
2 l.C r 0 

loCC, 
1 o05 
hi 0 

Ic lb 

282*3430 
277*7756 
2 73*4229 
260*2666 

Oc 9999 ->283*3428 
Oc 9523 -278*7114 
0* 9089 -274,3005 
0*8695 -270*0923 

^ 

20 c «. u 

20 C | C 

1*20 

1 o 2 £ 

265*2979 
261 o 4980 

0*8333 —266* 0723 
0*7999 -262*2263 


200*0 
2 C "tO 

1*20 
1 o 35 * 

257o 6569 
254*3637 

0*7692 -258*5425 
0*7407 — 255e C098 


tGCoO 

2C*C*cO 

1 * AO 
lo45 

£51*0089 

247*7835 

0*7143 -251*6179 
Oo 6896 —240a 3577 


200*0 
20 C 1 0 

1*50 

2 *QC 

244*6794 

218*9406 

0*6666 -245*2209 
0*5000 -219,2405 


200 cO 

?CC *0 

2*50 

3*C0 

199*9267 

185*1421 

Oo 4000 —200* 0767 
0*3333 -185*1897 


200 oO 

4*00 

163*34 31 

0*2500 -163*2595 

* 





i ! 


CIRCULAR 

HYDRAULIC 

1 

JUMP 

i 

FHOUDE 

NUMHCR 

RADIOS * 
RAT IO 

DEPTH RATIO ( Y1 /YO ) 
THREE ROOTS 

• 

CFO) 

1 C 0*0 

CR1/R0) 

1*0C 

ROOT 1 
140*9222 

ROOT 2 ROOT 3 

1 o 0000 — 1 4 1 o 92 2 1 


lOOoU 

lOCoO 

1*05 

1*10 

138*6547 

136*4930 

Oc 9524 -139*5908 
0*9091 -137*37 08 


100 »0 
lOOoC 

lolS ~ 

1 o 20 

134c 42 95 
132c 4566 

0*8695 -135*2535 
0*8333 -133*2311 


100*0 

lOCoO 

1*25 

1*30 

1 30* 5683 
128*7585 

O 08 OOO -1 31 0 2959 
0*7692 — l29o 4 442 


lOOoO 
)0C*0 _ 

1 o 35 
1*40 

127*0219 

125*3539 

0*74 07 -127*6680 
0*7142 -125*9628 

1 

lOCoO 
1 00 oO 

1*45 

1 o 50 

123*7499 

122 s 2061 

Oo 66 96 — 1 24* 3241 
0*6666 -122*74 77 


ICO* o' 

irooO 

2*00 
2 c 5C 

109c 3976 
99*9265 

Oo 5000 — l C9o 6976 
0*4000 -100» 0784 

i 

ICOeO 
100*0 

3*00 

CC 

92o5620 

61*6955 

0*3333 — 92e 6095 

0*2500 -8U6122 

! 

) 







CIRCULAR 

HYDRAULIC 

JUMP 

i 


FROUDE 

RA01LS 

DEPTH 

RATIO CY1/Y0) 

NUMHER 

RAT IO 

THREE ROOTS 


CFO 

(R 1/RO ) 

ROOT 1 

ROOT 2 

ROOT 3 

5C*0 " 

1*00 

70*2124 

1*0000 

-71*2124 

_ 50*0 

lo05 

69* C 949 

0* 9523 

-70*0311 

5C*0 

1,10 

60*0290 

0* 9090 

- 68 * 9063 

SC*0 

lc 15 

67 o 0109 

0* 8695 

-67*8350 

50 *0 

" 1*20 ~ 

66*0372 60 8333 

-66*81 17 






SO o 0 
SC* c 0 
sc? 6 

S( o 0 

50©* 
SC cO 
SO o ** 
SC O 0 
SjT r 0 
5c 1 0 


FROUDG 

HAD I6S 

DEPTH 

RATIO CYl/VC) 

numfer 

RAT ID 

THREE ROOTS # 

<FO ) 

( R l /R 0 ) 

ROOT l 

ROOT 2 ROOT 3 

25c0 

loCC 

34© 8589 

lc OOOO -35© 85 09 

25 a 0 

loOS 

34© 31 68 

Oo 9523 -35o2533 

__25o 0 _ 

110 

33©_7991_ 

Oo 9090 —3.4 ©67 63. 

2boO 

1d15 

33o 3C41 

Oo 8694 —34© 1 2 80 

25 rO 

_ 1 a 20 

32o 0 30 1 _ 

0©_8331 ^33o 6044. 

25 oO 

lo 25 

32o 37 50 

0© 7997 -33ol041 

_25 o 0 

1 o 3 0 

31©9399„ 

0#7689_ .—32® 62 55. 

25 oO 

lu 35 

3 1« 52 1 1 

Oo 7404 — 32o 1670 

. 2 _S _oJ)_ 

I o AO 

3 ljt> l l es_ 

Oo 71 39 — 31 o 7272 

25 o 0 

1 o AS 

30c 7309 

Oo 6893 -31 o 3043 

_ 25 o 0 

1 3 SO ' 

30©3576„_ 

Oo 6663 -30o0?88_ 

25o0 

2* CC 

27o2486 

Oc 4996 —27© 54 82 

.. 25©0 

_ _ 2o SO 

_24©_9391__ 

Oe 3996_ -25c 08.87. 

25 o 0 

3o00 

23c 1370 

Oo 3329 —23© 1 842 

25c 0 

AoOO 

HRuhI 

Oc 2496 —20© 3874 


1-25 

1 048 ___ 

. _ Oo 7999 

—65© 8332 

1 c 2V 

64* 21 07 

Or 7691 

-64c 8965 

1 v 35 

6 3©.3526_ 

Oc.74 06. 

— 63 ? 99tt6_ 

1 r 4C 

62o 5© OO 

Oo 71 42 

-63* 1369 

1 .i 4 5 

_ 61 c 7349_ 

0© 6896 

-62© 3091. 

lc50 

60o 9712 

O 06666 

-61© 51 28 

2 o C C 

_ 54o62 84 

_ _ Or 4 999 _ 

-54© 9283 

2 50 

49o 9320 

Oc 3999 

-S0o0819 

3*00 

46o 2748 

Oo 3332 

— 46 0 32 23. 

4c OC 

40©8750 

Oo 2499 

—40© 7915 





CIRCULAR 

HYORAULIC 

JUMP 



3W 


Cl R CULARH YOR A UL IC JUMP 


F TOUOE 

RADIUS 

DEPTH 

RATIO (Y1/Y0) 

NUM BER 

RAT 10 

THREE ROOTS 


_(F0 ) 

CR l/ROi. . . 

R OOT 1 

ROOT 2 

ROOT 3 

20 oO 

lo OC 

27o7887 

lo 0000 

-2 8*7887 

_ 2C 0 0 _ _ 

loCS 

__.27.o36 19 _ 

.0o9523 

-2.8*2980.. 

20 oO 

1 c 1 C 

26© 9539 

Oo 9089 

-27© 8316 

20 oO 

lo IS 

26o 5636 

Oo 8693 

-27o3875 

20 oO 

1 o 20 

26© 1 898 

Oo 8330 

-26© 9639 

20 oC 

1 325 

25o 831 2 

..Oc .7996 

-2.6* 5593.. 

20 ©0 

1 0 30 

25© 4870 

Oc 7688 

—26© 1 7 

. 20o0 

1 0 35 

2 5 o.l 56 2__ 

Oo 740.2 

jraso.ROi e_ 

20 oO 

1 © 4 c 

24o 8380 

0© 71 38 

-25o 4465 

_20©0_ 

lo4S 

24 0 5316 

0o6891 

-25o 1053 

2C~oO 

lo5C 

24© 2363 

0© 6661 

—24© 7774 

20 0 0 

2oC0 

21o 7?46_ 

0o4993_ 

-22o0739_ 

20c 0 

2© 50 

19o 94 27 

Oo 39%l 

-2 0© 0920 

k( oO 

3e CC 

1 8©51 18 

Oc 3326 

-1 6© 5587 

2o . C 1 

4 i CO 

1 6© 3930 

Oc 2494 

-1 60 30 90 


C IRCULAR 

HYORAULIC 

JUMP 

— 


395 


FROUOE 

RAD l C S 

DEPTH 

RATI 0 ( Y1 / 

YO> 

NU^.ntH 

RAT 10 

THREE ROOTS 


(HO > 

Ol l/RC ) 

ROOT X 

ROOT 2 

ROOT 3 

15w0 

loOC 

20o 71 91 

lo GCCO 

—.2 1 o 71 91 

X5oC 

X o 05 

20 o 40 77 

Oe 9522 

—2 X o 3437 

l£eO 

lolC 

20 0 ltS7 

Oo 9007 

-2 Oo 9872 

15o0 

iclS 

1 9o8243 

Oo 8690 

••2 Oo 6479 

. JSc«V. 

1 c 20 

19o.55.06_. 

0c.8327.__ 

-20o_3245 

1 Sc 0 

1 e 25 

X 9o 2Q79 

Oe 7993 

-20 o C157 

IS.,0 

1 c 3C 

1 9o 0355 

Oo 7664 

-19o72 05 

15oO 

1 a 35 

18* 792 7 

Co 7399 

>1 9» 4380 

l S v 0 

I c 40 

X6o5590 

Oe 71 34 

-1 9 0 1 671 

1 So c 

lc 4S 

1 6o 3339 

Oo 6887 ’ 

-18e9072 

lJSoO 

j £ 5C 

X eoll 68 

Oo 6656 _ 

-1 80.6575 

~ I5«r 

2o00 

X 6o 3027 

Oo 4 988 

-16o6015 

15*0 

2c 50 .. 

X4 9 94 0B 

.. _ Oo 3988 

-1 5o 0975 

X So O 

3o00 

1 3o 6693 

Oo 332 1 

-13 0 9357 

X So 0 

4oOC 

1 2o 3 1 76 

Oo 24 89 

-1 2o 2333 



CIRCULAR 

HYDRAULIC 

JUMP 


FROUDE 

RADIUS 

DEPTH 

RATIO (Yl/YO) 

NUMEER 

RATIO 

THREE ROOTS 


CFO) 

( R 1 /R 0 ) 

ROOT 1 

ROOT 2 

ROOT 3 

XC oO 

1 oOO 

1 3 o 65 1 0 

io oboe 

— 1 4 © 651 0 

lOoO 

lo 0 5 

1 3 o 4555 

Oo 951 9 

— 1 4 © 39 1 3 

10 eO 

lo 10 

1 3 o 2679 

Oo 9083 

—1 4 o 1 449 

lCcO 

lo 15 

1 Jo 0877 

Oo 6664 

— 1 3 o 91 06 

1— 

O 

O 

O 

1 0 20 

1 2 o 9 1 45 

Oo 831 9 ' 

— 1 3 o 6875 

*-* 

o° 

0 

X 0 25 

1 2 o 7 4 79 

Oe 7984 

— 1 3 © 4749 

10 oO 

1 0 30 

1 2 o 58 76 

Oo 7674 

— 1 3 o 271 7 

10 oO . 

1 r> 3 S 

1 2 o 4 331 

Oo 7388 

— 1 3 o 0773 

10 0 0 

O 

O 

0*4 

1 2 o 2642 

Oe 71 22 

— 1 2 o 891 1 

lOoO 

1 o 45 

1 2 o 1 405 

Oo 6875 _ 

-l_ 2 o 7126 

10 oft 

1 0 50 

12 oC.O 16 

Oo 6644 

- 1 2 o 54 1 2 

lOoO 

2 oOC 

1 Oo 8365 

Oo 4973 

“lie 1338 

1 C oO 

2 o SO 

9 o 9612 

Oo 3973 

-lOo 1084 

10 oO 

3 oCO 

9 o 2735 

Oo 3306 

- 9 o 31 84 

lOoO 

4 o 00 

80 2500 

Oo 24 74 

- 80 1641 



CIRCULAR 

HYDRAULIC 

JUMP 


FROUDE 

RADIUS 

DEPTH 

RATIO CYI'/YC) 

NUMOER 

RAT IO 

THREE ROOTS 


(FO) 

CR1/R0) 

ROOT I 

ROOT 2 

ROOT 3 

7 o 0 

XoOO 

9o 4 1 21 

lo 0000 

-1 Oo 41 21 

7cO 

1 o05 

9o 2869 

Oo 9514 

-1 Oc 2222 

7 c 0 

lo 10 

9o 1661 

Oe 90 74 

— 1 Oo 0423 

7 0 0 

lo 15 

9o 04 97 

Oo 8672 

. — 9* 87 X 4 

7 o 0 

lc 20 

8c 9373 

Oo 8304 

— 9o 7087 

7 © 0 

lc.25 

80 82 08 

Oc 7967 

“9o 554 C 

7r 0 

lc3C 

Or 7241 

Or 7655 

-9r4'*63 

1 

; 

0 

V 7.11 " 

—v 7: i* L* 

* x > 

• . >♦ 1 > 


Or 71 nt* 

-Do 1 2 DO 


U45 Co 430b Oo 6852 -9o(>004 


. > > 

7 o 0 
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?v 0 

7 C 0 
7^0 

1 *• SC 

2 c 00 

2? SO 

fio 3352 
7c 5648 
6 o 9775 

. Oc 662 C. 
0« 4945 
Oc.3.944. 

Oe 0761 . .. 396 

-.'0 8593 
-7» 1219 

7o 0 
__ . 7 o 0 

3c UC 

AcCO 

60 S 135 
So 81 98 

Oc 3279 
_£« 24.48 . 

-60 5557 
— 5o 73 1 3 




\ 

CIRCULAR HYDRAULIC 

JUMP 


FROUDE 

RADIUS _ 

..DEPTH.. 

RATIO (Y1/Y0) 

NUMLfc'R 

RAT IO 

THREE ROOTS 


(Ft ) 

(R 1/RC > 

ROOT 1 

ROOT 2 

ROOT 3 

4 o C 

loU(: 

5c 1 7 09 

lo 0000 

—60 1 7 89 

4oO 

1©C5 

5© 1265 

©0 94 94 

• 60 0598 

4 o 0 

lo 1 C 

5c 0748 

Co 9037 

—5© 9472 


icie 

5c 0 238 

Oe e 621 

-So 84 C5 

4 oO 

lc 2C 

4« 9738 

Oo 8243 

-So 7392 

4 oO 


4c 9248 

60 7896 

-So 6429 

4 c 0 

1 o30 

4o 8 76Q 

Oe 7577 

— 5©5512 

4c 0 

1 e 35 

4o 8300 

Oe 72 83 

—5© 4 63 T 

4 o 0 

lo4C 

4o 7843 

Oo 701 0 

— So 3801 

4o0._ 

1 o 4 5 

4o7398 

Oo 6758 

— 5o 3001 

4©0 

loSC 

4o 6963 

0o6522 

—5o 2235 

4 ©0 

2 o 00 _ 

4c 31 74 

Co 483? 

— 4o 6007 

4 ©0 

2 a SO 

4 o 0203 

Oo 3833 

— 4o 1 536 

4 o 0 

3o00 

3c 7816 

Oc 31 70 

— 3o 81 29 i 

AoO 

4©0G 

3o4 196 

0* 234 8 

— 3o 3210 1 



CIRCULAR 

HYDRAULIC 

JUMP 





FROUDE 

RADIUS 

DEPTH 

RATIO (Y1/Y0) 




NUM eER 

RATIO 

THREE ROOTS 





(FC) 

{ R l /R 0 ) 

ROOT 1 

ROOT 2 

ROOT 3 




2 oO 

loOC 

2o3723 

lo 0000 

— 3o 3723 




2 oO 

lo0£ 

2*3805 

0« 9379 

— 3o 3023 




2 ©0 

lo 10 

2e 3 839 

Oo 8837 

— 3o 2364 




2 0 0 

1 0 15 

2o 3 638 

Oo 8358 

— 3o 1 742 




2oG 

1 0 20 

2o 38 1 1 

0© 7930 

—3© 1 1 53 




2o0 

1 0 25 

2*3763 

0* 7545 

— 3o 0594 




2o0 

1 a 30 

2*3700 

Oe 71 97 

— 3o 0064 




2o0 

lc 35 

2o 3625 

Oo 6880 

— 2o 9560 




2 oO 

1 0 40 

2o 3541 

0» 6590 

— 2o 9079 




2c 0 

1 0 45 

2o 34 SO 

Oo 6323 

-2*8623 




_2 o_0_ _ 

1 0 50 

2o 3354 

Oo 6078 

— 2o 81 82 




2 Jo 

2oOC 

2*2292 

Oo 4366 

-2o 4658 




2e0 

2*50 

2*1265 

Oo 3390 

-2o 21 75 




2o0 

3o00 ” 

2o 041 1 

Oo 2757 

—2*0311 




2 oO 

4 oOO 

1*9018 

Oo 1 984 

—1 0 7668 





_ 00 WN STREAM Of HYDRAULIC JUMP (WITH S»OS fORCfcS) . 

FROUDjE RADIUS PEPTH R ATIO ( Y/Y U 

NUMBER RATIO THREE ROOTS 

(FI) (H/'RU ROOT 1 ROOT 2 ROOT 3 







397 


0o9S 

1 o 1 

- 1 o 85 29 

lo 1 01 4 

0* 702 8 

Ti 

0o<5C 

Ip 1 

Co 67C 7 

lo 1 503 

-1 17897 

J 

0oR5 

It 1 

0*6347 

lo 1236 

- 1 o 7270 


0 o 80 

lo 1 

0 o 596 1 

1*101 1 

“lo 664 9 


Or 7 t> 

Id 

Oo 552 5 

_ lo C823 

-1*6035 


Os 70 

lol 

0 o 50 7 5 

1 o 0667 

-1>54>9 


o« es 

lo 1 

0 o 4607 

1*0536 

-1 o 4833 

\ 

0 c fcO 

lol 

Oo 4 1 29 

1 o 0431 

-1*4240 


Or 5 5 

lo 1 

Oo 3645 

1*0343 

-1*3676 


Co 50 

lol 

Oo 31 63 

1*02 70 

-1*3120 


Oo 45 

lo 1 

Oo 26e 7 

lo 0209 

-1*2563 


Oo 4 1 - ' 

lot 

0 o 2224 

1* 0159 

-1*2071 


O o 3 5 

lo ! 

Oo 1781 

UOl 1 8 

-lol 586 


Oo 30 

lo 1 

Oo 1 366 

1*0064 

-1*1137 


f c 25 

lol 

Oo C9 87 

1 » C057 

-lo0732 


0 o 20 

lol 

Oo 0655 

1*0036 

— l • 037 8 


0» 15 

lol 

0*0379 

1*0020 _ 

—1 o 0087 


0«, 10 

lot 

0*0173 

lo C009 

-0*9659 


Oo 05 

lol 

0*0044 

1*0002 

-0* 9733 




OO WN STREAM 

OF HYDRAULIC JUMP 



FWOUOE 

RAO IL S 

DEPTH 

RATI 0 ( Y/Yl ) 


NUMBER 

RATIO 

THREE ROOTS 



(F 1 1 

Cft/Rl ) 

ROOT 1 

ROOT 2 

ROOT - 


1 o 00 

1*5 

1 o 3403 

Oo 4484 

—1 o 6638 


0o95 

loS 

1 o 3036 

Oo 4303 

-1*6089 


Oo 90 

US 

1 c 2690 

Oo 4 1 06 

-1*5546 


Co 65 

1 o 5 

- U5CC8 

1 o 2365 

C* 389 3 


Oc 60 

1 o 5 

- 1*4476 

lo 20 63 

0*3664 


Oe 75 

1 o 5 

- 1*3955 

lo 1 764 

Oo 34 2 0 


0 o 70 

1*5 

1*1529 

Oo 31 62 

-1*3441 


0 o 65 

1*5 

0*2891 

lo 1296 

-1*2937 


Oo 60 

1 o5 

0*261 0 

lc 1 085 

-1*2445 


Oo 55 

1 o 5 

0*2320 

lo 0897 

-1*1956 


Co 60 

1 o 5 

0*2025 

lo 072 9 

-1*1504 


Oo 4 5 

us 

0*1730 

1 o 0581 

-1*1061 


Oo 40 

1 • 5 

0* 1438 

1 o 0453 

-1*0641 


Oo 35 

.. . . _.lo £ 

Oo 1156 

1 o 0342 

— 1*0248 


Go 30 

1 o 5 

0*0888 

lo 0248 
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Cc?MAT<0/O (2.27.) . .. ... 

( 1 14 8 IE ,3 .6 ) . JOHN MILLER 
F.VCRALLIC JUMP IN RADIAL FLOW 

TO GENERATE A FAMILY OF CURVES OF FJ VS Yl/MJ FOR SEVERAL 

Y 1 = DEPTH OF FLOW AT THE END OF JUMP 

' RJ "s' RADIUS OF IMP l NG I NG JET 

FJ •= FROUDE NUM3ER OF. I MP1 NGI KG._JET 

VO = VELOCITY BEFORE THE JUMP 

VJ = VELOCITY OF IMPIN GIN G JET 

A = V 1 /RJ 

_S A_ =_ VO /V J 

initialization of jjata 

CATA b , C /- C © 6 7 2 * 0 o 2 4 5/ 


CO CC N = 


IFIN *EQo l) GO TO 1 

IF ( N »EQc_4) GO TOJ 

I F t N oEQ o 7) GO TO*l ~ - - - 

I F (N a EQ o 10) G O TO I 

WRITCCe.100) , 

JFO RM AT !✓//, T1 Z ,» VO/VJ* ,T24 » «V1/RJ« *T36 »«F J« ) 
WRITE! 6*101) . 

FORMAT ( • *• *_T 1 I , • * » T23 » * » »T35 > • 

GO TO 2 

WR ITE! C* 102 ) ’ 

FORMAT! • 1 • »T 12 . • VO/VJ* *T24 , • Yl/RJ* *T36 * »F J‘ > 
WR I TE ! 6 , 1 0 1 ) 


A = C*S 

00 10 M =1.14 

R=2 Co G*IA * * 2 oO )* E A 

S-( 2 2©C*! A**2o 0 ) *SA*C ) -!2» Q»A »C) -!3o0 *C/4« 0 
f = ( 8 oO * ( A ** 2 V 0 ) *( C* *2© 0 ) )V('4od*!C**2«0)/SA) 

= 4 o C*A*!C>i *2o 0 )/SA 

= CoC*! SA**2«0)/A 

= ! 20 05= SA*C/!A *»2oO) ) — 24p 0 ♦! SA**2oO) *C 

= 5»0*A4SA/C 
= 2C0o0 


/JL 


*A*S A ) 


K = l 

2C K = K + l 

d 1F!K-1C)50*5G*80 

9 50_R.ES = JR # JF J**3 ) )*!S» !F J* *! 2® 0*B) M + ! T* (F J* * (3 o 0 *0) ) )- 1U» ! F J*» 1 4oC 

4*3 ) ) )-! V*!FJ**(2oO+B) ) )+,! W*!F J**! 2o0*! 1 ®0+B> ) ) )+X 

cO CER = B*R*F J**!B-»l ©0)+ 2o Q» B*S* F J**!2oP *6-1 oO )+3® 0*B*T*F J *» !3® 0* E- 

ilikO )-4oG*D*U#F J**( 4®0*B~le0) -!2® (t+ 8) *V*FJ**!1 ©0 +B ) + 2©0 * ! 1 *0 )*W * 

SFJ»«! 1 » 0* 2 oC*B ) 

IF! ABS(RES‘)-0o"0bl ) 60*60*12 

l SING NEWTON’S M ETHOD CF NUMERICAL SOL UTIONS 

FJ g FJ- RES/DER 

CER =_ F IRST P ER I VAT I VE OF RE S 

GOTO 20 

80 WR 1 T E ( 6 « £ C ) SA * A »F J 

60 FORMA T ! • • , T 1~3 *F 4o 2 »T25 *F 3o 1 . Tl4 »'f6 0 2 ) 

A = A + C o5 

1C "CONT INUE ~ 

SA = SA + Co C2 
90 CONTINUE*' 

STOP 
EN C 
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